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INTRODUCTION=: 


The SIR POCKET GUIDE has been designed to provide quick referer 
to the commands and options available to the SIR* user. It also pr 
vides a pocket reference to commands and options in both the SPSS £ 
BMDP statistical systems. This guide assumes that the user is fe 
iliar with each of these systems and wishes to use this document or 
for a quick reminder on commands and options. This pocket guide 
not intended to be a users manual substitute for any of these system 
Throughout, this guide, items inclosed in [ ] are optional while { 
enclose lists of allowable specifications from which one may 
chosen. 
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CHARACTER SET 


A-Z, 0-9 may be used for naming variables, files and function: 
i +, ■ and E may be used for numerical representation. 

Blanks and commas serve as delimiters. Multiple blanks ai 
ignored except when in string constants enclosed in quotes. 

In COMPUTE commands: + add, - substract, * multiply, / dividi 
** exponentiation, = is equal to, ( ) grouping delimiters, ; 1 
separate multiple expressions appearing in a single C0MPU1 
command. 


In format control: /,*,X,T,I,F,E,A and picture clause ed: 
specifications (see the WRITE command for details). 

In logical expressions: LT less than, LE less than or equal tc 
GT greater than, GE greater than or equal to, NE not equal to, E 
equal to, with logical connectives AND, OR, NOT. 


*SIR, Scientific Information Retrieval; USERS MANUAL; Robinson, 
Anderson, Cohen and Gazdzik, SIR Inc. ,Box 1N0N, Evanston Ill., 
60204, 1979. 
































SIR - GENERAL RULES 


VARIABLE NAMES 

Variable names may be from 1 to 8 characters in length. The 
first character must be alphabetic. Names must be composed of 
letters and/or numbers only. They may contain no special char¬ 
acters or blanks. 

There can be a maximum of 4095 variables per record type. 


SUMMARY VARIABLE INITIALIZATION DURING RETRIEVAL 

The SIR system does not reinitialize the variable values of the 
summary record to missing values between each new summary record. 
Consequently, values may be carried forward from previous summary 
records if they are not overwritten purposely by the user. If 
variable initialization is required between each summary record 
the AUTOSET keyword must be used with the RETRIEVAL command ' 


FREE FORMAT COMMANDS 


Commands 

Command names must begin in column 1. 
Commands must be spelled correctly and 
completely. Command specifications may 
begin, following the command name.after 
one or more common delimiters(blanks or 
commas). 

Continuation lines 

Continuation lines can begin in any 
column after column 1 (i.e.,column 1 of 
continuation lines must be left blank). 

Sequencing 

If sequence numbers are present in col¬ 
umns 73-80, the NUMBERED YES command 
must be specified. 

Comments 

Comment lines require C or COMMENT in 
column 1. Text then follows one or 
more common delimiters. 


FIXED FORMAT COMMANDS 

SPSS Format As a subset of the above free format 

specifications, the SPSS control card 
format may be adhered to (i.e..commands 
in 1-15 and specifications in 16-72). 


DATA TYPES 


All non-string variables are handled as real valued in all cal¬ 
culations carried out by the SIR system. Data type options are 
important for flexibility in inputing data items, for their most 
efficient storage in the data file and for greater flexibility in 
data management. Sir supports the following data types: 

Type Examples Data file storage mode 


Integer 


Scaled 

integer 


-5, 375, 42 Packed into the smallest 

amount of storage space 
required for the largest 
integer that will fit in the 
number of columns specified 
for input of this variable. 

4.82, 3-97, 62.2 Converted to an integer for 
data file storage and recon¬ 
verted to decimal whenever 
retrieved from storage. 


Real 


4.75, 3.03 E-13 


One real computer word of 
storage for each data item. 









Date 


03/15/75,051273 


Stored as the integer numbe 
of days since Oct. 15, 1582. 


Time 0900, 09/00/32 

String "JOHN J. JONES" 

ARITHMETIC OPERATORS AND DELIMITERS 


Stored as the number of se 
conds since 12 midnight. 

Stored as a character strini 
variable with length 1 to 12! 
characters as specified by thi 
STRING LENGTH command. Defaul 
length is 20 characters. 


+ Addition 
- Subtraction 
* Multiplication 
/ Division 


l * Exponentiation 
( Open parenthesis 
) Closed parenthesis 
+ Concatenate two strings 


RELATIONAL OPERATORS AND LOGICAL CONNECTIVES 


EQ-Equals 

NE-Does not equal 

LT-Is less than 

LE-Is less than or equal to 

GT-Is greater than 


GE-Is greater than or equal tc 
AND-Both must be true 
OR-Either or both must be true 
NOT-Logical complement 


MATH FUNCTIONS 


ABS-Absolute value 
AC0S or ARCOS-Arc cosine 
AINT or TRUNC-Truncation 
ALOG or LN-Natural logarithm 
ALOG10 or LG 10-Logarithm to base 
ASIN or ARSIN-Arc sine 
ATAN-Arc tangent 


COS-Cosine 
EXP-Exponentiation 
RND-Rounded integer value 
SIN-Sine 

10 SQRT-Square root 
TAN-Tangent 

TANH-Hyperbolic tangent 


DATE AND TIME FUNCTIONS 
CDATE (x,f) 


CTIME (x,f) 

DATET (n1,n2) 

Columns 


Convert a date string to a Julian in¬ 
teger (the number of days since Oct. 15 
1582). 

COMPUTE DATE=CDATE 0 01/13/73' 

'MM/DD/YY') ’ 

Convert a time string x to the number 
of seconds since 12 midnight. 

COMPUTE SECONDS=CTIME 009.33 55' 

'HH/MM/SS') ’ 

Returns present date and/or time sub¬ 
string beginning at nl and ending at n2 
from the following 27 character current 
date/time string. 

Contents 


1-3 

9-5 

6-8 

9 

10-11 

12-13 

19- 17 
18-19 

20 - 21 
22 

23-29 

25 

26-27 


Day of the week (SUN, MON, TUE,...,SAT) 
Comma and blank 

Month of the year (JAN, FEB, ....DEC) 
Blank 

Day of the month (1 to 31) 

Comma and blank 
Year (e.g.,1977) 

Comma and blank 
Hour (1 to 12) 

Period 

Minutes (0 to 59) 

Blank 
AM or PM 
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sir - fu;jctio;is 


JULC 


(x) 


JULN (x,y,z) 


TIME (x) 


COMPUTE DATE=DATET (6,17) 
string 'DEC 25, 1978' left 
into the string variable DATE 
Christmas day 1978. 


places the 
justified 
if run on 


Converts the Julian integer x to 
14 1979 T datS 3trin S the form 


a 12 
'JAN 


WRITE 'THE 
(BIRTHDAT) 


PATIENT WAS BORN ON' JULC 


Returns an integer number of days from 
Oct. 15, 1582 to the date represented 
by the integer arguments x, y, and z 
which represent month, day and year 
respectively. 


Returns the integer number of seconds 
from 12 midnight to the time represent¬ 
ed by X. The argument x is an integer 
in the interval (0,240000) where the 
first two digits are hours, the second 
two minutes and the last two are se- 


STRING FUNCTIONS 
LEN (x) 


NUMBR (x) 

SBST (x,n1,n2) 

SRST (x1,x2) 


Returns the integer character length of 
the string variable x. 

IF (LEN(NAMEA) EQ LEN(NAMEB)) CNT=CNT + 


Converts the substring x to a binary 
number. 

Returns the substring of the string 
variable x beginning at character pos¬ 
ition n 1 and of length n 2 characters. 

Search string xl for the first occur¬ 
rence of string x2. Returns the pos¬ 
itive integer location where x 2 starts 
if it is found as a separate word. Re¬ 
turns negative integer starting 
location if x 2 is not found as a 
separate word. Returns zero if x2 is 
not found. 

IF (SRSTCNAME, 'JONES') GT 0) CNT=CNT+1 

This example increments the variable 
CNT by one if the last name 'JONES' is 
found as a separate word in the string 
variable NAME. 


WITHIN RECORD FUNCTIONS 
AM0D (x1,X2) 

CNT (xl,...,xn) 

FST (xl,...,xn) 

MAX (xl,...,xn) 


xl modulo x2. Returns the remainder of 
Xl/x 2 . 

Count the number of non-missing values 
among variables xl.xn. 

COMPUTE SUM=CNT(TEST1,TEST2,TEST3) 

Return the first non-missing value 
among xl,. . . ,xn.. 

Returns the maximum of the non-missing 
values among the variables x 1 ,...,xn. 
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MEAN (xl,,xn) 
MIN (xl,...,xn) 


Returns the mean of non-missing x’s. 


Return the minimum of the non-missii 
values among the variables xl,...xn. 

SRCH (v1,v2,v3) Search a table of increasing valu< 

variables beginning at variable vl ai 
ending at variable v2 for the locati< 
of the value of v3. 

If v3 matches an entry in the table,tl 
location of the matching entry is ri 
turned as a positive integer (countii 
from the beginning of the table). 

If no match is found, -n is return! 
where n is the location of the fir! 
value larger than v3. 

If v3 is greater than the largest valt 
in the table, -n is returned where n : 
one larger than the number of entri< 
in the table. 

Th.e indeterminate value is returned : 
the table vl to v2 is not in ascendii 
order. 

AGEGROUP = ABS(SRCH(GROUP 1,GROUP 10,AGE)! 

This example classifies AGE into one < 
ten age groups in the variable AGEGROl 

STDEV (x1,...,xn) Returns std. dev. of the non-missir 

x's. 

SUM (x1,...,xn) Returns the sum of the non-missing x'i 


ACROSS RECORD FUNCTIONS 

CNTR (x) Returns the number of non-missir 

values found for the variable x withj 
the record set. 

MAXR (x) Returns the maximum non-missing vali 

of the variable x within the recor 

set. 

MEANR (x) Returns the mean of the non-missir 

values of x within the record set. 

MINR (x) Returns the minimum non-missing valt 

of the variable x within the recor 

set. 

SUMR (x) Returns the sum of the non-missin 

values of the variable x within th 
record set. 

SUBTOTAL = SUMR (SYSBP) 

STDEVR Cx) Returns the std. dev. of th 

non-missing x’s within the record set. 

FSTR (x) Returns the first non-missing value fo 

the argument x from among values in th 
set of records processed. 


LSTR (x) 


Returns the last non-missing value fro 
among the set of records. 





SIR - SYSTEM WIDE COMMANDS 

SPECIAL FUNCTIONS 
COUNT (n) 


EXISTS (v) 


RAND (x) 


VARLAB (x) 


VALLAB (x) 


REQUIRED COMMANDS 


Returns 
records 
a case. 


the count of the number of 
presently in record type n for 


Determines whether 
for v has a missing 
if v has a missing 
value, 1 otherwise. 


the present value 
value. Returns 0 
or indeterminate 


^E XISTS(weig HT) EQ 1) SUMWT=SUMWT + 


This example checks 
variable WEIGHT has a no 
before attempting to add 


to assure the 
n-missing value 
it to SUMWT. 


Returns a 
interval 
porates a 
will not 
quence. 


P se Y do random number in the 
The function incor- 
random starting value and 
generate a repeatable se- 


Given the name of a 
function returns a 
containing the short 
has no label,a zero 
returned 


variable x, this 
character string 
label of x. If x 
length string is 


Returns a character string containing 
the value label of the current value of 
the variable x. A zero length string 
is returned if no value label exists 
for the current value. If the current 
value of x is missing, VALLAB returns a 
missing label. 


Most commands, in each SIR command set, are optional. The fol¬ 
lowing commands, however, must always appear in the specified 
command set. 

command set commands 


System Initialization FILE NAME, GET FILE, TAPE FILE or DISK 
RESTORE and PASSWORD 


Case Definition N of CASES, RECS PER CASE, RECTYPE COLS 

CASE ID and MAX REC TYPES 

Record Definition RECORD SCHEMA and VARIABLE LIST/INPUT 

FORMAT or DATA LIST 


SYSTEM WIDE COMMANDS================================== ======== - 

The commands in this section can appear throughout the SIR command set 
as indicated. 

COMMENT text 

C text 

This command allows placing comment material in the SIR command 
set for explanatory purposes only. No action results from this 
command and it may appear as often as necessary throughout the 
SIR command set. 

COMMENT THIS TEXT IS LISTED WITH THE COMMAND SET. 

DOCUMENT text 

Documentary text which is to be listed with LIST SCHEMA options. 
Allowed locations are prior to the first RECORD SCHEMA command in 
the schema case definition command set and between the RECORD 







SCHEMA command and the VARIABLE LIST command in each record type 
definition command set. 


DOCUMENT 


PRINT BACK 


THIS TEXT IS LISTED WITH LIST SCHEMA OPTIONS 
AND SERVES TO DOCUMENT THE STUDY. 

( v 

YES 

<NO 

CONTROL 

FORMAT 

REPEAT 

NOREPEAT 

Z 1 


This command controls the listing of the SIR command set. If no 
PRINT BACK command is used or if the YES option is used, SIR 
produces a listing of the user's input commands. In addition, 
SIR prints a variable-format correspondence table for each INPUT 
FORMAT or DATA LIST it encounters during a schema definition run. 
SIR also lists any lines of SIR code generated by DO REPEAT 
commands. 


The NO option surpresses the entire listing. The CONTROL option 
lists commands only. The FORMAT option lists the format corres- 
pondence table. The REPEAT and NOREPEAT control the listing of 
lines of code generated by DO REPEAT commands. 


nl[ ,n2 


3 


This command causes nl lines to be skipped in printing the com¬ 
mand set. If n2 is specified, n2 lines must remain on the page 
after nl lines are skipped. If there are insufficient lines 
remaining on the page, page eject occurs. 


TASK NAME text 


This command resets the subtitle line which is printed at the top 
of each page of output from a SIR run. Page eject occurs in the 
SIR output each time this command is encountered and the title, 
as specified by the RUN NAME command, and the new subtitle are 
printed before the next command is processed. 

TASK NAME THIS LABEL CAN BE UP TO 59 CHARACTERS LONG. 
STRING LENGTH n 

This command declares the length in characters of every string 
variable whose initial declaration follows this command until 
another STRING LENGTH command is encountered. This command can 
appear as often as necessary to specify different string variable 
lengths. 20 characters are assigned to each string variable by 
default and a maximum length of 125 characters is allowed. 


STRING LENGTH 36 
COMPUTE 


INPUT MEDIUM 


CHARVAR = "THIS VARIABLE CAN HOLD 36 CHARACTERS". 

[cards "1 

< filename/ 

Specifies the source of any data to be read into SIR during a 
run. If CARDS is specified, the data must immediately follow the 
data read command in the command set itself and be terminated by 
a ?Tj eri u of_recorcJ marker. If a file name is specified, the data 
Wl11 from a l° ca l file found with that name. More than 

one INPUT MEDIUM command can occur in a SIR command set to allow 
input from multiple data files as necessary. This command must 
appear in the command set before the command which causes the 
data to be read. 


GET FILE 
PASSWORD 
INPUT MEDIUM 
READ INPUT DATA 
INPUT MEDIUM 
UPDATE REC 
FINISH 


JSTUDY 

«*»»»* 

DATA 1 

ACCEPT, ERRFILE = ERR0RS1, ALL 
DATA2 

ACCEPT, ERRFILE = ERR0RS2, ALL 
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SIR - SYSTEM INITIALIZATION COflMANDS 


DO REPEAT stand-in variable name = variable list/ 

[stand-in variable name = variable list/...] 
one or more replica commands 

END REPEAT 

This command provides a means of specifying sets of repetitious 
commands in a shorthand fashion. 

DO REPEAT K=TEST1 TO TEST10/ 

IF (EXISTS(K) EQ 1) COUNT = COUNT + 1- 

SUM = SUM + K 

END REPEAT 

This example accumulates a count and a sum of the non-missing 
test scores among ten test score variables. 

FINISH 

The last command in the users command set must be FINISH. 

ERROR LIMIT /NONE 

WARNING LIMIT (n 

A limit can be specified to the number of error or warning mes¬ 
sages respectively generated by SIR. NONE imples no messages. Up 
to 50 messages are printed by default. 

ERROR LIMIT 10 


SYSTEM INITIALIZATION C0MMANDS = = = = = = = = = = = = = = = = = = = = = = = = :: = = = = - =; = = = - = = --. 

The commands described in this section may appear only once at the 
beginning of the user's SIR command set. 

EDIT 

This command informs SIR that the purpose of this run is to check 
the syntax of the command set only. No commands are actually 
executed by an edit run. This command must appear at the very 
beginning of the command set when used. 

FILE NAME name »» 

To be used with initial SIR creation run only, this command 
specifies the name to be assigned to the SIR files for this data 
base. The name can be from 1 to 6 characters long with the first 
character being alphabetic. 

GET FILE name 

This command replaces the FILE NAME command for SIR runs on a 
previously defined data base. The name must correspond to the 
name specified in the FILE NAME command used when the data base 
was originally created. 

TAPE FILE name *» 

When SIR processing is to be done directly from tape, this com¬ 
mand is used in place of the GET FILE command. In this mode, 
cases are processed sequentially directly from tape and not 
stored to disk. This command assumes the data base has been 
placed on tape by the TAPE STORE utility. The TAPE FILE command 
is not implemented for SIR Version 1.1. Implementation of this 
command is scheduled for SIR version 1.2. 




DISK RESTORE 


name, FILENAME = local file name 
,LOADING=n [/NOREWIND] 


If the SIR files reside on tape, this command can be used ir 
place of the GET FILE command to restore the SIR files to disk 
and assign them to this run. Name is the SIR file name, FILENAME 
is the tape file name. The LOADING factor, n, is a number 
between 0 and 1 that indicates roughly how much of the amount of 
data expected has been input to the file thus far. If not 
provided, n is assumed to be .5. 


**The appropriate one of these four commands is required for every 
batch SIR run. 


NUMBERED j NO j> 

If this command is present with the YES option, columns 73-80 of 
the command lines are reserved for line numbering. If not 
present the NO option is assumed by default and columns 73-80 are 
available for command specifications. 


JnoejectI 

PAGESIZE j^n J ,Columns per line 


Page control for the listing of the SIR command set is controled 
with this command, n specifies the number of lines to occur on 
each page of output. If NOEJECT is specified, no page ejection 
is done. Default is 60 lines/page. Cols/line may be within the 
interval 72-136. 


PASSWORD password 


(Required) 


The password is a 1-7 character string which must be provided for 
every access to the SIR system. The password is set at initial 
data base creation by the password given with this command on the 
creation run. 


RUN NAME text 


This command sets a main title line of up to 59 characters which 
will be printed at the top of each page of output from SIR run. 

SECURITY read password, write password 


This command is required using the read/write passwords assigned 
to the user by the study manager if read and write security has 
previously been imposed on the data base to be accessed. If no 
security has been previously imposed, with READ SECURITY and 
WRITE SECURITY commands, the SECURITY command is unnecessary. 

PASSWORD BSTUDY 
SECURITY JOHNR,JOHNW 


READ SECURITY (security level) password [(security) password] 
WRITE SECURITY (security level) password [(security) password] 

Assigns security levels (0 lowest, 30 highest) to specified pass¬ 
words. These passwords may be assigned to users for use in 
SECURITY commands when accessing SIR. If these commands are not 
used to establish security on a data base, all users obtain high¬ 
est security level access to that data base by simply providing 
the correct password on the PASSWORD command. 


SCHEMA CASE DEFINITION C0MMANDS=================== =================== - 

The case definition commands are used when defining a new SIR data 
base. These commands supply information about the nature of the data 
base as well as general information that applies to each case in the 
data base. 

CASE ID name (Required) 
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SIR - SCHEMA CASE DEFINITION 


Specifies the name of the variable which will identify the case 
It must be present with this same name in the variable 1 = 1 - 
every record type defined later and be consistent in type 


cord type defined later 
throughout the data base. It may be of type integer reai^or 
character string. This variable will, by default be included in 
the CIR to avoid its storage with each record in the data file 
The values of this variable must be unique for each case in the 
SIR data base. If the CASE ID is of type integer or real it mav 
take on negative or positive values. ’ Y 

CASE ID IDNUMBR 


COMMON LIST 


variable list 


Variables named in this command will be included in the common 
information record (CIR). The CIR is a single record created for 
each case. During SIR data retrievals, the CIR contents are 
always available regardless of which record type is being 
processed within a case. B 


COMMON LIST AGE,SEX,ETHNIC 

MAX KEY SIZE n 

Specifies the maximum length in bits that will be required for 
each SIR index key. This command is necessary only when a record 
type will be entered after the initial data base creation which 
will require a larger index key size than any record type enter¬ 
ed initially. The size n may be determined by the formula: 

SIZE = A + keysizes 

where A = log2 (// of different record types) + 1 
and 

Keysize = 2+log2 (largest possible value) for 
integer keys 

Keysize = 60 for floating point keys 

Keysize = 6*(max. II of characters) for string keys 

The keys for each record type are the CASE ID, record type number 
and any sort variables. The largest SIZE calculated among all 
record types is the value to be used with the MAX KEY SIZE 
command. 


MAX INPUT COLS number 


This command is required when one or more of the record types in 
the data base has more than 80 columns per input line. If it is 
omitted, the assumed record length is 80 columns. 

This command specifies the length of the largest input line of 
any record type in the data base. Once specified, each input 
line of each record in the data base is viewed as being this 
length. 

The number specified is always rounded up to the nearest 10 units 
The maximun allowable number of input lines for a record type may 
be calculated by dividing the number specified by this command 
into M095 and truncating the result to an integer. 

MAX REC COUNT n 

Specifies the maximum number of records that will occur in any 
record type within any case in the SIR data file. A maximum of 
1023 is assumed by default. 

MAX REC COUNT 300 

MAX REC TYPES n (Required) 

This command specifies the maximum number of record types that 
will eventually be specified for this data base. There is es¬ 
sentially no limit to the number of record types that can be 
defined. This number is normally set somewhat higher than the 
number of record types assumed for the study initially. 
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MAX REC TYPES 


30 


N OF CASES n (Required) 

Specifies the maximum number of cases that eventually might be 
contained in this data base. If the number is not known, a high 
value should be specified. 

N OF CASES 5000 

SEQUENCE COLS beginning columnt , ending column] 

For those record types which will require more than one data card 
for input, this command specifies the column(s) of the data cards 
which will contain the card sequencing number. Only record types 
requiring multiple cards per record are checked in these columns. 


SEQUENCE COLS 79,80 

SEQUENCE COL 80 

RECS PER CASE n 

Specifies the average number 
each case in the completed data 
and the number specified in the 
up tables for index keys. 


(Required) 

of records that are expected for 
base. The product of this number 
N OF CASES command allows setting 


RECS PER CASE 150 


RECTYPE COLS 


beginning numbert.ending number] (Required) 


This command specifies the column(s) in which the record type 
identification number will be found on each incoming data record. 


RECTYPE COLS 6,7 


TEMP VAR LIST 


variable list 


This command identifies a set of variables which will be used as 
scratch or temporary variables in the process of creating other 
variables later using IF or COMPUTE commands. Variables named in 
this list will not be stored in the data file. CIR variables, 
including the CASE ID, cannot appear in this list. 

TEMP VAR LIST TEMPI TO TEMP25 


SCHEMA RECORD DEFINITION COMMANDS=========r=========================== 

The commands in this section are used to describe the contents and 
structure of individual record types. There will be one set of record 
definition commands for each record type in the study. 

During a Schema Definition run, the Record Definition Commands that 
are used must be placed after the Case Definition Commands. Also the 
Record Definition Commands of each record type must be grouped to¬ 
gether. Different record types, however, do not need to be arranged 
in sequential order nor do they all need to be defined in a single 
run. 

Within a set of Record Definition Commands, there is a specified or¬ 
dering of commands that must be followed. The precedence groups among 
these commands are listed below. Within a group, the commands may 
appear in any order. 


group 

command 

group 

command 

1 

2 

RECORD SCHEMA 

SORT RECORDS 

5 

COMPUTE 

IF 


SEQUENCE CHECK 

REC SECURITY 

MAX REC COUNT 

6 

VALID VALUES 
VAR RANGES 


DOCUMENT 

7 

MISSING VALUES 
VAR SECURITY 




SIR - SCHEMA RECORD DEFINITION 


VARIABLE LIST-INPUT FORMAT 
or DATA LIST 

DATE VAR LIST 
TIME VAR LIST 
SCALED VAR LIST 
INTEGER LIST 


VAR LABELS 
VALUE LABELS 


ACCEPT REC IF 
REJECT REC IF 


ACCEPT REC IF 
REJECT REC IF 


(logical expression) 
(logical expression) 


For data editing at input, logical consistencies which must be 
true between variable values on incoming records (or values on 
incoming records with data already in the CIR record) can be 
specified with these commands. The data record is rejected if 
the logical expression specified is false for ACCEPT REC IF and 
rejected if true for REJECT REC IF. 

MPUTE variable=arithmetic expressiontjvariables 

arithmetic expressiont;...]] 

OR 

variable list(index)=arithmetic expression!! • ] 

OR 

variablesvariable list(index) [;...] 

The COMPUTE command has all of the capabilities of the SPSS 
COMPUTE and also provides both the ability to include multiple 
expressions in a single COMPUTE separated by semi-colons and an 
indexing capability to select a specified member for computation 
from a list of variables. The index capability is available onlv 
for COMPUTE commands used in SIR retrievals. J 


COMPUTE 

COMPUTE 

COMPUTE 


COUNT = COUNT + 1 

AGEGRP=ABS(SRCH(AGE 1,AGE5,AGE)) 

SUMS01 TO SUMSQ6(AGEGRP) = SUMSQ TO SUMSQ6(AGEGRP' 
+ SYSBP*#2 

COMPUTE AGEGHP=ABS(SRCH(AGE1,AGE5,AGE)):SUMSQ1 TO SUMSO( 

(AGEGRP)=SUMSQ1 TO SUMSQ6(AGEGRP) + SYSBP**2 

TA LIST (II of lines in the record type)/ line number 

variable list-columns [(variable type)] 

[variable list-columns [(variable type)]] 

[/line number...] 

This command provides an alternative to the VARIABLE LIST and 
INPUT FORMAT commands. Variable type may be A or I for character 
or integer. Floating point variables with fixed decimals are 
indicated by an integer in place of the variable type indicating 
the number of places to shift the decimal from its assumed right 
justified location. If variable type is omitted, SIR assumes the 
variable is floating point with decimal included or rieht 
justified if not included. B 

DATA LIST 


DATE VAR LIST 


(3)/1 

CASE ID 1-5 (A), AGE 7-9 (1), BIRTHDATE 20-25(A) 

variable list(date format)[/var list(format)] 

Variables named in this command will be converted to Julian in¬ 
tegers and stored in the data file as the number of days the date 
represents from Oct 15, 1582. 3 

DATA VAR LIST BIRTDATE (’MMIDDIYY’)/MARRIED(•DDIMMIYYYY 1 ) 

’s included in the date format indicates that those columns 
, ignored at input (e.g.,-MMIDDIYY' reads 02/14/76 as 
021N76. A variable included in the DATE VAR LIST must be read as 
a character string with an A format. 

(logical expression) computation 
[;computation..] 


The IF 
commands 


command provides the ability to do conditional COMPUTE 
m that the computation is done only if the logical ex- 
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pression specified is true, 
tical to the COMPUTE command. 


Computation capabilities are iden- 


IF (EXISTS(SYSBP) EQ 1) COUNT = COUNT + 1 

INPUT FORMAT (format specification) 

This command is required immediately following the VARIABLE LIST 
command and specifies the format for reading the data from the 
satisf V the variable list. Valid format items 
X,I,F,A,T,(,) and ,. Repetition factors may also be used. Free 
format input is not supported. The SPSS FIXED keyword is 


(I't,1X,A25,5I3,/,T6,2(F5.2 1 IX)) 


INPUT FORMAT 

This example tells SIR to read the first four columns as an 
integer, skip one column, read 25 columns as an alphebetic string 
and five variables as 3 column integers, skip to the next input 
line, tab to column 6 and read two fields five columns wide as 
decimal numbers with two decimal places to the right. These last 
two fields are seperated by a single blank. 


INTEGER LIST 


variable list (min, max)[ /var list (min, 
max)/.. . ] 


Variables included in this command will have an integer identity 
for data file storage and output. For calculation purposes they 
will be considered real valued and will be truncated for storage 
aad °, u ;?i! t listin 8- Variables specified in the VARIABLE LIST or 
DATA LIST commands cannot be included in this command since they 
are already specified as to type. Consepuently, this command is 
for use with variables created by COMPUTE and IF commands. 

INTEGER LIST SCORE(1,500)/INDEX(0,20) 

MAX REC COUNT n 


This command allows the user to specify the maximum number of 
records that will be found in a particular record type for any 
oase ,, the data base. The number is set at the MAX REC COUNT 
specified for the case if this command is not included in the 
record definition set. 


MISSING VALUES 


variable list (missing value list)[ /var list 
(missing value list)/...] 


Up to a maximum of 3 separate values may be specified for each 
variable to signify missing data. The keyword BLANK may be 
included to indicate that blanks imply missing. If a missing 
valae • not 3 P eciflad for a given variable,BLANK will be the 
default missing value. If other values are specified as missing 
values for a variable, then a blank field on the input data 
record will be read as zero. 


MISSING VALUES SEX TO INCOME(BLANK)/AGE(98,99) 

Missing values may also be specified for a string variable. For 
example, 


MISSING VALUES BIRTHDAY( "99/99/99 ") , STATECZZ") 


If the 
value, 
missing 


string variable contains more characters than the missing 
then a match between the left most characters and the 
value will be read as missing. 


RECORD SCHEMA number[,name] 


This command must appear . at the beginning of the record 
definition command set for each record type. The record type 
number specified for a record type will normally be expeeted^o 

?hp ea RPT??n? a rm= n00ming J data line ln the oolomn(s) specified in 
the RECTYPE COLS command. To enter data without the record type 
number on the input lines, see the data input commands. The name 
is an optional 1-8 character label that can be given to the 
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excessive amounts of computer time then required for data entry. 


VALID VALUES SEX(1,2)/SCORE(1,5,7,9) 

Categorical character variables (variables read as characte 
strings but which are to be treated as discrete integer variable 
within SIR), will be stored in SIR as sequential integers i 
included in a VALID VALUES command. 

VALID VALUES SEXCM’.'F') 

VALUE LABELS SEX ("M") MALE ("F") FEMALE 

The value 1 will be stored in the data file for male and 
the value 2 for female. Notice, however, that when 
specifying MISSING VALUES or VALUE LABELS for the variable 
SEX, the character values are refered to. 

VALUE LABELS variable list (valuel) labell (value2)label2,. 

[/var list (valuel) label!,...] 

Up to a 20 character label may be specified for some or all o 
the categories of each variable. 

VALUE LABELS SEX{"M")MALE("F")FEMALE/OCCUPTN(1)BLUE COLLAR 

(2)HOUSEWIFE(3)WHITE COLLAR(4)0THER 

During retrievals the value label may be obtained for printing i 
reports by use of the VALLAB function. 

VAR LABELS variable name, label[/variable name, label/... 

Variable labels may be as long as desired by the user. Only th 
first 40 characters or up to the first •;’ whichever occurs firs 
will be transmitted to an SPSS save file or the code book listin 
when the REGULAR listing option is requested. Semicolons in th 
label also cause skipping to a new line in the DETAILED cod 
book listing and can, therefore, be used for format control i 
the label. 

VAR LABELS VITALCPY, PATIENTS VITAL CAPACITY IN CC; 

THE VITAL CAPACITY IS MEASURED; 

AT THE BEGINNING AND END OF THE EXAM/ 
HEARTRT, PATIENTS HEART RATE IN BEATS PER MIN. 

During retrievals, a variables label may be obtained by use o 
the VARLAB function. 

VAR RANGES variable listtmin, max)[/var list...] 

Values falling outside the ranges specified by this command ari 
rejected at data input. 

VAR RANGES WEIGHT (30.00, 160.00)/HEIGHT (20,80) 

VAR SECURITY variables listtread level, write level) 

[/var list_] 

Read and write access levels are assigned as specified from ( 
(lowest) to 30 (highest) to those variables included in the lisl 
In order to gain access to a variable, the user's security level- 
word specified in the SECURITY command must be at least as higl 
as the variable's security level. A level of zero is assigned tc 
each variable by default. 

VAR SECURITY NAME, ADDRESS, PHONE, S0CSEC (30,30) 

VARIABLE LIST name [name ...] (Required) * 

All variables to be defined within this record type must be namec 
by use of either this command or the DATA LIST command. Name: 
consist of 1 to 8 alphanumeric characters and must begin with : 
letter. The shorthand notation VAR 1 to VAR 10 will create ter 
variables VAR!, VAR2,...,VAR 10. 
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IDNUMBR, NAME, ADDRESS, VR TO V30 


SIR - SCHEMA LISTING 

VARIABLE LIST 

•either a VARIABLE LIST or a DATA LIST command must be included with 
each schema definition command set. 


SCHEMA LISTING COMMANDS=============================================== 

The schema listing commands provide the user with a flexible set of 
options for producing reference documentation for the study. The 
listings made by the LIST SCHEMA command will fit within 8 1/2 x 11 
inch page format and are suitable for binding. 

LIST SCHEMA keyword listC r til (var iabl e list) ,rt#...] 

The keyword list consists of one or more of the options 

STRUCTURE, LABELS, REGULAR or DETAILED. 

The STRUCTURE option produces a report describing the case and 
record structure of the SIR data base. Any DOCUMENT commands 
defined during data schema definition are also listed in this 
report. 

The LABELS option generates a report containing the variable 
names and short variable labels for the record types and vari¬ 
ables specified. 

The REGULAR option results in a listing containing variable 

names, variable type, position in the input record, short vari¬ 
able labels, range or valid values, read and write security 

levels and missing value designators. 

The DETAILED codebook listing contains all of the information in 
the REGULAR listing with the addition of the full variable label 
and all value labels. 

When variable information is printed, each variable is tagged 
with a label indicating its origin and location within the 
record. Labels are of type Vnnn, Cnnn or Tnnn. The nnn is a 

number indicating the location of the variable within its record. 
V, C and T represent record variables, CIR variables, and 
temporary variables, respectively. 

A SHORT or LONG keyword may be used to begin the keyword list 
(except SHORT is not allowed with STRUCTURE). The SHORT key¬ 
word suppresses the listing of headings and DOCUMENT output. LONG 
is default. The SHORT keyword is useful for listings at 
terminals. 

LIST SCHEMA SHORT, LABELS, REGULAR 1,3,7,99 

LIST SCHEMA LABELS 1 (SEX TO EDUCATN) 

LIST SCHEMA STRUCTURE, LABELS, DETAILED 0,1,2,3,9 

A zero included in the record type number argument list produces 
a separate listing of the CIR record schema. 

WRITE SCHEMA FILENAME = filename [/NOREWIND] [RECTYPES = 

record type list] 

Writes a card image code book schema definition command set to 
the specified file. 

LIST STATS 

Lists current general information on the data base including: 
last update, last retrieval, number of accesses, records in each 
record type, etc. 


SCHEMA MODIFICATION = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = 

SIR provides the means of modifying various case and Record definition 
specifications as required. If data has not yet been entered into the 


16 







data base, code book definition sets may be re-entered as many times 
as necessary with corrections. Once data is in the data base, any 
change in the schema definition which might affect the data base will 
cause the record type involved or possibly the entire data base to b 
"locked" so that it cannot be accessed. 


unlock a locked record type, 
uted in order: 


the following three commands must be 


1. TAPE STORE 

2. PURGE SIR FILE 

3. DISK RESTORE 


Note: Data Base Administrator security is required to lock or unlock 

a locked record type. 


Often changes in the schema definition are most easily accomplished by 
re-entry of the complete schema definition set for the entire code 
book or at least for an entire record type. The SIR Interactive 
Editor provides several facilities that make such record schema 
modification rather simple. For example to add two variables RDATE 
and SDATE to record type 3»the following procedure could be followed: 

X,WRITE SCHEMA FILENAME = ABC/RECTYPES = 3 

This writes the entire set of record type 3 schema commands to 
ABC. 


READ,ABC 

This reads the contents of ABC into the SIR Editor. Using the 
Editor commands, The text is now modified to include RDATE and 
SDATE. 


This executes the current text which re-defines record type 3 as 
desired to include the new variables. If data is in the data 
base, record type 3 will now be locked and the three step 
procedure described above will be required to unlock it. 

The only commands which may not be changed in an existing SIR data 
base are CASE ID and MAX INPUT COLS. To change either of these, the 
user must completely re-enter the data base from raw data. The WRITE 
SCHEMA and SIR FILE DUMP commands, however, are helpful should this 
ever become necessary. 


Changes in the variable definitions in data bases which already have 
data in them can have major implications on the data already in the 
data file. Consequently, the user is advised to read carefully the 
section on schema modification in the SIR Users Manual and proceed in 
this area very carefully. 





MODIFY SCHEMA rectype number or name[,new name] 

This command allows the user to modify those portions of an 
existing record schema that have no effect on the structure of 
the data file without re-entering the entire schema. The new 
name field is optional. If a new name is entered, it replaces 
the old record type name. 

The following set of commands are the only ones that are allowed 
to follow the MODIFY SCHEMA command. They must be entered in the 
order listed here: 


SEQUENCE CHECK 
REC SECURITY 
DOCUMENT 
VAR SECURITY 
VAR LABELS 
VALUE LABELS 
ACCEPT REC IF 
REJECT REC IF 
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SIR - DATA ENTRY/UPDATE 

Changes entered for VALUE LABELS, ACCEPT REC IF and REJECT REC IF 
in the MODIFY SCHEMA completely replace those presently entered. 
This means that all value labels must be entered for the 
variables to be changed and that if the ACCEPT REC IF or REJECT 
REC IF command is included, the entire set of correct consistency 
checks must be re-entered. 

{ compute') 

VALID 
RANGES 
MISSING 
BOOLEAN J 

The purpose of the CHECK command is to provide the user with the 
means of reverifying existing data against the new specifications 
entered during schema modification. The various keywords will 
lead to reverification of data in the data base on all COMPUTE, 
IF, VALID VALUES, VAR RANGES, MISSING VALUES, ACCEPT REC IF and 
REJECT REC IF commands in the modified schema respectively. Data 
is not reverified unless this command is present. 

If used, the CHECK command must be placed immediately before the 
VARIABLE LIST or DATA LIST command in the revised schema set. 


DATA ENTRY/UPDATE COMMANDS = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = :: = = = = = = = = = = = 

Once the SIR data base has been defined in a Schema Definition run, 
the user can input, change or delete data in the data base using the 
commands described in this section. 

READ INPUT DATA RECTYPE = number or name 

[/NOSEQ][/NONEW][/ACCEPT] 

[/ERRFILE = filename [,NOREWIND][,ALL]] 

Data read from the location specified in the latest INPUT MEDIUM 
command will be added to the data file. Both new and replacement 
data are allowed. ACCEPT causes records with invalid data items 
to be stored to the data file with the invalid data items re¬ 
placed by missing values. The following exceptions will always 
cause complete record rejection: 

1. If the error is in a consistency check 

2. If the error is in a record identifier 

3. If a line is missing or out of order in a multiple line record 

By default all records with error values are rejected before data 
file entry. The ERRFILE keyword provides a file where copies of 
rejected records are written. The NOREWIND option prevents ERR¬ 
FILE from being rewound. ALL causes every record with errors to 
be written to the error file, even those written to the data file 
due to the presence of ACCEPT. 

The RECTYPE keyword implies that all of the data lines in the 
data input file are of the same type and do not, therefore, need 
the record type number on the data lines. (Note - the user must 
still specify a RECTYPE COLS for the record at schema definition 
time) . 

The NOSEQ option instructs SIR to bypass sequence checks for 
multiple line records. The NONEW keyword prevents the creation 
of new cases during a data entry run. 

READ INPUT DATA ACCEPT ERRFILE = ERRORS, ALL 

ADD REC RECTYPE = number or name 

[/NOSEQ][/NONEW][/ACCEPT] 

[/ERRFILE = filename [,NOREWIND][,ALL]] 

Add new data only, no overwriting of old data allowed. 

REPLACE REC RECTYPE = number or name 

[/NOSEQ][/ACCEPT] 

[/ERRFILE = filename[,NOREWIND][,ALL]] 






Replace old data only, no new data allowed. 


UPDATE REC [ERRF1LE = filename[,NOREWIND][.ALL]] 

[/NOCOMP][/NOBOOL] 

[/MISSCHAR = character] 

[/RECTYPE = number or name] 

C/NOSEQ][/ACCEPT] 

Insert updating data items into existing records. Updating dati 
lines need only contain identfiiers and the updated data items ii 
the proper columns. The new value(s) will replace the ol< 
value(s) . 

By default, SIR will perform all COMPUTE commands defined in thi 
schema definition for this record type that are appropriate foi 
the data lines coming in from the input data file. If the NOCOMI 
keyword is specified, none of the input record COMPUTE command: 
are executed. 


By default, SIR also performs all of the consistency check: 
defined for this record type on the incomming data. If th< 
NOBOOL keyword is included, ACCEPT REC IF and REJECT REC II 
commands are not executed. 


The MISSCHAR option allows a character to be specified which will 
be recognized as a missing value indicator for a data item. It 
the character specified (blank not allowed) is found in the left 
must possition of a data item field, that data item is set tc 
missing in the data base. 


MODIFY REC rt/l, case id value, sort id], sort id2 ,... 

The specified record is made available for modification. COMPUTE 
and IF statements following the command act on this record. 

MODIFY REC 3 , 1324 , 4 , 1 

COMPUTE SCORE = 5 

IF (SEX EQ 1) MARSTAT = 3 


MODIFY REC IF rt II (logical expression) 


Makes all records in the specified record type, for which the 
logical expression is true, available for modification as il¬ 
lustrated in the MODIFY REC command. The logical expression may 
reference variables in the named record type and the CIR. 


MODIFY CASE case id value 


Allows modification of the CIR record of the specified case. 

MODIFY CASE IF (logical expression) 

Allows modification of all case CIR records for which the lo¬ 
gical expression is true. Only the CIR variables may be refer¬ 
enced in the logical expression. 

DELETE REC rt/), case id value, sort idl, sort id2,... 

Purge the identified record from the data file. 

DELETE REC IF rt/) (logical expression) 

Purge all records in the identified record type which satisfy the 
logical expression. The record type and CIR variables may be 
referenced in the logical expression. 

DELETE CASES case idl, [case id2,...] 

Purge the specified cases from the data file. 

DELETE CASE IF (logical expression) 

Purge all cases satisfying the stated logical expression. Only 
CIR variables may be referrenced in the logical expression. 





SIR - BATCH RETRIEVAL 

DELETE CASE IF (CASEID GT 1242 AND CASEID LT 1401) 

EVICT REC [RECTYPE = name or number] 

[/NOSEQH/EVICTCIR] 

[/ERRFILE = filenamet.NOREWIND]] 

This command provides a means of removing the data records found 
on the input data file from the data base. If a set of data 
records were added to the data base which the user later wishes 
to remove, this command provides the necessary tool for that 
purpose. 

The data records on the INPUT MEDIUM file must contain (at a 
minimum) the case id value, the record type number (unless the 
RECTYPE option is used) and values for all the sort ids, if any 
are called for. All lines of multiple-line records must be input 
with complete sequencing information. Any other data on the 
records is not looked at by SIR. The EVICTCIR tells SIR to 
eliminate a case completely from the data base (ie..remove the 
CIR record) when all of the other records have been evicted. 


BATCH RETRIEVAL C0MMANDS================ 


A SIR retrieval program consists of one or more "blocks" of retrieval 
commands. Within each block, retrieval operations are performed on 
records of a single type in the data file or on summary records. 

Each block in a retrieval program ends with the beginning of a new 
block or with an END statement. When END statements are used as block 
terminators, it is permissible to insert certain statements between 
blocks. 

The general structure of a SIR retrieval program is as follows: 


0 or 1 
0 or 1 
0,1 or more 
0 or 1 
0,1 or more 
0 or 1 
0,1 or more 
0 or 1 


RETRIEVAL block 
FOR EACH CASE block 
PROCESS REC blocks 
FOR EACH REC blocks 
PROCESS REC blocks 
SUMMARY blocks 

procedures (statistics, reports, etc.) 
AFTER RETRIEVAL block 


The RETRIEVAL block is used for defining data tables, initiallizing 
summary variables, printing report preamble information and defining 
report heading and footing information. 


The FOR EACH CASE block is used to limit the number of oases accepted 
in a retrieval. Without this block, SIR assumes that all the oases in 
the data base are to be processed. 

The PROCESS REC and FOR EACH REC blocks are used to select records 
from the data base and create summary record variables from these 
records. The SUMMARY block is used to perform final operations on the 
variables in the summary record and to make final determinations about 
whether or not specific summary records will be used by the procedures 
which follow. 


The type of retrieval performed by SIR depends on the presence or ab¬ 
sence of the FOR EACH REC block in the retrieval program. 

When a FOR EACH REC block is present, a summary record is created for 
each record of the type specified on the FOR EACH REC command subject 
to the selection criterion (SELECT REC IF or REJECT REC IF) command 
included in the block. 


If there is no FOR EACH REC block in the retrieval program, a single 
summary record is created for each case selected from the data base. 
In this type of retrieval, the summary record for a specific case may 
be suppressed by using the SELECT CASE IF or REJECT CASE IF command. 


The SUMMARY block signals the end of summary variable creation for the 
current retrieval. All the variables that appear anywhere in the re- 
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trieval program become summary variables and appear in the summary re¬ 
cord in the order that they are created in the retrieval program. 

T h ® AFTER RETRIEVAL block is used for writing report endini 
information K 1 


SIR provides many different 
appear in the program blocks, 
mands grouped by type: 


kinds of retrieval commands which aa< 
The following is a list of these com- 


COMMAND TYPE NAME 


Case Control 
Record Selection 

Summary Record 
Selection 

Summary Variable 
Initialization 


SELECT CASE IF 
REJECT CASE IF 

SELECT REC IF 
REJECT REC IF 
SELECT FIRST 
EXIT BLOCK IF 

SELECT REC IF 
REJECT REC IF 
PERFORM PROCS 

PRESET 


Summary Variable 
Definition 


Report Generation 


AFTER RETRIEVAL 


COMPUTE 

IF 

MOVE VAR LIST 
RECODE 

ASSIGN MISSING 
VAR LABELS 
VALUE LABELS 
VAR RANGES 

WRITE 
WRITE IF 
HEADING 
FOOTING 
PAGE EJECT 
BLANK LINES 


This statement initiates a retrieval 
located at the end of the retrieval 
cedure options. COMPUTE, IF, HEADING. 
BLANK LINE(S) AND PAGE EJECT commands 
will be executed once at the end of the 


program block which may be 
command set after all pro- 
FOOTING, WRITE, WRITE IF, 
which follow this command 
retrieval. 


AFTER RETRIEVAL 

COMPUTE 

WRITE 


ASSIGN MISSING 


AVERAGE = GRTOTAL/GRCOUNT 
*,'THE AVERAGE PERFORMANCE=' 

AVERAGE (F8.3) 

variable list (missing value)[/var list(value) 


bv S1 COMPUTF L al ?r S whioh to be assigned to variables created 
0M ^ UT ^ ° r commands when indeterminate results occur. The 
assigned value will be used whenever the variable is written out 

for a a variabL t0 the°d e f Ur ?; If "° ASSIGN MISSING b°™>and is used 

lor a variable, the default missing value is BLANK. 


ASSIGN MISSING SCORE, AVESCORE(99)/TOTAL(999) 


BLANK LINE 
BLANK LINES 


Inserts one or more blank lines in a 
occurrence of BLANK LINE skips one 1 
LINES gives the number lines to be 
eject, the remainder of n is ignored. 


retrieval output 
ine. n specified 
skipped. If n c 


list. Each 
with BLANK 
auses page 
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SIR - BATCH RETRIEVAL 


COMPUTE see SCHEMA RECORD DEFINITION COMMANDS 


PERFORM PROCS [(logical expression)[computation!] 


May appear within a SUMMARY block. Normal output of the summary 
record is suspended and a single copy of the summary record is 
written out for each PERFORM PROCS command encountered (if a log¬ 
ical expression is present, it must first evaluate to true). 
Makes possible output of multiple copies of a summary record with 
computations controlling the content of each one. 


SUMMARY 

DO REPEAT I = T1 TO T5/ 

PERFORM PROCS (EXIST(I) EQ 1) DRUG = I 

END REPEAT 

END 


END 


Terminates a retrieval block. May be omitted in cases where the 
end of the retrieval block is obvious. Required if the user de¬ 
sires to place commands between two blocks. 


o 


EXIT BLOCK IF (logical expression) 


This command may only appear in a PROCESS REC or FOR EACH REC 
block. It causes termination of record processing within a 
retrieval block when the logical expression is true. Execution 
then jumps out of the block if the command is in a PROCESS REC 
block and jumps to the. next case if the command is in a FOR EACH 
REC block. 


PROCESS REC 
SELECT REC IF 
EXIT BLOCK IF 
MOVE VAR LIST 
END 


A 

(SEX EQ 2) 

(VISTDATE GT JULN(01,01,77) 
BLOODPRS TO VITALCPY 


FOOTING output specification 

The output specification will be written at the bottom of each 
page of output from a retrieval prior to page eject. The output 
specification is as described in the WRITE command. 


FOOTING 
FOR EACH CASE 


100X, 'PAGE.'.PAGE 


{ COUNT = totalt ,increment[,start]] 

LIST = case idl[,case id2case idn 
THRU case idm]...] 

SAMPLE = fraction [ .seed] 

Located after the RETRIEVAL command and prior to any other re¬ 
trieval program blocks, this command provides information about 
which (or how many) cases are to be processed by this command 
set. 

FOR EACH CASE COUNT = 25 

FOR EACH CASE COUNT = 300, 1, 3050 

FOR EACH CASE LIST = 931, 1122, 2A00 THRU 2851, 5001 

FOR EACH CASE SAMPLE = .01, 1785 

f WITH ] r -s 

UNTIL |UNTIL| 

FOR EACH REC rt f [ [ < THRU Msort id list) ] Msort id list)] 

name FROM THRU 

(jlFTERJ ^ J 

Signifies the beginning of a master processing block and 
initiates reading of each record within the specified record 
type. A summary record is produced for each record in the record 
type which is selected by the SELECT REC IF command, if any, 
present in this processing block. Using the specified keywords 
and sort id list, reading within the record type can be started 
and stopped at specified locations based upon values of the sort 
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variables. The sort variables must be specified in the sort i 
list in the order they were defined at schema definition time b 
the SORT RECORDS commands. 

FOR EACH REC 3 

FOR EACH REC 3, WITH(3,2) 

FOR EACH REC 3, FROM( VISTNUM . 1 1 THRU f VTSTN1IM n 1 


3, WITH(3,2) 

3, FROM(VISTNUM,1)THRU(VISTNUM,5) 
output specification 


The output specification list will be written at the top of es 
page after each page eject. The output specification 
described in the WRITE command. 


HOVE VAR LIST 


'CASE REPORT’ ,80X, 'PAGE.'.PAGE 
see SCHEMA RECORD DEFINITION COMMANDS 
(SEX EO 2 AND AGE GE 21) INDEX = INDEX+1 
variable list 


Values for the variables m the variable list are moved directlv 
v^7=m the i r e ? ord deln S processed into the summary record. 
Variable labels, value labels, missing values, etc. in effect or 
these variables are made available for SIR procedures which 
□ , 8 comraand raa y appear only in a PROCESS REC or FOR EACH 

nhL DlOCk. 

MOVE VAR LIST AGE, SEX, EDUCATN 


Causes page ejection in the retrieval output listing. If n is 
remain ln’the a page? JeCt ” U1 °° 0Ur ° nly if feWer than » lines 

,ESET variable listtassignment value list)[/var list 

(values)/...] 

val ? i f dles t0 specified values. The frequency with 
the variables are preset is determined by the location 
PRESET command in the retrieval set; 

—If located in the RETRIEVAL block, the variable 
listed are set at the beginning of the retrieva 
and unchanged by SIR throughout the remainder 
of the retrieval. 

—If located in a FOR EACH CASE block, the vari¬ 
ables listed are reinitialized at the beginning 
of each new case. 

^^ a PROCESS REC or FOR EACH REC block, the 
named variables are re-initialized once each 
time the process block begins execution. 

a SUMM ARY block, the variables initializec 
at the beginning of the retrieval are re-ini- 
tialized at the beginning of each new summary 
record. 

PRESET T0TAL1 TO TOTALIO (0)/SUM1, INDEX (0,1) 

headinv ESET Tf CO th^ nd mu ? t de located immediately after the block 
“ e iistf t^ assignment ^ 
assignment^value^list? ti3fied ' MISSINGisa valid member of the 


(sort id list)] [ 


(sort id list)! 


Ehe beginning of a processing block which initiates the 

before a FOR^iU^Fr in "f* Spei ?i fied r^ord type. If located 
R EACH REC command, or if none is present, the PROCESS 
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REC block Is executed once for each case. If located after a FOR 
EACH REC block, it is executed once for each summary record 
created within that case. 


PROCESS REC 
SELECT REC IF 
HOVE VAR LIST 
END 

PROCESS REC 


3 

(AGE GE 21 AND AGE LE NO) 
EDUCATE TO INCOME 

3,WITH(5,1) 


PROCESS REC 


2,AFTER(FAMMEMBR,VISTDATE,1) 
UNTIL(FAMMEMBR,VISTDATM,4) 


RECODE var listtvalue list = new value)...(value list= 

value) [var list...] if the recode does not 
change the recoded variable type or 

new variable=old variabletvalue list=new value) 
(value list=new value) [/new variable=old vari¬ 
able (CONVERT)...) if the recode results in a 
change in variable type. 


This command allows specifying new variable values and collapsing 
multiple values into single values. Since direct recoding of 
character to numeric values would be attempting to change a vari¬ 
able from character to numeric, it is necessary to recode charac¬ 
ter values into a new variable which is numeric in type as in¬ 
dicated by the second form of the command above. 

RECODE INCOME (0 THRU 5000=1) (5000 THRU 10000=2) 
(10000 THRU HIGHEST = 3) 


RECODE 


NEWCLAS=SOCCLAS('L0W’=1) ('MED'=2) C’HI'=3) 


REJECT CASE IF (logical expression) 

SELECT CASE IF (logical expression) 


These commands may appear in any SIR retrieval block except 
RETIEVAL and AFTER RETRIEVAL and they may be placed between 
blocks. If the SELECT (REJECT) CASE IF logical expression is 
false (true) then execution for that case stops and processing of 
the next case begins. 


REJECT REC IF (logical expression) 

Located within a PROCESS REC or a FOR EACH REC block, this 
command determines which records from the record type being read 
are to be processed by the commands within the block. Records 
are rejected if the logical expression is true. 

REJECT REC IF (SEX EQ 2) 


RETRIEVAL [AUTOSETt,UPDATE]] 

This command signals the beginning of a retrieval block and must 
be the first command in the retrieval command set. The AUTOSET 
keyword causes SIR to initialize the values of all summary 
variables in the summary record, that are not named in a PRESET 
command, to missing at the beginning of the retrieval and each 
time before a new summary record is begun. It AUTOSET is not 
used, SIR does no automatic summary variable initialization and 
values from previous summary records will carry forward unless 
initialized by the user. The UPDATE keyword allows modification 
of data in the data base by the retrieval. 

SELECT REC IF (logical expression) 

Located within a PROCESS REC or a FOR EACH REC block, this 
command determines which records from the record type being read 
are to be processed by the commands within the block. Records 
are accepted if the logical expression is true. 
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PROCESS REC 
SELECT REC IF 
COMPUTE 
END 


7 

(SEX Ed 1 AND EXISTS (SIZE) EQ 1) 
SUM 1 = SUM1 + SIZE ; CNT = CNT + 1 


SELECT FIRST (logical expression) 

Located in place of a SELECT REC IF, this command causes 
selection of only the first record for which the logical 
expression is true. 


SELECT FIRST (BLOODPRS GT 1R0) 

fINCLUDE I 

SUMMARY [ < L = variable list ] 

[exclude] 

This command signifies the beginning of a processing block which 
completes and passes the summary record to the procedures which 
follow. Final transformations may be performed on variables 
within the summary record or additional variables may be created 
by including necessary commands in this block. By default, 
unless multiple PERFORM PROC commands are included, a single copy 
of the present summary record is passed on to the procedures 
which follow each time this block is executed. By default, the 
summary record will contain every variable defined or moved by 
the retrieval program in the order defined. The EXCLUDE option 
allows indicating variables not to be included in the summary 
records passed to the procedures. INCLUDE specifies the set 
which is to be included and in what order. 


SUMMARY 

COMPUTE 

WRITE 

END 


EXCLUDE = TEMPI TO TEMPI 0, AVECOST 

AVECOST = COSTOTL/NVISTS 

'THE AVERAGE COST IS',AVECOSTtF8.2) 


VALUE LABELS 


variable list(value 1) label Kvalue 2)label 2. 
[/variable list (value 1)...) 


This command permits the user to specify labels for individual 
values of one or more summary variables. The command can be 
placed anywhere in the retrieval command set. 


VAR LABELS variable 1, label 1 [/variable 2, label 2...] 

This command allows the user to define variable labels for com¬ 
puted summary variables. 


VAR RANGES variable list (low,high) [/variable list...] 

This command allows the user to define a legal range for computed 
summary variables. Values falling outside the specified range 
will be assigned missing values. 

WRITE output specification 

WRITE IF (logical expression) output specification 

These commands provide a means of writing information to the 
output list device during a retrieval. Values may be printed in 
free format or under either format or edit specifications: 

Format specifications: 


Iw - Print the whole portion of an integer or floating point 
variable in a field w columns wide with no decimal point. 
•Fill if unable to fit. 

Fw.d-Print an integer or floating point variable right just¬ 
ified in a field w columns wide with d digits to the 
right of the decimal point. *Fill if unable to fit. 

Aw - Print a character string left justified within a field w 
columns wide. Right most characters dropped if unable to 
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Edit Specification: 


9 blank fill to the left 
Z Zero fill to the left 

decimal point in a column {one per edit) 

, comma in a column (more than one allowed) 

$ floating dollar sign (one per edit) 

* asterisk protected field (e.g. for check printing) 


The following example illustrate some uses of the edit 
specification: 


Value 


Edit Specification Output 


136.38 

('99999.9') 

136.38 

CZZZZZ.9' ) 

136.38 

C$9999.9') 

136.38 

('$ZZZZ.9') 

150 

CZZZZZ') 

9 7)) 2.50 

('999,999.99' 

136.38 

C$999,999.99 

136.38 

( i »*« ***.ggi 

97*12.50 

(i»** *»* _ ggi 


bbl 36.1) 
00136.1) 
b$ 136.)) 
$0136. i) 
00150 

b$9,7*)2.50 
bbb$136.38 
****136.38 
**9,71)2.50 


Line, page and column control: 


each / implies carriage return line feed. 

*at the start of the list causes page eject. 
nX causes spacing of n columns 


nT causes tabbing to column n 


Special functions allowed: 


DATE prints 
TIME prints 
PAGE prints 


present date, 
present time, 
present page of output. 


WRITE i)OT, 
WRITE IF 
WRITE 
WRITE IF 


'THIS REPORT WAS PRINTED ON', DATE, 'AT', TIME 
(SEX EQ 1) 'MALE RESULTS = ' MEN(F8.2) 

'THE ANSWERS ARE', 3X, VAR5(F8.2), 3X, VAR6(F8.2) 
(SEX EQ 2) ’FEMALE RESULTS =’,WMEN( 1 99999.99’) 


INTERACTIVE EDITOR C0MMANDS = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = 

The SIR Interactive Subsystem provides the facilities for performing 
interactive information retrieval from the user's data base. During a 
SIR Interactive run the user can: 

- input a set of SIR retrieval commands and correct and edit them 
with the interactive Editor, 

- execute the retrieval commands receiving the results at the 
terminal or directing them to local files, or offline printer, 

- save the retrieval commands in the SIR data base stored 
procedures file for later use, 

- execute previously defined sets of commands from the stored 
procedures file, 

- specify run-time parameters for a given retrieval, 

- have the SIR retrieval program interact with the user by 
prompting for and accepting input during execution, 

- change selected data base variables using COMPUTE and IF state¬ 
ments (provided the UPDATE option is included with the 
RETRIEVAL command). 
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The SIR Interactive Subsystem is invoked by executing SIR with the 
following command: 

SIR, IA. 

This command indicates that the SIR interactive version is to be 
executed. In addition to IA, the following parameters may be 
specified in the SIR command: 

FN = SIR file name 

PW = Password 

RS = read security 

WS = write security 

PS = Page size (lines) 

WD = page width (cols.) 

EP = named procedure to be executed immediately 

Required parameters will be prompted for by SIR if not provided in the 
SIR command. PS and WD will default to 60 and 136 respectively if not 
included in the command. The user may wish to set these parameters to 
a size appropriate for his terminal. If PS and WD are included in the 
SIR command without values specified, they default to 2N and 80 
respectively. See the SIR EXECUTION OPTIONS section for more details 
and information on the ENTER C/H/N/S/BLANK CR message. 

In the description of commands that follow, the comma and blank can be 
interchanged,i.e. either one can be used to separate the fields of an 
Editor command. No other characters are legal as field seperators. 
Conversely, any special characters can be used as text delimiters (on 
the SEARCH and REPLACE commands) except for comma and blank. 

Where line It s are referenced in a command, the * symbol may be used in 
place of a number to designate the last line in the present text. For 
example; 

L,» 

results in the last line of the current text file in the SIR Editor 
being listed. 

ACCEPT ,local file name (AC) 

Copies a local numbered text file, created previously by the 
editor or WRITE command, into the editor. RACCEPT rewinds the 
local file before copying. RRACCEPT rewinds before and after 
copying. 

ACCEPT, TEMPFIL 

ADD [.first line#[,increment]] (AD) 

This command informs the editor that new text is to be entered at 
the terminal. SIR will then prompt with line Its starting at the 
specified first line 0 and incrementing thereafter by the 
increment specified. Increment is 10 by default. If no line It 
is specified, SIR starts at 10 if the editor presently contains 
no text or the last existing line If plus 10 if there already is 
text. An '=' in response to a prompt informs the editor to ter¬ 
minate the add mode. Note: the editor will overwrite and inter¬ 
leave with existing lines. 

ADD, 10, 5 

Single lines of text may also be entered by typing a line number, 
a period and the desired line of text. For example, to insert a 
line between 70 and 80 of the current text, the user could type: 

>75.SELECT REC IF (WID*HGT GE 25) 

If the user specifies a line number that already exists in the 
text, the new line replaces the existing one. 

CLEAR (CL) 

Clears all current text from the editor work area. 
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COPY .source line, destination line (CO) 

COPY .source begin, source end, destination 

COPY .source begin, source end, destination, incr 

This command produces copies of text from one location at another 
location. The original text is left undesturbed. 

>COPY ,110, 370 
>C0PY ,270, 370, * 

>C0PY ,50, 290, 9000, 3 

In the first example, line 110 is copied into line 370. In the 
second example, lines 270 through 370 are copied after the the 
present end of the text. (when the * is specified, the 
destination line becomes the present last line plus the 

increment). In the third example, lines 50 through 290 are 

copied to lines 9000, 4003, 4006, etc. 

CWRITE .local file name (CW) 

Saves the current text without line numbers. A file written by 
this command is in card image format. The command RRCWRITE 
leaves the resulting file in rewound form. 

DELETE .line II (D) 

DELETE .first line II, last line II 

Removes the designated line(s) from the editor text. When the 

second form of this command is used, all lines in the interval 

are delete inclusive. 


DELETE, 214 
DEL, 181, 185 


FINISH 

END 

QUIT 


(F) 

(E) 

(Q) 


Terminates execution of the SIR Query Subsystem session. 


GET procedure name 


(G) 


Causes the text of the named procedure to be copied from the 
user's SIR procedure file into the editor work area. Any text 
presently in the work area is first cleared. 


GET, PATLIST 

HELP ( H ) 

HELP , command name 


Provides the user with helpful information on SIR Query commands. 
The first form lists a short description of every command while 
the second form lists a description of the specified command. 


LIST (L) 

LIST .line II 

LIST .first line II, last line II 

LIST .NAMES 


Lists all of the current editor work area text, a specified line 
of text, an interval of lines or a list of the names of the 
defined procedures currently in the user's SIR procedure file. 


LIST ,215 

LIST ,9,* 

L .NAMES 


move .source line, destination line (M) 

MOVE .source begin, source end, destination 

MOVE .source begin, source end, destination, incr 


This command moves one or more lines of text from one part of the 
users workspace to another. The original or source lines of text 
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are deleted. The arguments used with this command are identical 
to those used with the COPY command. 


OUTPUT ,local file name ( 0 ) 

Placed just before a RUN command, this command diverts all output 
from the execution to the named local file. Output defaults to 
the terminal again as soon as the execution resulting from the 
RUN command is completed. 

PAGE .number (PA) 

This command sets the page size for output at the terminal where 
the number given is the number of lines to be printed per page. 
The page size can be changed whenever the user wishes. Paging 
can be turned off completely by typing zero for the page size. 

After a page has been printed, the system types: 

ENTER C/H/N/S/BLANK CR 

A response of 'C' continues this listing without further paging 
'H' prints help on the command while 'N' or 'S' both terminate 
further printing. 'N’ completes the present execution without 
further output printing while 'S' terminates all processing 
immediately and returns control to the editor. A blank space 
followed by a carriage return tells the system to continue 
printing the next page. 

PREAD .file name (PR) 

This command reads a file containing one or more SIR procedures 
as previously written out by the PWRITE command. RRPR rewinds 
before and after reading. 

PURGE .procedure name list (PU) 

Deletes one or more procedures from the user's procedure file. 

PURGE, PATLIST .SORTST, REPORT! 

PWRITE .filename (PW) 

This command creates an external, non-SIR file containing a copy 
of all the procedures presently in the procedure file of the 
current data base. RPW rewinds the file before writing on it 
while RRPW rewinds the file both before and after writing on it. 

.local file name (REA) 

a local card image unnumbered file into the editor work 
Line numbers are automatically supplied beginning at 10 
and incremented by 10. If text already exists, the new lines are 
appended at the end and incremented by 10. 

READ, CARDFIL 

REPLACE ,d textl d text2 d (REP) 

REPLACE ,d textl d text2 d.line II 

REPLACE ,d textl d text2 d, first line II, last line II 

Replaces textl by text2 in all lines, in a specified line only 
or within a specified interval of lines. Any special character 
can be used as delimiters 'd'. 

The first form of the REPLACE command, which replaces text in all 
lines, should be used carefully. Since each line of text in which 
a change is made is listed at the terminal, it is suggested that 
the page size be set to 1 when using this command. In this way 
the results of the REPLACE command can be observed in the first 
line modified before instructing the system to continue with 
changes in all of the text lines, (see the PAGE command.) 


READ 

Copies 
area. 
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REPLACE «JONIS*JONES* 

REP 'THIS IS'THESE ARE',39 
REP/HALP/HELP/,10,400 

The three forms of the REPLACE command described above can also 
contain column anchors which make the search for the text! string 
more efficient. The general forms are: 

REPLACE ,C1 d textl d text2 d 

REPLACE ,C1,C2 d textl d text2 d 

In the first form, the replacement is only preformed on those 
lines which contain textl starting in column Cl. In the second 
form, the replacement occures if text! begins anywhere between Cl 
and C2 inclusive. 


RESEQUENCE [.starting line III, increment]] (RES) 


Resequences all the line numbers. If given without an argument, 
the text is resequenced by increments of 10 starting at 10. 


RESEQUENCE 385,10 

RETURN ,local file name 

Returns the named local file. 
REWIND .local file name 

Rewinds the named local file. 


(RET) 


(REW) 


RUN 

RUN 

RUN 

RUN 


(RU) 

.(parameter .parameter ,...) 

.procedure name 

.procedure name (parameter ,parameter ,...) 


Causes execution of current editor work area text or a named pro¬ 
cedure from the SIR procedure file. Parameters may be supplied 
as required by the procedure. If a procedure is run from the 
procedure file, present work area text is first cleared. Use of 
the OUTPUT command prior to the RUN command will divert the 
output from the execution to a local file. By default, output 
from the execution returns directly to the terminal. See the 
CALL command for details on the parameter list. 

RUN ,(1,32) 

RUN, PATLIST (1292,181)0) 

RUN 

SAVE .procedure name[ ,0K] (5A) 


Specifies that the current text is to be stored into the SIR 
procedure file from where it can later be retrieved for further 
modification or run as an SIR user defined procedure. The 

procedure name can be any name up to 7 characters in length 
starting with an alphabetic character. The OK keyword, if 
present, allows a present file with the same name to be 

overwritten. If OK is not present, the system will prompt for it 
before overwriting an existing file. 

SAVE .PATLIST 

SEARCH ,d text d (SE) 

SEARCH ,d text d, line II 

SEARCH ,d text d, first line II, last line II 

Searches for the specified text in all of the text, in the line 

indicated or within the interval given. Any special character 

can be used as delimeter 'd' except blank and commas. Each time 
a line is found containing the text specified, the line is 
printed at the terminal. 


SEARCH, *J0NIS* 

SE 'THIS IS' ,115,60 
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Column anchors can be used with any of the three forms of the 
SEARCH command. The general forms are: 

SEARCH ,C1 d text d 
SEARCH ,C1, C2 d text d 

If the first form of SEARCH is used, SIR will only list those 
lines in which the text starts in column Cl. If the second form 
is used, SIR will list the lines if they contain the text string 
anywhere between columns Cl and C2, inclusive. 

TAB .column 1 [.column 2,...] (T) 

Defines one or more tab columns to be used when lines of text are 
added. The tab character is the ampersand Whenever the tab 

character appears in a line of text being input, the remaining 
text in the line is placed starting at the next tab stop. The 
default tab stop is 16. 

TAB ,16, 28, MO 

WRITE .local file name (W) 

Saves the current numbered lines of text as a local temporary 
file. It can be read in again with the ACCEPT command. The 
command RRWRITE produces a rewound file. 

X .command 

The X command is designed primarily to perform the SIR utility 
routines and any SIR stand alone command that can be typed in a 
single line. The command is executed immediately without 
disturbing the text presently in the users workspace. 

X, LIST SCHEMA LABELS 3 

X, LIST STATS 

EP control card parameter 

The EP control card parameter allows immediate execution of a 
specified procedure upon entry into SIR. Using this parameter in 
the SIR control card as part of a computer system command file 
allows a user to execute SIR and enter a specified procedure by 
merely typing a single word. For example, the following control 
card would' cause SIR to immediately execute the procedure NEWP 
from the users procedure file: 

SIR, IA, PS, WD, EPrNEWP, FN=INFA, SP. 

PW=BABS, RS=RBABS, WS=WBABS. 

$ ,command. 

The $ command allows operating system commands to be executed 

from within SIR. 

$, FILES. 

Note: the $ command statement must be terminated by a period. 


SPECIAL RETRIEVAL COMMANADS============== = = = = = = = = = = = = = = = = = = = = = = = = = = = = = 

SIR provides the user with several extensions to the batch retrieval 
capabilities described under the section entitled BATCH RETRIEVAL 
COMMANDS. These commands are especially designed for use with the SIR 
interactive subsystem but,where appropriate, may also be used in 
batch. They provide: 

- a CALL facility for nesting procedures within other procedures 
together with a full parameter substitution (macro) capability, 

- two functions for providing interactive prompting for input 
from executing procedures and extensions to other retrieval 
commands to allow full use of these functions, 
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- commands to allow direct creation of a single case or a single 
record from within a retrieval. 


CALL procedure name!(parameter , parameter ...)] 

Allows calling a stored procedure from within another procedure. 
The text of the called procedure is placed within the text of the 
calling procedure. If a parameter list is included, the called 
procedure must be coded with a number in brackets in each 
location the respective parameter is to be substituted (numbered 
from the left in the parameter list). 

CALL CASENAME (I, CASE ID) 

The rules for forming the parameter list in the RUN and CALL 
commands are: 

1) The parameter list consists of numbers or character strings 
separated by commas. 

2) The entire parameter list must be enclosed in parenthesis. 

3) A list of items will be treated as a single parameter if it is 
enclosed in parenthesis. For example; (15, EXAM 1, (A1,A2,A3)) 
has 3 parameters; 15, EXAM1 and A1, A2, A3. 

N) If a parameter is enclosed in quotes, the quotes are considered 
as part of the parameter. 

5) Quoted strings containing commas must be in parenthesis to 
be considered a single parameter. 

6) Null parameters may be transmitted by two commas in succession 
or by empty left and. right parentheses. 

7) A parameter list cannot extend beyond one line and it cannot 
result in the line into which a parameter is substituted 
extending beyond its present maximum line length. 

CASE IS case id 

Causes the case with the specified id value to be created or 
retrieved in the data base. UPDATE is required on the RETRIEVAL 
command if a case is to be created. 

RETRIEVAL UPDATE 

COMPUTE ID = NREAD CENTER CASE NUMBER') 

CASE IS ID 

COMPUTE SEX = NREAD CENTER SEX(M = 1 F=2)') 

DELETE CASE 

Deletes the case presently retrieved. All records belonging to 
the case, including the CIR, are deleted. The user must have 
security clearance higher than that of any record or variable in 
the case and have used the UPDATE option in the RETRIEVAL 
command. 

DELETE RECORD 


Deletes the record currently retrieved from the data base. The 
user must have used the RETRIEVAL UPDATE command and have higher 
security clearance than that of the record or any variable in it. 

RETRIEVAL UPDATE 
PROCESS REC 2 

SELECT REC IF (BGNDATE LT JULN(4,1,73)) 

DELETE RECORD 
FINISH 

EXIT 

Causes termination of the procedure presently executing and 
returns control to the interactive editor. 


IF 


(logical expression) 


computation 
JUMP label 
EXIT 
WRITE 
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["delete CASE 7 
[delete recordJ 

Extends the batch IF command to allow any of the listed arguments 
in addition to computation. Multiple computations, separated b> 
are allowed as defined prevoiusly for the IF command. 

IFCCOUNT EQ 1) ANSWER=SPREAD('DO YOU WANT IT LISTED') 

IFCBIRTDATE EQ DATE) JUMP EARLY 
IFCCOUNT GT 50) EXIT 

JUMP label 

Causes execution to branch to the command with the specified 
label name. Labels can be applied to any command (other than a 
DO REPEAT or DO END command). A label is a 1-8 alphanumeric 
character name terminated by the symbol which precides the 

command being labeled. A label can appear in a line by itself or 
the convention LABEL : CONTINUE may be used. 

LABEL 1 : COMPUTE SEX = SREAD CENTER SEX') 

IF(SEX EQ 'MALE 1 OR SEX EQ ’FEMALE 1 ), JUMP LAB2 
PRINT ’SEX MUST BE MALE OR FEMALE’ 

JUMP LABEL1 
LAB2:- 

The following rules apply to use of the JUMP command for program 
branching control in a retrieval: 

a branch can always be made anywhere within a retrieval block. 

a branch can be made from within a block to outside of the 
block. 

- a branch cannot occur from outside a block directly into the 
block but only to the block header (e.g. PROCESS REC, FOR EACH 
REC, SUMMARY, etc.) 

NREAD (’prompt text’) 

A function which issues the specified prompt to the terminal 
along with a ’?’. A numerical response is then read from the 
terminal in free format and placed into the command from which 
the function was executed. 

10 COMPUTE VALUE = NREAD (’WHAT IS THE FIRST VALUE’) 
frt# i 

RECORD IS k >[,(sort idl, sort id2,...)] 

[ name J 

Identifies a specific record to be created or retrieved within 
the present case. If a record is to be created, RETRIEVAL UPDATE 
is required. 

RETRIEVAL UPDATE 
CASE IS 127 

RECORD IS 3,(1976,OCT) 

COMPUTE SALES = 1 A A.6;PROFIT = 9.5;UNITS = 1552 
FINISH 

In this example, a new record type 3 record is created in the 
case with id 127 and with sort variable values 1976 and OCT. 
This new record will contain the indicated values for SALES, 
PROFIT and UNITS. 

SREAD (’prompt text’) 

A function which issues the specified prompt to the terminal 
followed by a ’?’. A single character string response is then 
moved from the terminal into the command in which the function 
call is located. 
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10 LABEL 1 : COMPUTE NAME=SREAD('ENTER LAST NAME') 


SUPPLIED PROCEDURES============================ = = = = = = = = = = = = = = = = = = = = = = = 

The SIR procedures are an important part of the SIR retrieval 
capability. They enable the user to: 

obtain simple descriptive statistics on the summary variables 
(FREQUENCIES, PLOT, CONDESCRIPTIVE), 

interface directly to powerful statistical packages (SPSS SAVE 
FILE, BMDP SAVE FILE), 

- create a new SIR file from the summary records (SIR SAVE 
FILE), 

- create a printer or card image file from the summary records 
(WRITE RECORDS), 

produce a complex, hierarchical report from the summary 
records (REPORT). 


The procedure commands appear within the body of the user's retrieval 
command set between the SUMMARY and AFTER RETRIEVAL blocks. If the 
SUMMARY command has an EXCLUDE or INCLUDE option specified, only the 
remaining summary variables are available to precedures unless 
otherwise specified by a VARIABLE clause in the procedure. 

BMDP SAVE FILE FILENAME = filename/ 

VARIABLES = /ALL \ / 

[.variable list/ 

SORT = variable list / 

SAMPLE = proportion [,seed]/ 

BOOLEAN = (logical expression) / 

Creates a BMDP save file of the summary records. 

BMDP SAVE FILE FILENAME = BMDPRN/SAMPLE = .1 

VARIABLES = ^ ALL \ / 


CONDESCRIPTIVE 


ALL "l 
variable list] 
le name/ J 


WEIGHT = variat 
SAMPLE = fraction [.seed]/ 

TITLE = 'text'/ 

BOOLEAN = (logical expression)/ 
f ALL I 

STATISTICS = J keyword listj / 

Provides a listing of any or all of the following statistics: 


VAR....variance 
MAX....maximum value 
MIN....minimum value 
MEAN...the mean value 
STDE... standard error of 
the mean 

STDV...standard deviation 


SKEW....skewness 
KURT....kurtosis 

CV.coef of variability 

Cl.95% confidence interval 

about the mean 

WCOUNT..weighted count of non¬ 
missing values. 


CONDESCRIPTIVE 


FREQUENCIES 


VARIABLES = INCOME TO SOCSCR/ 
BOOLEAN = (GROUP EQ 1) 

STATISTICS = MAX, MIN, MEAN, STDV 


INTEGER = 


= {all X 

[variable listj 


(low,high) 


ALL 

[variable 

[variable list(low,high)...]/ 
GENERAL = /ALL 1 (/) of cells)[ 

[variable listj 

CONTINUOUS =/ALL X (# of cells, 

[variable listj low,high) 

(ALL 

INTERVALS = '(variable list 
WEIGHT = variable name/ 


[variable list...]/ 


t> (inter 
J [ . . . 


terval values) 

]/ 











SAMPLE = proportion/[,seed]/ 
BOOLEAN = (logical expression)/ 
STATISTICS =/ALL \ / 

I keywordsj 
TITLE = 'text' / 

HISTOGRAM / 


Provides frequency tables and optionally histograms on desire 
variables. The statistics option list available is that of th 
CONDESCRIPTIVE procedure with the addition of: 

MED . median 

MODE .mode 


FREQUENCIES 

FREQUENCIES 


GENERAL = SEXC2), SOCCLASS(6)/ 
HISTOGRAM 

CONTINUOUS = BLOODPRS (8,60,200)/ 
BOOLEAN = (ADULT EQ 1) 


PLOT 


SCATTERGRAM 

or 

LINEPLOT 


TITLE 

TITLEX 

TITLEY 

XDIV 

YDIV 

SYMBOLS 

SAMPLE 

BOOLEAN 


dependent variable list (low,high) WITH 
independent variable (low,high) BY control 
variable (low, high)/ 

dependent variable list (low, high) WITH 
independent variable (low,high) BY control 
variable Clow, high)/ 

'plot title'/ 

'x-axis title'/ 

'y-axis title'/ 

number of x-axis divisions/ 

number of y-axid divisions/ 

symbol list/ 

proportion[,seed]/ 

(logical expression)/ 


STATISTICS = ALL or NONE/ 

DEVICE = PRINTER (number of columns, number of lines 

lines per page)/ 


Provides two way scattergram or line plots of two variables wit 
the option of stratifying the plots on a third variable. A dif 
ferent symbol will be used for each stratum. Statistics printe 
when requested include the following results of a simple linea 
regression of the y-axis variable on the x-axis variable: 


correlation coefficient, r 

r-squared 

significance of r 

- standard error of the estimate (std dev of 
residuals) 

intercept, its std error and significance 
slope, its std. error and significance 

PLOT SCATTERGRAM = AGE(20,80) WITH SYSBP(60,180) 

BY SEX (1,2)/ 

TITLE = 'SYSTOLIC BP WITH AGE BY SEX'/ 

TITLEX = 'PATIENTS AGE'/TITLEY = 'SYSTOLIC BP'. 
SYMBOLS = M,F/ 

STATISTICS = ALL 


SIR SAVE FILE FILENAME = filename/ 

VARIABLES =/ALL \ / 

[var listj 

SAMPLE = proportion [,seed]/ 

BOOLEAN = (logical expression)/ 

RECTYPE = number [.name]/ 

Creates a new SIR file from the summary records. To allow chang¬ 
ing information in the new SIR file from that found in the sourct 
SIR file code book file, the following additional keywords art 
available to indicate desired characteristics of the new SII 
FILE. 







SIR - REPORT GENERATOR 


SPSS SAVE FILE 


NOFCASES 

RECSCASE 

CASEID 

MAXTYPES 

MAXKEYSZ 

MAXRECS 

SORTIDS 


INTEGER 


= number of expected cases/ 

= estimated average records per case/ 

= name of case id variable 
= maximum number of record types/ 

= maximum 9 of bits in the key/ 

= maximum records of any type/ 

= defines a list of sort id variables 
when there are to be two or more 
records per case. 

= variable list (min, max, scale), 
[variable list (min, max, scale)...] 
defines the list of variable to be 
stored as integers in the new SIR 
file. All variables not named in 
this command will be stored as real 
numbers. 


FILENAME = filename/ 


VARIABLES 

SUBFILES 

SORT 

WEIGHT 

SAMPLE 

BOOLEAN 


='variable list / 

= varname (valuel) namel(value2) 
name2.../ 

= variable [(D)][,variable...] 

= weight variable / 

= fraction [ ,seed] 

= (logical expression)/ 


Creates an SPSS save file from the summary records. 


SPSS SAVE FILE FILENAME = SPSSRUN/SUBFILES=SEX (1)MALE(2)FEMALE 

WRITE RECORDS FILENAME r name/ 

VARIABLES = variable list/ 

FORMAT = (format specification)/ 

SAMPLE = proportionf,seed] 

BOOLEAN = (logical expression)/ 

SORT = variable [(D)] [.variable [(D)],...]/ 
MISSCHAR = character 


Writes the summary records to the named file. The format specifi¬ 
cations may include: 


Floating point field, w columns wide, d columns to the 
right of the decimal point. Right justified. 

Integer field, w columns wide. Right justified. 

Character string field, w columns wide. Left justified. 

Skip n columns on the output record. 


Tab to column n before printing the next value. Tabbing 
may only be done to the right. 


Text in the FORMAT specification may be enclosed in quotes or be 
specified by a^hollerith designation (eg NHDATA). A slash in- 
of a new line. A repeat multiplier may be used 
expressions may be nested to 10 deep. By 


dicates the sta 
and parenthetical 
default, when a variable has a missing value, the specified field 
11 with blanks. The MISSCHAR specification provides a 
alternative character to be used in filling 


filled 

user specified 
missing value fields. 


REPORT GENERAT0R=== 


The SIR report generator provides the user of SIR with the tools need- 
ed to produce complex, hierarchical reports. It is a SIR procedure 
and works in a manner similar to the other SIR procedures, i.e. it 
operates on one summary record at a time, in the order they are pro¬ 
duced by the retrieval unless a SORT clause is specified. The REPORT 
procedure should be used when the report-generating facilities built 
.““ e S ? R retrieval are inadequate. Consequently, the report 
procedure will be needed whenever the report requires sorting, 
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multiple levels of headings, sub total and grand total computatic 
output, and the ability to print different kinds of output fc 
different records. 

The only data that may be used in a report are the summary recor 
variables and the variables defined within the report itself usir 
COMPUTE and IF statements. The user may not refer to CIR variables c 
any variables in the various records of the SIR data base directl 
from within a REPORT procedure. The users retrieval program, whic 
creates the summary records available to the report, should t 
designed so that any CIR variables or record variables required in th 
report are transferred to the summary record. 

The special variables PAGE, DATE and TIME are always available to th 
report procedure. 

In a report, we want certain actions to take place when the values o 
specific variables or logical conditions change. We call thes 
changes breakpoints. Usually the action we desire at breakpoints ar 
printing subtotals and grand totals. The structure of the report i 
determined by the ordering of the records from which the report is t 
be created on the breakpoint variables. The user decides in advanc 
which variables constitute the breakpoint variables and orders th 
summary records appropriately on these variables by use of the SOR 
clause in the REPORT command. 

REPORT SORT = variable!(D)][,variable[(D)]... ]/ 

BOOLEAN = (logical expression)/ 

FILENAME = local file name/ 

SAMPLE = proportion/! .seed]/ 

PAGESIZE = number of lines per page!.columns 
per line]/ 

The FILENAME clause is required, all other clauses are optional 
The default for the PAGESIZE command is 60, 136 respectively. Th 
number of columns per line must lie within the range 72-136 an 
the number of lines per page cannot exceed 66. The variable 
appearing in the BOOLEAN statement must be summary variable 
only. The SIR REPORT procedure command is followed by a comman 
set which produces the desired report. In addition to th 
following special report commands, HEADING, FOOTING, WRITE, WRIT 
IF, COMPUTE, BLANK LINE(S), PAGE EJECT and IF commands may b 
used as required. 

AFTER REPORT 

Begins a block of commands placed at the very end of the repor' 
procedure to be executed once after the report is complete. 

AFTER REPORT 

WRITE 'THE REPORT IS COMPLETE... THERE WERECOUNT(15),'CASES 
WRITTEN' 

END 

AT END 

Begins a set of commands to be executed each time a higher level 
sort variable changes in value and before execution is allowed tc 
return to the higher level. It is located after the LEVEL! 
heading and computation commands. It provides for printing sub¬ 
totals at breakpoints. If there are several AT END sets ir 
various levels of a branch, they are executed in reverse order. 
Summary variables should not be referenced in AT END sets direct¬ 
ly as they will contain the new value which generated the break¬ 
point. If it is necessary to print the value of a summary recorc 
variable as it exists proir to the breakpoint, it must be first 
saved in a new variable by a COMPUTE statement prior to the A1 
END set. The saved value can then be referenced in the at fmt 


2,FAMMEM 
MEMTOT=0 

AOT,'TOTAL 1 , 98T,COST('$9999.99') 


set. 

LEVEL 
COMPUTE 
AT END 
WRITE 
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SIR - REPORT GENERATOR 

BEFORE REPORT 

This block of commands may be used, immediately after the REPORT 
command to provide an introductory section to a report. This 
block is executed once at the very beginning of the report before 
any summary records are read. No variables from the summary 
record can be referenced in this block as they are not yet 
available to the REPORT procedure. BEFORE REPORT is designed for 
printing a cover page for the report and for initializing global 
variables for the report. 


Begins a set of commands for forming the body of the report. Each 
time a particular branch containing a DETAIL set is executed the 
commands in the DETAIL set are executed. A DETAIL block is 
usually located at the innermost level of each branch of a report 
although it may be located within any level. The DETAIL section 
will usually contain commands for printing individual lines of 
reports and for incrementing subtotals and grand totals. 


DETAIL 

COMPUTE 

WRITE 

COMPUTE 


VDATE=JULC(VISTDATE) 

6T,VDATE(A12),3 IT,PROCDR(IA),50T,C0ST(F6.2) 
GRTOT=GRTOT+COST;MEMTOT=MEMTOT+COST 


END 

END LEVEL 
END REPORT 


To improve readability of the user's REPORT procedure, an END 
command may appear at the end of BEFORE REPORT, AFTER REPORT 
blocks or at the end of LEVEL sections within the body of the 
report. The END cannot be used to separate the three sections of 
the LEVEL-headings, DETAIL and AT END. Optionally, the 
statement END LEVEL may be used to end a LEVEL section and END 
REPORT may by used to terminate the report set. 


EXIT 

EXIT IF (logical expression) 


Terminates the report. Can be used in the ON ERROR level to 
terminate the report in case of error. 

JUMP TO LEVEL level identification 

Directs execution from one branch path into another at a higher 
level where necessary report commands have already been defined. 
Allows re-joining separate paths to avoid multiple sets of du¬ 
plicate commands. 

LEVEL level number .breakpoint variable name 

(breakpoint logical con¬ 
dition) 


Identifies the levels in the nesting structure of the report. 
Nesting levels are indicated hierarchically by the level number. 
Levels are identified numerically 1 highest, 2 next...,k, lowest 
or innermost. Multiple branches at the same level are identified 
by decimal numbers m.n(eg. 2.1,2.3,2.4,etc.). Multiple branch 
modes in different branches which are identical may be labeled 
with up to 3 characters for identification(2.1 A,2.1B. 
2.1 ABC.etc .) . 


The breakpoint varible name must be the name of a summary record 
variable. A change in the value of this variable triggers a 
breakpoint action to occur at the given level. All breakpoints 
nested within the given level also take place. 


The breakpoint logical condition is a standard SIR logical ex¬ 
pression. A change in the value of this logical expression from 
false to true, or vice versa triggers a breakpoint at the given 
level and all levels nested within it. When a logical condition 
is used on a LEVEL command, the level number must be of the form 
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m.n to indicate one of multiple possible branches at the sami 
level within a report. 

Any computation or output commands which appear immediately after 
the LEVEL command and which are not contained within an AT END or 
a DETAIL command set will be executed only when a breakpoint 
occurs at that level or at a higher level. 

LEVEL 1,FAMID 

LEVEL 2,FAMEM 

LEVEL 3 . 1 ,(FAMEM EQ 1 ) 

LEVEL 3.2,(FAMEM EQ 2) 

LEVEL 3.2A,(STATUS LT V) 


The ON ERROR command provides a special LEVEL command which can 
be used in multiple branch reports where the LEVELS commands 
present in the report do not cover all possible conditions that 
could occur. If a condition occurs that is not covered by at 
least one of the branches, SIR would normally abort the report 
execution. The ON ERROR command will be executed instead if it 
is present. If a level number is included, it must be of the 
form m.n and applies only to errors occurring at level m. If no 
level number is included, it is a global error recovery block and 
applies to errors encountered at any level not covered by an ON 
ERROR block at a specified level. If used without a level 
number, this set can appear anywhere in the body of the report. 

ON ERROR 3.2 

”jj*TE 'ERROR AT LEVEL 3 . . . CASE ’ , FAMID(16 ) 

END 


1 PAGE level number 

Tells the system to which level to return for printing headings 
when a page eject occurs. When a PAGE EJECT command is en- 
™“ tered °*" a " automatic page eject condition occurs, and an ON 
PALE command has been included in the report, the following 
actions take place: 1 

- The page footing is printed if a FOOTING 
command was included. 

Page eject is performed. 

- The HEADING command, if any, is executed. 

- Starting at the level specified in the PAGE 
command, the current branch of the report is 
followed from the outermost level and all 
WRITE, WRITE IF and BLANK LINES statements 
which are not in DETAIL and AT END sections are 
executed. All other commands are ignored. 

- When the statement that caused the page eject 
is reached, the report processing resumes in 
the normal manner. 

If there is no ON PAGE command in the report, only the first 
three steps above will be taken before processing continues. 


SIR UTILITIES============= 


P rov ide capabilities for maintaining, archiving and 
SECURITY 8 ? at w l?ase ; . To run a utility, appropriate PASSWORD, 
stuuHiii (with data base administrator passwords) and GET FILE com¬ 
mands are required. 
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SIR - UTILITIES 


TAPE STORE FILENAME = local filename/NOREWIND/ 

NEWFN = newfile name/ 

NEWPW = new password/ 

This utility moves the entire data base tie. code book, detail, 
procedures and data file) to the named local file. This file 
will normally be a magnetic tape although there is no restric¬ 
tion on the physical location of the resulting file. If the 
NOREWIND option is specified, the file written will not be re¬ 
wound at the end. This utility does not purge the existing SIR 
data base. This utility is also used to unlock a locked data 
base. 

PURGE SIR FILE 


This utility deletes the SIR data base from the disk. It will 
typically be used after the data base has been stored on magnetic 
tape by use of the TAPE STORE utility. 


SIR MERGE 


FILENAME 

PASSWORD 

SECURITY 

RECTYPE 

IDTRUNC 

UPDATE 

BOOLEAN 

RENAME 


source SIR filename/ 
source SIR password/ 
source SIR read password/ 
source //, host II, host name/ 

YES or 
NO / 

ADD or 
REPLACE / 

(logical expression)/ 

(old var name, new name), [(...)]/ 


Merges a single record type of a source SIR file into a host SIR 
data base which may contain multiple record types. 


The first four clauses FILENAME, PASSWORD, SECURITY and RECTYPE 
are required and in the order shown. The SECURITY commands mu 3 t 
provide level 30 password access to the source SIR data base. The 
RECTYPE clause indicates the present record type number of the 
source record type and the number and name of the record type 
that it will be in the host SIR data base. IDTRUNC is used to 
truncate string case id's which may be longer than those allowed 
in the host SIR data base. UPDATE is used when the source record 
type is being merged with an existing record type in the host. 
If UPDATE is not specified, both adding new records in the host 
and overwriting old records is allowed. 


If the source record type is merged into an existing record type 
in the host, the following restrictions must be met: 

The source and host record types must have 
identical sort id lists. 

Each variable in the source record must have 
the same type and length as the host. Missing 
values defined in the host are used in all 
cases. 


VAR RANGES and VALID VALUES are ignored. 

- Only the VAR LABELS and VALUE LABELS of the 
host remain as the labels in the host after the 
merge. 

The RENAME clause allows renaming source variables to avoid dup¬ 
licate names with those already in other record types within the 
host SIR data base. 


The BOOLEAN can reference variables in the source SIR data base 
only. 


SIR SUBSET 


FILENAME 

RECTYPES 

BOOLEAN 

SAMPLE 


filename/ 

rt II or name, (logical exp.)[,rt// or 
name, rt// or name, (log . exp .),... ]/ 
(logical expression)/ 
fraction [ ,seed]/ 


AO 



or COUNT 
or LIST 


= total[,increment [.start]]/ 
= list of case identifiers/ 


Creates a new SIR file which is a subset of the original SII 
file. Specific record types or parts thereof can be extracted 
The BOOLEAN can apply to the CIR record only. The RECTYPE: 
clause specifies the list of record types that are to be include: 
in the new subset SIR data base and any logical filters that ari 
to be applied in deciding which records from the specified recorc 
types are to be included in the subset. If a logical expressioi 
is included for a record type, it can only reference variables ii 
that record type or the CIR. The COUNT, LIST and SAMPLE clause: 
are as defined for the FOR EACH CASE command. 


SIR FILE DUMP 


FILENAME 

SEOUENCE 

RECTYPES 


BOOLEAN 
SAMPLE 
or COUNT 
or LIST 

COMPVAR 


filename/[NOREWIND/] 
begin column, end column / 
rt II or name,(logieal expression) 
[,rt# or name.rt// or name (logical 
expression),...]/ 

(logical expression)/ 
fraction [ ,seed]/ 
total, increment, start/ 
list of case identifiers / 

NO 

YES / 


Produces a dump of the SIR data file to a card image file which 
is as closely as possible a reproduction of the data records 
orginally entered. An abreviated code book schema definition 
command set is first written to the specified file followed by 
the requested data recordsffor a complete schema command set, use 
the WRITE SCHEMA command). Computed variables created at date 
entry or by SIR MERGE will conform to the following rules wher 
written to cards: 


--Computed variables, if numerical, will be output on subsequent 
cards at the end of the record type to which they refer. They 
will be written in FlO.n format in columns 1-72. E conversion 
will be used as necessary if the number will not fit using 
FI .n. 


—Computed variables of character type will be output on subse¬ 
quent cards following the numerical computed values. 

--The specification COMPVAR = NO suppresses the output of 
computed variables and should be used if only a reproduction of 
the original data is desired. 


Note that any variables added to the SIR data base by an SIR 
MERGE are considered computed variables. 


SIR FILE LIST 


RECTYPES 



or name,(log.exp)[,rt# or 
or name,(log.exp.),...] 


name 


ORDER 

BOOLEAN 
COUNT 
or LIST 

or SAMPLE 



(logical expression)/ 
total!.increment!.start]]/ 
case id [ .case id ,...,[case id 
THRU case id ]] / 

fraction [ ,seed] / 


Produces a sequential listing, record by record of the data 
values present within each selected record within the requested 
record types of the desired cases. The ORDER clause controls the 
order in which variables are printed within each record. ALPHA 
implies alphabetical ordering and VARNUM implies variable number 
ordering. VARNUM is the default. 


SIR EXAMPLES===========: 


CODE BOOK SCHEMA DEFINITION EXAMPLE 






SIR - EXAMPLES 

This code book schema definition command set defines a two record 
type SIR data base made up of demographic and clinic visit 
records on patient cases. 


RUN NAME 
C 

FILE NAME 
PASSWORD 
READ SECURITY 
WRITE SECURITY 
C 

TASK NAME 
N OF CASES 
RECS PER CASE 
MAX REC TYPES 
CASE ID 
RECTYPE COLS 
SEQUENCE COL 
DOCUMENT 

COMMON LIST 
C 

TASK NAME 
RECORD SCHEMA 
DOCUMENT 

MAX REC COUNT 
VARIABLE LIST 

INPUT FORMAT 
DATE VAR LIST 
VAR LABELS 


VALUE LABELS 

VALID VALUES 
MISSING VALUES 
VAR SECURITY 
C 

TASK NAME 
RECORD SCHEMA 
DOCUMENT 

MAX REC COUNT 
SORT RECORDS 
DATA LIST 


DATE VAR LIST 
TIME VAR LIST 
SCALED VAR LIST 
VAR RANGES 
REJECT REC IF 
TASK NAME 
LIST SCHEMA 


EXAMPLE DATA BASE DEFINITION RUN 
***SYSTEM INITIALIZATION COMMANDS 
EXAMPLE 
AAAA 

(30)MASTERR,(15)USER1R,(55USER2R 
(39)MASTERW,(15)USER1W,(5)USER2W 
***CASE DEFINITION COMMANDS*** 

CASE STRUCTURE DEFINITION 
500 
150 
10 

PATNTNUM 

5,6 

80 

THIS EXAMPLE ILLUSTRATES COMMANDS FOR THE DEFINING 
OF A STUDY DATA SET CONSISTING OF TWO RECORD TYPES. 
PATNTNUM,BRTHDATE,SEX,VISTDATE 
***RECORD TYPE 1 DEFINITION SET*** 

RECORD TYPE 1 -DEMOGRAPHIC DATA- 

1, DEMOGRPH 

DATA FOR THIS RECORD IS COLLECTED AT THE PATIENTS 
FIRST CONTACT WITH THE CLINIC. 

1 

PATNTNUM,NAME,ADDRESS,BRTHDATE,DOCTOR,AGE,SEX, 
OCCUP.EDUCTN,INCOME 

(IN, 1X,A25,5X,A110,/ ,5X,A6,2I2,3I1,I6) 

BRTHDATEC MMDDYY) 

PATNTNUM ,THE PATIENTS ID NUMBER/ 

NAME ,THE PATIENTS NAME/ 

INCOME .FAMILY INCOME 

OCCUP (1)BLUE COLLAR C2)HOUSEWIFE (35WHITE COLLAR 
(N)OTHER 
SEX(1,2) 

NAME TO INCOME (BLANK) 

NAME, ADDRESS(30,305 

***REC0RD TYPE 2 - CLINIC VISIT DATA*** 

RECORD TYPE 2 - PATIENT CONTACT 

2, CONTACT 

THIS RECORD IS COLLECTED EACH TIME THE PATIENT MAKES 
A HEALTH CONTACT WITH THE CLINIC. 

300 

VISTDATE,VISTTIME 
FIXED(1)/1 

PATNTNUM 1-R (I),VISTDATE 5-10 (A) 

VISTTIME 11-111 (A) .WEIGHT 15-19 (I) 

HEIGHT 20-21 (I) 

VISTDATE(MMDDYY) 

VISTTIME(HHMM) 

WEIGHT (-2) 

WEIGHT(30.00,160.00)/HEIGHT(20,80) 

(VISTDATE LT BRTHDATE) 

CODE BOOK LISTING REQUESTS 
STRUCTURE, LABELS,DETAILED 


BATCH RETRIEVAL 


This retrieval command set produces an SPSS SAVE FILE made up of 
clinic contact data retrieved from the above SIR data base on 
patients over 21 years of age at entry who are overweight. 


RUN NAME 

GET FILE 

PASSWORD 

SECURITY 

RETRIEVAL 

PRESET 

C 

PROCESS REC 
C 

MOVE VAR LIST 


EXAMPLE SUMMARY RECORD RETRIEVAL 
EXAMPLE 

USER 1R, USER1W 
RECCNT(0) 

***READ THE DEMOGRAPHIC RECORD*** 
1 

***MOVE DESIRED VARIABLES*** 
BRTHDATE TO INCOME 


1(2 





END 

C 

REJECT CASE IF 
C 

FOR EACH REC 
C 

SELECT REC IF 

COMPDTE 

MOVE VAR LIST 

END 

SUMMARY 
SPSS SAVE FILE 
AFTER RETRIEVAL 
WRITE 

FINISH 


*»*DECIDE IF CASE IS DESIRED 
(AGE LE 21) 

***READ CONTACT RECORD*** 

2 

»**IF OVERWEIGHT,CREATE SUMMARY REC*** 

((SEX EQ 1 AND WEIGHT GT 85) OR 
(SEX EQ 2 AND WEIGHT GT 68)) 

RECCNT=RECCNT+1 
VISTDATE TO INCOME 

EXCLUDE RECCNT 
FILENAME=OVERWT 

'***SUMMARY COMPLETE',/,5X,RECCNT(I5)RECORDS 
WRITTEN 1 


SIR REPORT EXAMPLE 


Consider a data file made up of family clinic records. A report 
is desired which lists procedures done by Dr. Jones for his 
patient families, listing each precedure, the date it was done 
and its cost. The report is to be organized by family broken 
down by father, mother and children. A sample from the desired 
report is as follows: 


ANNUAL PATIENT COSTS - DOCTOR JONES 1/15/78 PAGE NO 1 

FAMILY NAME ANDREWS 


ADDRESS 1632 S 

28TH ST SEATTLE 



HUSBANDS NAME 

JOHN 



DATE 

PROCEDURES 


COST 

FEB 13,1977 

3906 


16.00 

DEC 1,1977 

1312 


6.20 



TOTAL 

$N7.60 

WIFES NAME 

SUE 



DATE 

PROCEDURES 


COST 

MAR 5,1977 

N218 


25.96 

OCT IN, 1977 

3762 


125.76 



TOTAL 

$151.72 


In the data, file, the CASE ID variable identifies individual 
familes and its name is FAMID. Within the families, individual 
members are identified by a sequential number (i.e.,father = 01, 
mother = 02, children = 11, 12, etc.) stored in the variable 
FAMMEM. Finally, each office visit is recorded by date in the 
variable VISTDATE, which is stored as a julian integer. We must 
first order the records by families (FAMID), then by family 
members (FAMEM) within cases and finally select only those visits 
for the year 1977 and for which JONES was the doctor. The fol¬ 
lowing SIR command set is sufficient to produce the desired 
report: 


RETRIEVAL 
FOR EACH REC 
SELECT REC IF 

MOVE VAR LIST 
PROCESS REC 
MOVE VAR LIST 
REPORT 

BEFORE REPORT 
HEADING 

LEVEL 

WRITE 

COMPUTE 
AT END 


VISIT 

(DOCTOR EQ 3 AND (VISDATE GT JULN 
(01,01,77) AND LT JULN (12,31,77))) 
VISTDATE, PROCDR, COST 
DEMOGRPH 

FAMID, FAMMEM, LNAME, FNAME, STREET, CITY 
SORT = FAMID,FAMMEM,VISTDATE/ 

FILENAME = REPORT1 

'ANNUAL PATIENT COSTS - DOCTOR JONES' 

NOT, DATE, 'PAGE NO..'.PAGE 
1,FAMID 

'FAMILY NAME',20T,LNAME(A20)/’ADDRESS',20T, 
STREET, CITY 

GRTOT = OjHEADING = 'CHILDREN' 


N3 




SIR -EXAMPLES 


WRITE 

LEVEL 

COMPUTE 

AT END 

WRITE 

LEVEL 

WRITE 

WRITE 

LEVEL 

ON PAGE 

DETAIL 

COMPUTE 

WRITE 

COMPUTE 

LEVEL 

WRITE 

JUMP TO LEVEL 

LEVEL 

IF 

COMPUTE 

WRITE 

WRITE 

JUMP TO LEVEL 
ON ERROR 
WRITE 


EXIT 

END REPORT 


20T, 'GRAND TOTAL 1 , 35T, GRTOT('$9999.99') 

2, FAMMEM 
MEMTOT = 0 

NOT, 'TOTAL' ,N8T,MEMTOT C$9999.99 ') 

3.1(FAMMEM EQ 1) 

5T 'HUSBANDS NAME', 22T, FNAME(A10) 

8T, 'DATE', 26T, 'PROCEDURES', 52T, 'COST' 

N 

2 

VDATE = JULC (VISTDATE) 

6T VDATECA12), 31T, PROCDR(IN), 50T, COSTCF6.2) 
GRTOT = GRTOT + COST; MEMTOT = MEMTOT + COST 

3.2, (FAMMEM EQ 2) 

5T, 'WIFES NAME', 22T, FNAME(AIO) 

N 

3.3, (FAMMEM GT 2) 

(HEADING EQ 'CHILDREN') WRITE 5T, HEADING(A8) 
HEADING = 1 ’ 

22T, FNAME(A10) 

8T, ’DATE', 26T, 'PROCEDURES', 51T, 'COST' 

1 | 

3. A 

'ERROR IN CASE', FAMIDQ6),/ 

'FAMILY MEMBER', FAMMEM(I2), 

'VISIT DATE', JULC(VISTDATE) 


INTERACTIVE RETRIEVAL 


Using the same SIR data file from which the REPORT example 
above was done, suppose we wish to be able to interactively 
request a list of all visits by date, procedure done and cost 
for any patient in the data base. We wish the program 
interaction and report to look like the following: 

RUN PROCLIST 


THIS PROGRAM LISTS A REPORT ON ALL 
PROCEDURES DONE-TO-DATE ON ANY 
PATIENT IDENTIFIED BY FAMILY ID 
AND MEMBER ID 

WHAT IS THE FAMILY ID? 178N63 
INDIVIDUAL ID? 01 


REPORT ON JOHN SMITH It 178M63 JAN. 1,1978 

ADDRESS 1832 1 NTH S.W. 

SEATTLE WN 
938N1 

PHONE (206)-8N2-1853 


DATE 


PROCEDURE COST 


FEB 13,1977 3906 

OCT IN,1977 4219 


16.00 

9.00 


TOTAL 25.00 


CONTINUE? YES 

WHAT IS THE FAMILY ID? 


The following SIR command set will provide the above interactive 
report when save as proclist in the SIR procedure file. 

10 RETRIEVAL 

20 WRITE IX,/."THIS PROGRAM LISTS A REPORT ON ALL",/, 

30 "PROCEDURES DONE-TO-DATE ON ANY",/, 



40 "PATIENT IDENTIFIED BY FAMILY ID",/, 

50 "AND MEMBER ID",/ 

60 LI : COMPUTE IDENT = NREAD ("WHAT IS THE FAMILY ID") 

70 ;INDID = NREAD ("INDIVIDUAL ID”) 

80 ' ;COSTT =0.0 

90 CASE IS IDENT 

100 PROCESS REC DEMOGRPH, WITH (INDID) 

110 WRITE /."REPORT ON",2X,FNAME, LNAME, 

120 2XIX, IDENT (17) , NOT , DATE 

130 WRITE 5T,"ADDRESS",15T,STREET,/, 

1 40 15T,CITY,STATE,/,15T,ZIP,/, 

150 5T,"PHONE",15T,PHONE,/ 

160 WRITE 4T,"DATE",20T,"PROCEDURE", 

170 35T, "COST",/ 

175 PROCESS REC VISIT .WITH UNDID) 

180 COMPUTE PDATE = JULC(VISDATE);COSTT = COSTT + COST 

190 WRITE PDATE(A 12), 22T,PROCDR(IN),33T, 

200 COST(F6.6) 

210 END 

220 WRITE 30T, '-',/, 25T,'TOTAL ' , 

230 33T,COSTT ('$999.99') 

240 COMPUTE TEST = SREAD ('CONTINUE') 

250 IF (TEST EQ 'YES’ OR 'Y’) JUMP LI 


SIR EXECUTION OPTIONS = = = = = = = = = = = =,= = = = = = = = = = = = - = = = = = = = = = = - = = = = = = = = = . = = = 
The SIR command has the general form: 

SIR, p1,p2,...,pN 

where the p's stand for parameters listed below. These parameters 
may appear in any order within the command. 

Some of the parameters have default and alternate default values. 
When the parameter does not appear in the list, the default value 
is used. When it appears with no value specified, the 
alternative default is used. When it appears in the form p1= 
value, the value given is used for the parameter. 


Parameter 

Type 

Default 

Alternal 

I = 

input file 

Batch 

INPUT 

COMPILE 

L = 

output file 

Batch 

OUTPUT 

OUTPUT 

IA 


Interactive 

_ 


FN = 

SIR file name 

Interactive 


' - 

PW = 

password 

Interactive 

_ 


RS = 

read security 

Interactive 



WS = 

write security 

Interactive 

_ 


PS = 

page length (lines) 

Interactive 

60 

24 

WD = 

page width (cols.) 

Interactive 

136 

80 

EP = 

procedure file name 

Interactive 




In the batch mode, 

the command 




SIR. 





is the same as: 





SIR, I=INPUT, L=OUTPUT. 




In the interactive 

EDITOR mode, the 

command: 



SIR, IA. 





is the same as: 

SIR, IA, PS=60, WD=136. 


The IA parameter must be used to obtain the SIR Editor sybsystem. 
It is suggested for terminal usage that the command: 

SIR, IA, PS, WD. 

be used to provide the alternate default values of 24 lines/ page 
and 80 columns/line, appropriate for many CRT terminals. 

When using SIR interactively, any required parameter (eg. PW, RS, 
WS, etc...) not appearing in the SIR execution command will be 
requested via an interactive prompt. Consequently, the more 
parameters the user includes in his execution command statement. 







SIR - EXECUTION OPTIONS 


the fewer introductory prompts he need contend with to get going 
in SIR. 

With this in mind, the use of an executable catologue procedure 
can reduce execution of SIR to simply typing a single user chosen 
word after completing logon procedures. For example, suppose the 
catalogue procedure looks like the following: 

GRABCSIR) 

SIR, IA, FN = FSTUDY, PS=10,WD=72,SP.AR, SP. 

PW = JON, RS=DBAR, WS=DBAW 

Provided the user has stored these lines as a catalogued 
procedure file called SIREX, he may bring up SIR on the study 
FSTUDY with the single command: 

-SIREX 

GOOD MORNING 
WELCOME TO SIR 1.1 
SIR EDITOR READY 


During SIR interactive runs, SIR prints the following "pause 
message" at the end of each page: 

ENTER C/H/N/S/BLANK CR > 

C response tells SIR to continue the listing without further 
pausing. It is equivelent to setting the page size to zero for 
this listing. Blank fallowed by carriage return will continue 
printing the next page after which it will again pause. The "N" 
and "S" specifications are useful when the output from a'SIR run 
(as opposed to an Editor listing) is being printed at the 
terminal. Both "N" and "S" responses terminate further output at 
the terminal and return control to the SIR Editor. However, "N" 
allows the SIR run to be completed (without further listing at 
the terminal), while "S" terminates the SIR run immediately with 
no further processing occuring whatever before control is 
returned to the editor. 

On most computer systems, the execution of SIR may be interupted 
from its present task and returned to the SIR Editor mode by 
simply typing the BREAK key on the keyboard. When the user types 
BREAK, the messages: 

SIR JOB REPRIEVE 

SIR EDITOR READY 


are typed on the terminal. All execution has been terminated and 
the user may proceed with another task. This capability should 
be used with caution and as a 'last resort' as it unconditionally 
terminates whatever is presently being done. It is useful, 
however, when the user finds he has begun something he does not 
wish to continue with. The user should check with his local site 
about alternatives to the BREAK key that might be in affect on 
his system. 

When additional files are to be used by SIR during a run (eg. ,if 
data is to be read), these files must be attached as local files 
under the name referenced for them in the SIR run before 
executing SIR (or the $ command may be used, in SIR, to attach 
them). 

The four SIR data base files (code book, detail, procedure and 
data) are saved automatically by SIR as permanent files when they 
are created during the SIR creation run. Any other files created 
by the SIR run (eg., an SPSS save file or an output list file) 
will exist after SIR has been terminated as local files under the 
name given to them in the SIR run. They must be properly 
disposed of (listed, saved, etc.) by the user prior to logoff or 
they will disapear. 
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SPSS GENERAL RULES 


The contents of the control field must begin in column 1 The 
control field is columns 1-15. The specification field is columns 
16-80. 

File names and variable names may have a maximum length of 8 
characters, the first of which must be alphabetic. They must be 
composed of letters and/or numbers. No special characters or 
embedded blanks may be used. 

The following words may not be used as variable names. 


If sequence numbers are present in columns 73-80 of the control 
statements, NUMBERED YES must be specified. 

The maximum number of variables per system file is 1000. Th^ 
maximum number variables per archive file is 5000. 

The maximum number of subfiles per file is 100. 

The HISSING VALUES statement can specify a maximum of 3 values per 
variable. K 

The ASSIGN MISSING statement can specify only one missing value per 
variable. ^ 

SUBFILE,SEQNUM, CASWGT values cannot be altered through COMPUTE. IF 
and RECODE statements or their starred counterparts. 

Valid formats on the INPUT FORMAT statement are: A, F, P, S, T, V, 

The N of CASES statement must specify: 

the exact number of cases for CARD input. If the exact number 
is not known use UNKNOWN option with an END OF RECORD card 
inserted at the end of the card input data. 

Quotation marks must be used with alphanumeric values in: 

a) the control cards. 

b) the data if the FREEFIELD form of the INPUT FORMAT control 
card is used. 

The maximum field width of an alphabetic variable is A10. However 
only the first 8 characters of an alphabetic value will be printed. 


SPSS, Statistical Package for the Social Sciences, Second" 
Edition, Nie, Hull, Jenkins, Steinbrenner and Bent, 
McGraw-Hill, 1975. 






SPSS - CONTROL COMMANDS 

SPSS CONTROL COMMANDS == 


CONTROL WORDS 


ADD CASES 
ADD DATA LIST 


ADD SUBFILES 
ADD VARIABLES 
ALLOCATE 
ASSIGN MISSING 
COMMENT 
COMPUTE 1 
‘COMPUTE / 
COUNT 1 
‘COUNT / 
DATA LIST 


DOCUMENT 
DELETE SUBFILES 
DELETE VARS 
DO REPEAT T 
END REPEAT/ 
ECOLOGY 

EDIT 

ERASE DOC 
FILE NAME 
FINISH 
GET ARCHIVE 


GET FILE 
IF 1 

•IF / 

INPUT FORMAT 


INPUT MEDIUM 


KEEP VARS 
LIST ARCHINFO 

LIST CASES 
LIST ERRORS 

LIST FILEINFO 


MERGE FILES 


MISSING VALUES 
N OF CASES 

NUMBERED 

OPTIONS 
OSIRIS VARS 
PAGESIZE 


PARAMETERS 


n 

FIXED(nrecs)/recnum varlist cols (format) 
(vartype) varlist... 
BINARY(nrecs)/recnum varlist varpos 
varlist.... 
subfile name(nl)... 
varlist 

TRANSPACE=number of bytes 
varlistCvalue) 
any text 

computed variable s 
arithmetic expression 
varnamescriterion varlist 
(value list) 

|FIXED[nrecs]/recnum varlist starting cols 
. [ending cols (vartype)] varlist... 
LBINARY(nrees)/recnum varlist varpos varlist 
any text 

subfile name list 
varlist 

stand-in varnamesvarlist/ 
any comment 
ON 
OFF 

check syntax only 
remove documentation 
(file name, file label) 
terminate SPSS run 
FILEsf1 , VARIABLES::varlist 


fl 


FILE= f 5 , VARIABLES:: ALL 
file name 

(logical expression) computed variables 
arithmetic expression 
FIXED(format list) 

<FREEFIELD 

Lbinary (format iist)J 
CARD I 

Jdisk L 
Jtape r 
LotherJ 

{ var list 
VARLIST 1 
COMPLETE/ 

CASES=n/VARIABLES=varlist 
to get a listing of SPSS error messages 
(30 printed pages) 

COMPLETE ^ 

VARLIST 
SORTVARS 
VARINFO 
LABELS 

SUBDIRECTORY 
DOCUMENTS 

FILES = f1 , VARIABLESsALL 


one or more 


FILE= f5,VARIABLESsALL 

S varlist(missing value list)/ 

unknown]" 

YES 1 
NO / 

option number list 

{ variable number list 
number of lines'l 
NOEJECT j 
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PRINTBACK 


PRINT FORMATS 
RAW OUTPUT UNIT 
READ INPUT DATA 
READ MATRIX 
RECODE 1 
•RECODE J 
REJECT IF 
RENAME VARS 

REORDER VARS 
RUN NAME 
RUN SUBFILES 


SAMPLE 

SAVE ARCHIVE 


Jno 1 

^CONTROLf 


(FORMAT J 

varlistCvalue)/... 
file number 


varlistCvalues=newval) 
varlist (CONVERT)/... 

(logical expression) 
varlist=namelist/ 
varlist=namelist... 
varlist 

run label of up to 6M c 
[(subfile list) (...)... 

•{each 
[all 

sampling factor 
[compressed] [file name,label/] 
FILE= fl, VARABLE=vrlist/ 


SAVE FILE 
SEED 


SELECT IF 
^SELECT IF 
SORT CASES 

STATISTICS 

SUBFILE LIST 
TASK NAME 
VALUE LABELS 
VAR LABELS 
VARIABLE LIST 
WEIGHT X 
*WEIGHT J 

WRITE CASES 


WRITE FILE INFO 


FILE= f5, VARIABLES=ALL 
[compressed] [file name,file label/] 
[STANDARD [ 

<a large odd integers 
[’alphabetic string'J 
(logical expression) 


rlist 


{<!i} 


Listl> 


(statistics number lis 
\.ALL 

subfile nameICnl) 
label 

varlistCvalue!)label 1. 
varnamel, var label/... 
variable list 
variable name 


'(format list)varlist 
OR 

BINARY varlist 

(options 1 with WRITE CASES ^ 
gives listwise deletion of 
..missing data.) _ 

VARIABLESrvarlist 


any one or comb, 
FILENAME 
VARLIST 
SUBFLIST 
NCASES 
[MISVALS 

charJbcd 1 

lEBCDj 

NUMBERED 


PRNTFMTS 

VARLABS 

VALLABS 

DOCUMENTS 

ALL 


/ 


SPSS STATISTICAL PROCEDURES 


AGGREGATE 


GROUPVARS; varlist/ 
VARIABLES; varlist/ 


AGGSTATS; 


any one or comb. 
SUM 

VALIDN 

MEAN 

SD 

MAX 

MIN 

SKEW 

KURT 

PCTGT (value) 
PCTLT (value) 
PCTBTN (low,high) 


C RMISS=value]/ 


of, 







SPSS - ANOVA 

Program to group cases on the basis of one or more 
aggregating variables. 

AGGREGATE GROUPVARS=CITYID/ 

VARIABLES=FREQUOTE,NCONTACT,CAMPAIGN/ 
AGGSTATS=VALIDN,MEAN,SD 

OPTIONS: 1.Includes missing data. 

2. Deletes missing data,1istwise. 

3. Grouping variables included in aggregated file, 
t.Outputs a compositional aggregated file. 

5. No end-of-record after each case in output 
binary file. 

6. An end-of-file to be written after the data on 
the binary output file. 

7.Suppresses the AGGREGATE output data file. 

8.Suppresses the printed output. 

9.Causes AGGREGATE output to be written on a 
formatted BCD output file. 


STATISTICS: 

1.Group identification number, total number 
cases and values of grouping variables. 

2.Statistics 1, plus exact contents of each 
in aggregate file. 

3.Statistics 2 for the first 10 cases only. 


of 

c ase 


ANOVA dependent varlist BY independent varlist(min,max) 

WITH covariate list/dependent varlist.../... 

One-to-five-way analysis of variance and covariance for 
factorial designs assuming a fixed-effects model. 


ANOVA INCOME BY EDUC,RACE,SEX(1,3) WITH AGE,FATHINC 

OPTIONS 8 

STATISTICS 1 


OPTIONS: 1. 

2 . 

3. 

N. 

5. 

6 . 

7. 

8 . 

9. 

10 . 


Includes missing data. 

Suppresses value labels. 

Ignores two-way and higher interactions 
among factors. 

Ignores three-way and higher interactions 
among factors. 

Ignores four-way and higher interactions 
among factors. 

Ignores five-way and higher interactions 
among factors. 

Processes covariates concurrently with main 
effects for nonraetric factors. 

Processes covariates after main effects for 
nonmetric factors. 

Indicates regression approach. 

Indicates heirarchical approach. 


STATISTICS: 

1. Multiple classification analysis.(availa¬ 
ble with option 9.) 

2. Unstandardized partial regression 
coefficients for the covariates. 


BREAKDOWN 


VARIABLES=varlist!low,high)/ 

TABLES=varlist BY varlist BY... (Integer mode) 

OR 

TABLES= varlist BY varlist BY... (General mode) 

OR 

VARIABLES= var1ist(low,high)/ 

CROSSBREAK=varlist BY varlist... (Crossbreak mode) 


Summary statistics among subgroups of the cases in 
the data file. 
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BREAKDOWN VARIABLES=ACHTEST(LO,HI) 
APTEST1 TO APTEST5 (0 , A ) 
MOTIVE/TABLES=ACHTEST 
BY APTEST1 TO APTEST5 
BY MOTIVE 


OPTIONS: 2 

STATISTICS: 1 


OPTIONS:1. 

2 . 

3. 

A. 


Includes missing data. 

Excludes missing data for dependent 
variables. 

Suppresses all labels. 

Prints in tree-diagram format. 


(CROSSBREAK Facility) 


1-3 Same as BREAKDOWN 
A. Not used. 

5. Deletes cell frequencies. 

6. Deletes the sums. 

7. Deletes the standard deviations. 

8. Deletes value labels but prints 
variables labels. 


STATISTICS: 1. One way analysis of variance table. 

2. Test of Linearity table. (1 also required) 

CROSSBREAK Facility 


1-2 Not used. 

3-12 Same as CROSSTABS 1-10 statistics. 


CANCORR VARIABLES=varlist/ 

RELATE=(n,F,t) varlist WITH varlist/ 
[CANVARsmdrp ] 

4 or W[ RELATES...] 

[CANVAR J 

Canonical Correlation analysis program. 


CANCOR VARIABLESsATTITUDI TO ATTITUD7 

BEHAV1 TO BEHAV6/RELATESATTITUD1 
TO ATTITUD7 WITH BEHAV1 TO BEHAV6/ 
CANVARs.30/ 


OPTIONS: 2,10 

STATISTICS: 3,A 


OPTIONS: 


1 . 
2 . 
3. 

A. 

5. 

6 . 


18. 

19. 

20 . 


Includes missing data. 

Deletes missing data, pairwise. 

Computes a covariance matrix instead 
of a correlation matrix. 

Matrix input. 

Computes matrix about the origin. 

Continues computations if either variable 
set is linearly dependent. 

Continues computations if a matrix is 
not positive definite, (with option 2) 

Outputs a simple correlation matrix file. 

Uses VARIABLE LIST with matrix input. 

(with option A) 

Calculates weighted estimates of missing 
data in computing canonical scores. 

(with option 2) 

Not used. 

Reads in mean and s.d. with correlation matri 
Output canonical and variate scores file. 
Substitutes means for missing data when corre 
lation or covariance matrix is being computed 


STATISTICS: 1. 

2 . 

3. 


Means and standard deviations. 
Correlation matrix. 

Correlation matrix only if non-compu- 
table coefficients are encountered. 


SPSS - CONDESCRIPTIVE 


4. Canonical variate scores. 



J varlist 1 

CONDESCRIPTIVE 

[ ALL J 

Univariate 

summary statistics on continuous variables. 

CONDESCRIPTIVE 

INCOME TO SESINDEX.SCALEI, 

SCALE3,ORGMEMBS,FAMSIZE 

OPTIONS: 

1 , 3 

STATISTICS: 

ALL 

OPTIONS: 1. 

Includes missing data. 

2. 

Suppresses variable labels. 

3. 

Outputs Z-score file. 

i). 

Prints a reference dictionary. 

STATISTICS: 1. 

Mean. 

2. 

Standard error. 

3-4. 

Not used. 

5. 

Standard deviation. 

6. 

Variance. 

7. 

Kurtosis. 

8. 

Skewness. 

9. 

Range. 

10. 

Minimum. 

11. 

Maximum. 

12. 

Sum. 

13. 

Coefficient of variation. 

14. 

95% C.I. for sample mean. 


CROSSTABS VARIABLES=varl1st(low,high)/ 


TABLES=varlist BY varlist BY... 

OR 

TABLES=varlist BY varlist BY... 

Two-way and 

n-way joint frequency tables program. 

CROSSTABS VARIABLESrRACE(0,1) ITEM1 TO ITEM10/ 

TABLES = RACE BY ITEM 1 TO ITEM10 

STATISTICS: 

ALL 

OPTIONS: 1. 

Includes missing data. 

2. 

Suppresses all labels. 

3. 

Deletes row percentages. 

4. 

Deletes column percentages. 

5 . 

Deletes total percentages. 

6. 

Deletes value labels but prints 
variable labels. 

7. 

Includes missing values in tables, 
not in statistics. 

8. 

Prints row variables in reverse order. 

9. 

Prints an Index of the tables produced. 

STATISTICS: 1. 

Chi-square / Fisher's exact test if cases 
<21/Yates' corrected Chi-square for all 
2x2 tables. 

2. 

Phi for 2x2 tables, Cramer's V for larger 
tables. 

3. 

Contingency coefficient. 

4. 

Lambda, symmetric and asymmetric. 

5 . 

Uncertainty coefficient, symetric and 
asymmetric. 

6. 

Kendall's tau b. 

7. 

Kendall's tau c. 

8. 

Gamma 

9. 

Somer's D, symmetric and asymmetric. 

10. 

Eta. 

11. 

Pearson's R. 
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DISCRIMINANT GROUPS: 


[varname 1 (high,low)/ 

[subfiles/ 


/TOLERANCE=value/ 


VARIABLES: varlist/ 

ANALYSIS: varlist [(level)].../ 

'DIRECT 
WILKS 
MAHAL 

METHOD MAXMINF 

MINRESID| 

RAO 

MAXSTEPS:value/FIN=value/FOUT:value/ 

PIN:value/POUT:value/FOUTsvalue/VINsvalue/ 
FUNCTIONS:ng,ip,sig/ 

Tequal 


PRIORSs < SIZE 

(_value list 

Discriminant analysis program 


RUN SUBFILES (A,B,C) (D,E) 

DISCRIMINANT GR0UPS=SUBFILES/VARIABLES:V1 TO V20/ 

ANALYSIS=V1 TO V20/METHOD=RAO 


OPTIONS: 10,13,16 

STATISTICS: 5 


OPTIONS: 1. 

2 . 

3. 

A. 

5. 

6 . 

7. 

8 . 

9. 

10 . 
11 . 

12 . 

13. 
1 A. 

15. 

16. 

17. 

18. 

19. 


STATISTICS: 1. 

2 . 

3. 

A. 

5. 

6 . 

7. 

8 . 

9. 


Includes missing data. 

Includes cases with missing values 
during classification. 

Suppresses step-by step output. 

Suppresses stepwise summary table. 

Prints classification results table. 
Prints discriminant scores and classi¬ 
fication information. 

Prints a single plot of cases. 

Prints a separate plot for each group. 
Suppresses classification phase features 
for cases initially unclassified. 

Prints territorial map. 

Prints unstandardized discriminant 
function coefficients. 

Prints classification functions. 

Rotates the discrimination functions. 

Uses individual group covariance matrices 
for classification. 

Outputs matrix materials file. 

Allows matrix input. 

Outputs discriminant scores file. 

Outputs file of membership probabilities 
for all groups. 

Outputs file of actual group and classi¬ 
fied group numbers. 

Means. 

Standard deviations. 

Polled within groups covariance matrix. 
Pooled within groups correlation matrix. 
Matrix of pairwise F ratios. 

Univariate F ratios. 

Test for equality of group covariance 
matrices. 

Group covariance matrices. 

Total covariance matrix. 


FACTOR 


VARIABLES:varlist 


[TYPE = 

DIAGONAL: 


' PA1 
PA2 
RAO 
ALPHA 
IMAGE 
BYPASS 
vl,v2,. ./ 


/ 
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NFACTORS= value / 
MINEIGEN; value/ 
ITERATE;value / 
STOPFACT=value/ 
fVARIMAX 
QUARTIMAX 
ROTATE= < EQUIMAX 
OBLIQUE 

[norotate 

DELTA= value / 
FACSCORE[=mdrp ]/...] 



Factor analysis program. 

FACTOR VARIABLES;VAR01 TO VAR30/ 

TYPE=PA1/ROTAT£=EQUIMAX/... 


OPTIONS: 2 

STATISTICS: ALL 

OPTIONS: 1. 

2 . 

3. 

A. 

5. 

6 . 

7 . 

8. 

9. 

10 . 

11 . 
12 . 
13. 
14-23. 
24. 


STATISTICS: 1. 

2 . 

3. 

4. 

5. 

6 . 

7. 

8 . 


Includes missing data. 

Deletes missing data, pairwise. 

Allows correlation matrix input. 

Allows factor matrix and communalities input. 
Outputs simple correlation matrix. 

Output factor matrix and communalities file. 
Outputs factor-score coefficient matrix file. 
Outputs means and standard deviations. 

Uses VARIABLE LIST with matrix input. 

(with option 3) 

Produces weighted factor scores when missing 
data is encountered, (with option 2 and 7) 
Sequence the output factor scores (with option7) 
Reads in mean, s.d. with correlation matrix. 
Output factor scores. 

NOT USED. 

Outputs fator scores in (4E20.14) format 
(unsequenced, i.e. option 11 not specified) 
or in (20X.3E20.14) format(sequenced,i.e. 
option 11 specified. 

Means and standard deviations. 

Correlation matrix. 

Inverse and determinant of correlation matrix. 
Communalities, eigenvalues and proportion of 
total and common variance. 

Initial-factor matrix. 

Rotated-factor matrix and transformation matrix. 
Factor-score coefficient matrix. 

Plot of rotated factors (with orthogonal 
rotations.) 


FREQUENCIES 


GENERAL; 



INTEGER;var1ist(low,high)[/varlist.../...] 


One-way frequency distributions with descriptive statistics. 


FREQUENCIES INTEGER;NHELP(0,7)INCOME TO EDUC INCOMER 

(0,99) 

GENERAL;INCOME,AGE,RACE,SEX,ITEM1 TO ITEM10 


OPTIONS: 


1. Includes missing data. 

2. Deletes value labels. 

3. Prints tables in 8-1/2 x 11 format. 

4. Outputs a permanent print file. 

5. Prints all tables in condensed format. 

6. Prints all tables greater than a page in 
condensed format. 

7- Deletes frequency tables. 



Prints histograms. 


STATISTICS: 


9. 

Prints a 

reference dictioi 

1. 

Mean. 


2. 

Standard 

error. 

3. 

Median. 


A. 

Mode. 


5. 

Stand ard 

deviation. 

6. 

Variance. 


7. 

Kurtosis. 


8. 

Skewness. 


9. 

Range. 


10. 

Minimum. 


11. 

Maximum. 


12. 

Sum. 


13. 

Coefficient of Variation. 

1 A. 

95J C.I. 

for sample mean. 

r variable name 

/ ] 


[CONSTANT 

[MATRIX matrix-identifier [=]value list /] 

[STORE store-identifier [=] value list /] 
[SYMMETRIC symmetric-identifier [ = ] value list /] 
[RESTRICT restrict-identifier[=] value 


store-identifier 

or [ 

matrix-identifier > 

or 

constraints-identifier 


| store-identifier 
i matrix-identifier 


/ 


[CONSTRAINTS constraints-identifier [ = ] constraints-list/1 
[HYPOTHESIS hypothesis-identifier[=] hypothesis-list/] 
[MIXED mixed-identifier [=] value 
[store-identifier 1 [store-identifier 1 fstore-identifier 


r l ima 


■{matrix-identifier ^ jmatrix-identifier {matrix-identifier 

FmiiTTnu aintS ^ den ^| fi fr ^4 r symmetric-identifier 

EQUATTON equation-identifier [=] variable name 
[WITH] variable list / 

rrSS T J, n ? m ? ] [=] E 0 UA TION [ = ] equation-list 
(IDENTITY [r] equation list j 

Undog [ = ] variable list EXOG[=] variable listj 
MIXED [=] mixed-list [RESTRICT [=] restrict-list] 

JFTEST^[=] restrict-list1 [N03SLS] [RESIDCOV] [COEFCOV] 

\rmonly J 

[CONSTRAINTS [=] constraints-list] 

[HYPOTHESIS [ = ] hypothesis-list] / 

Generalized and three-stage least squares. 


GUTTMAN SCALE scale narae=varname(division point)... 
Guttman scale analysis. 


GUTTMAN SCALE 
STATISTICS: ALL 


MEMBSCAL=ALTMEM(1) 
OCCMEMO) RECMEM(I) 


Includes missing data. 

Suppresses variable labels. 
Suppresses ordering of variables. 


STATISTICS: 1. 


Inter-item and part-whole correlation 
coefficient. 

Coefficient of reproducibility. 

Minimum marginal reproducibility. 
Percent improvement achieved by Guttman 
scale. 
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SPSS - JFACTOR 


5. Coefficient of scalability. 



OPTIONS: 1. 

2 . 


3. 

4 . 



Includes missing data. 

Deletes missing data,pairwise. 

Allows correlation matrix input. 

Allows factor matrix and communalities 
input. 

Outputs simple correlation matrix. 

Output factor matrix and communalities file. 
Outputs factor-score coefficient matrix file. 
Outputs means and standard deviations. 

Print normalized residual correlation matrix. 


STATISTICS: 1. 

2 . 

3. 

4. 

5. 

6 . 

7. 

8 . 

10 . 


Means and S.D. 

Correlation matrix. 

Inverse and determinant of correlation 
matrix. 

Communalities, eigenvalues and proportion 
of total and common variance. 

Initial-factor matrix. 

Rotated-factor matrix and transformation 
matrix. 

Factor-score coefficient matrix. 

Plot of rotated factorsCwith orthogonal 
rotations.) 

Print normalized residual correlation matrix. 





MANOVA <dep var list> BYCfactor list> WITH <covariate list> 
f POLYNOMIAL (metri 

TRANSFORM(variable list) = 1 REPEATED/ 


netric) / 

Lx) / J 


RENAME 

ANALYSIS 

PARTITION (factorname) 

CONTRASTt factorname) 


PROTATE 
PROTATE(n) 

PROTATE(MAXEIGEN=m) 


[^SPECIAL (matrix 
= (newnamel,newname2 ...)/ 

=<dep var list > WITH <covar list >/ 
= (df1,df2, ...)/ 


{ DEVIATION/ 
DIFFERENCE/ 

SIMPLE/ 

POLYNOMIAL(metric)/ 
SPECIAL(matrix)/ 


(significance level) / 


NOROTATE/. 

VARIMAX/ 

EQUIMAX/ 

QUARTIMAX/ 


( NOROTATE/ 

VARIMAX/ 

EQUIMAX/ 

QUARTIMAX/ 

= <list of effects in design>/ 

■iate Analysis of Variances 


MANOVA WEIGHT 1 BY SEX(1,2),DRUG(1,3)/ 

PARTITION(DRUG)=(1,1)/ 

CONTRAST(DRUG)=POLYNOMIAL(1,2,3)/ 

DESIGN =SEX,DRUG(1),DRUG(2),SEX BY DRUG(1) 
+ SEX BY DRUG(2) 


OPTIONS: 1. 


5. 

6 . 

7 . 

8 . 

9. 


10 . 


11 . 
12 . 
13. 
IN. 

15. 

16. 
17. 

STATISTICS: 1. 

2 . 

3. 
N . 
5. 


Include missing data. 

Print parameter estimates only. 

Obtain last effect in design as residual. 

Omit grand mean from model. 

Print estimates of parameters for summary 
tables only. 

Print factor's last parameter as "negative 
sum" . 

Input of with-cells variance-covariance 
matrix, cell means, and cell counts. 

Output of within-cells variance-covariance 
matrix, cell means, and cell counts. 

Input of within-cells correlation matrix, 
standard deviations, cell means, and cell 
counts. 

Output of within-cells correlation matrix, 
standard deviation, cell means, and cell 
counts. 

Output of error correlation matrices. 

Output of error variance-covariance matrices. 
Suppress special orthogonal processing. 
Perform computations using observations 
model. 

Calculate inner products in double precision. 
Suppress printing of parameter estimates. 
Print minimal multivariate output. 

Print observed cell means, deviations, and 
counts. 

Print variances-covariances and correlations 
of the estimated parameters. 

Print design matrix. 

Print error correlations. 

Print hypothesis sum-of-squares and cross- 
products . 
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SPSS - flULTIPLE RESPONSE 


6. Print step-down F-tests table. 

7. Print averaged F-test table. 

8. Print raw discriminant function coefficients. 

9. Print standardized discriminant function 
coefficients. 

10. Print estimates of effects for canonical 
variates. 

11. Print canonical correlations. 

12. Print tests for homogeneity of variances. 

13. Print test for homogeneity of variance- 
covariance matrices. 

11|. Print sum-of-squares and cross-products 
matrix for each cell. 

15. Print variance-covariance matrix for each 
cell. 

16. Print correlation matrix for each cell. 

17. Print principal components of error 
correlation. 

18. Print principal components of error 
covariance matrix. 


MULT RESPONSE GROUPS=group name [group label] (variable list 
(value list)) group name [group label] 
(variable list (value list)).../ 

VARIABLES::variable list (value list) variable 
list.. ./ 

FREQUENCIES::item list/ 

TABLESsitem list BY item list [BY item list..]/ 
item list BY item list .../ 


MULT RESPONSE GROUPS=PROBS NATIONAL PROBLEMS MENTIONED 
(PR0B1 PR0B2 PR0B3O ,9)) 

MAGROR READERSHIP OF NATIONAL MAGAZINES 
(TIME REPUBLIC NEWSWEEK STONE (2))/ 

VARIABLES = EDUC(1,3)/ 

FREQUENCIES = PROBS MAGROR/ 


OPTIONS: 


STATISTICS: 


NONLINEAR 


MODEL 

[DERIVATIVES 

PARAMETERS 


TABLES = PROBS BY EDUC MAGROR 

1. Include missing values. 

2. Causes the ’groupwise' deletion of missing 
values for multiple dechotomies. 

3. Causes 'groupwise’ deletion off missing 
values for multiple response groups. 

4. Causes suppression of value labels for other 
than multiple dichotomies. 

5. Causes the percents and totals in crosstab 
tables to be based on responses rather than 
on respondents. 

6. Use 8Q-column on-line output format. 

11. Forces a "three-up" condensed format for 
frequency tables involving other than a 
multiple dichotomy. 

12. Causes the condensed format for frequency 
tables to be used for all items other than 
multiple dichotomies. 

1. Include row percents in crosstab tables. 

2. Include column percents in crosstab tables. 

3. Include 2-way table total percents in cross¬ 
tab tables. 


VARIABLES =dependent variable list WITH 
indep. list, 

NB= no. of coefficients. 

YHAT= arithmatic expression 
G(I)= arithmatic expression ] 

[B(1)= value,... B(NB)=value] 

[BL(1)=value,... BL(NB)rvalue] 

[BU(1)=value.... BU(NB)rvalue] 



[BNAMEO)=value,. BNAME(NB)=value] 
[FIX(1)=value, ... FIX(NB)=value] 


{ MAXIT 
TOL1 
TOL2 
TOL 3 
TOLA 


=parameter value 


Nonlinear regression 


NONLINEAR 

MODEL 

DERIVATIVES 


PARAMETERS 


OPTIONS 

STATISTICS 


VARIABLES= y WITH x, NB=3 
YHAT= B(1) + B(2) / (x+B(3)) 
G (1) = 1. 

G(2)=1./(x+B(3)> 

G(3)= -B(2)/(x+B(3))**2 
BC1)=0. 

B(2)=10. 

B(3)=2. 

3 

4,7,8 


OPTIONS: 


1. Include missing data. 

2. Use Gauss' method in performing regression. 

3. Suppress derivative checking. 

4. Output formatted for 80-columns per line. 


STATISTICS: 


1. Print derivative checking report. 

2. Print singular values and vectors of Z Z 

at the initial values of the parameters and 
at the final parameter values. 

3- Print a detailed iteration description. 

4. Print iteration summary. 

5. Print approximate variance-covariance matrix 
of parameters. 

6. Print normalizing elements and correlation 
matrix of parameters. 

7. Print 95% C.I. about the parameters assuming 
a linear hypothesis. 

8. Print the results of exploration. 

9. Print estimated values of dependent variable 
and residual table. 


NONPAR CORR /varlist WITH varlistY 

\yarlist J 

Spearman and/or Kendall rank-order correlation 
coefficients. 


NONPAR CORR RESDYTH,INCOME,NACT,EDRESPON,OCLEURES 


OPTIONS: 


1. Includes missing data. 

2. Deletes missing data, listwise. 

3- Two-tailed test of significance. 

4. Outputs correlation matrix file. 

5. Kendall correlations. 

6. Kendall and Spearman correlations. 


STATISTICS: none, no additional statistics available. 


NPAR TESTS test name [parameter]=varlist [parameters] / test name. 


BINOMIAL = varlist 


("cutting point 1 

l or r 

Rvalue 1, value 2j 


CHI-SQUARE= varlist[(low,high)] / 
IEXPECTED= fi...fn/] 


(uniform! 

K-S ^ NORMAL > 
[POISSONJ 


[parameter]=varlist / 


59 





SPSS - ONEWAY 


RUNS 


( MEAN ‘ 
MEDIAN 
MODE 
value , 


MCNEMAR 

SIGN 

WILCOXON 

COCHRAN 

FRIEDMAN 

KENDALL 

MEDIAN 


=varlist 

=varlist 

=varlist 

=varlist/ 

=varlist/ 

=varlist/ 

=(value) = 


[WITH varlist]/ 
[WITH varlist]/ 
[WITH varlist]/ 


dep. varlist BY 
indep. variable 
(value 1,value 2) 



=dep. varlist BY indep. variable 
(value 1, value 2)/ 


Non-parametric statistics. 


NPAR TESTS 

K-S 

= INCOME BY SOCSCL(1,2)/ 

OPTIONS: 

1 . 

2. 

Include missing data. 

Listwise deletion of missing data. 


3. 

Special handling of variable lists for 
paired variable tests. 


4. 

Use random sampling if not all cases will 
fit into the maximum memory available. 


5. 

Output will be restricted to 75 printed 
columns. 


6. 

Suppress variable labels. 

STATISTICS: 

1 . 

Mean, maximum, minimum, SD, no. of valid 
cases for mentioned variables. 


ONEWAY dependent varlist BY independent variable 
(min,max)/ 

POLYNOMIAS=n/CONTRAST*n1,n2.. ./ 

RANGES* testname (alpha) where testname* 

One-way analysis of variance. 

ONEWAY DEPVAR BY INDEPVAR(1,6)/ 

CONTRAST*1 ,-.5,-.5/ 

POLYNOMIAL=3/ 

RANGE5=LSD(.10) 


LSD 
DUNCAN 
^ SNK > 
TUKEYB 
TUKEY 
LSDMOD 
SCHEFFE 


STATISTICS: ALL 


OPTIONS: 1. 

2 . 

3. 

4 . 

5. 

6 . 


Includes missing data. 

Deletes missing data, listwise. 

Suppresses variable labels. 

Outputs a file with number of cases, mean, 
and standard deviation for each category. 
Not used . 

Uses first 8 characters of value labels of 
independent variable as group labels. 

Uses input of category frequencies, means 
and standard deviations, not raw data. 

Uses input of category frequencies, means, 
pooled variance and degrees of freedom for 
the pooled variance. 

Uses harmonic mean for ALL group in range 
test. 


STATISTICS: 1. Number of cases for each category, means, 

standard deviations, standard errors, 
minimum, maximum and 955 confidence interval 
for the mean. 

2. Both fixed and random effects measures. 
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3. Homogeneity of variance statistics. 

't. Print random effects (MODEL II) estimate of 
the between component variance. 


PARTIAL CORR varlist WITH varlist BY varlist(order values)/... 

matrix type varlist BY varlisttorder values)/..!' 

Partial correlation coefficient program. 


PARTIEL CORR EDUCATN WITH INCOME BY FATHRED, FATHRINC(l) 


OPTIONS: 1. 

2 . 
3. 
A . 

5. 

6 . 

7. 

8. 
20 . 


Includes missing data. 

Deletes missing data, pairwise. 

Two-tailed test of significance. 

Allows matrix input. 

Outputs simple correlation matrix. 

Uses VARIABLE LIST with matrix input. 

(with option A) 

Omits degrees of freedom with significance. 
Prints only nonredundant partials. 
Substitutes mean for missing data when 
computing the correlation matrix. 


STATISTICS: 1. 

2 . 

3. 


Zero-order correlations. 

Means and standard deviations. 

Correlation matrix only if non-computable 
correlations are encountered. 


PEARSON CORR varlist WITH varlist 
varlist 

Pearson product moment correlations for pairs of 
variables. 


PEARSON CORR P0P60 WITH MEDSCH TO PTTERTRY 

OPTIONS: 3 

STATISTICS: ALL 


OPTIONS: 1. 

2 . 

3. 

R. 

5. 

6 . 
7. 

20 . 


Includes missing data. 

Deletes missing data, listwise. 

Two-tailed test of significance. 

Outputs correlation matrix file. 

Omits number of cases and significance. 
Prints only nonredundant coefficients. 
Compute correlation matrix using temporary 
means. 

Substitute mean for missing data when 
computing correlation matrix. 


STATISTICE: 


1. Means and standard deviations. 

2. Cross product deviations and covariance. 


REGRESSION fBACKWARDl 

[METHOD; 1 FORWARD }■ ] 

[stepwise J 

VARIABLES; varlist/ 

REGRESSION; dependent variable 

(NSTEPS, FIN, TOL, FOUT) WITH 
independent variables(inclusion level)/ 

[RESIDUALS=mdrp ] 


Stepwise multiple linear regression program. 


OPTIONS: 2,A 

STATISTICS: ALL 


OPTIONS: 1. 

2 . 

3. 

R. 


Includes missing data. 

Deletes missing data, pairwise. 
Suppresses variable labels. 
Allows matrix input. 
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SPSS - RELIABILITY 


5. 

6 . 

7. 

8. 


9. 


10 . 
11 . 
12 . 


1 


3. 


19. 

20 . 

21 . 


Means and standard deviations with matrix 
input. (with option 9) 

Suppresses step-by-step output. 

Suppresses the summary table. 

Outputs correlation matrix file. 

Uses VARIABLE LIST with matrix input. 

(with option M). 

Sequence residuals file. 

Output of Residuals. 

Write y' and residuals on alternative 
output file. 

Suppress the plotting of cases with missing 
data . 

Suppresses printing of axes on plots of 
residual vs. y' 

Output of means and standard deviations. 
Variance/covariance matrix output. 

Output normalized regression coefficients 
(Bs) . 

Output unnormalized regression coefficients 
(Bs) . 

Force regression through the origin. 
Substitute mean for missing data. 

Print t, not F in step-by-step output. 


STATISTICS: 1. 

2 . 

3. 

it. 


6 . 

7. 


9. 

10 . 


Correlation matrices. 

Means, standard deviations and number of 
cases. 

Correlation matrix only if non-computable 
correlations are encountered. 

Plot residuals against sequence of eases in 
a file. 

Statistics for residuals Von Neumann Ratio, 
Durbin-Watson statistics and a runs test are 
computed. 

Plot standardized residuals against 
standardized y' values. 

Print correlation matrix with number of cases 
for each coefficient. 

Print estimated variance-covariance matrix 
of the unnormalized regression coefficients. 
Compute 95$ C.I. for unnormalized regression 
coefficient (B). 

Compute 95$ C.I. for normalized regression 
coefficients. (BETA). 


RELIABILITY 


VARIABLES= varlist/ 
SCALE(label) = scale list / 
[ALPHA I 

. SPLIT | 

M0DELS= YGUTTMAN l 

STRICTPARALLELf 
[PARALLEL J 


FORMATS = INPUT(i/j), 0UTPUT(k/l) 
VARIABLESrVarlist /... 


Item and scale analysis program 


OPTIONS: 


6 . 


9. 

10 . 

11 . 
12 . 


Includes missing data. 

Not used. 

Suppresses labels. 

Allows input of covariance matrix. 

Allows input of standard deviations 
and correlation matrix. 

Allows input matrix to be triangular. 
Allows input of means. 

Outputs covariance matrix file. 

Uses VARIABLE LIST with correlation 
matrix input. 

Stops after writing a covariance matrix on 
BCD out. 

Outputs item means on BCD out. 

Outputs triangular matrices as single 
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vectors. 

13. Outputs matrices as single vectors 
rather than by rows. 

I 1 !. Space saver solution. 

15. Automatic deletion of items with zero 
variance. 

16. Continue processing with zero variance 
iterns. 

17. Friedman's X . 

18. Cochran's 0. 


STATISTICS: 1. 

2 . 

3. 

11 . 

5. 

6 . 


9. 

10 . 


Item means, standard deviations. 
Covariance matrix. 

Correlation matrix. 

Scale mean, scale variance. 

Mean, SD, minimum, maximum and range 
item means. 

Mean, SD, minimum, maximum and range 
item variance. 

Mean, SD, minimum, maximum and range 
interitem covariance. 

Mean, SD, minimum, maximum and range 
interitem correlation. 

Prints total statistics. 

Analysis of variance. 

Tuckey test for additivity. 

Hotelling T-squared. 


of 

of 

of 

of 


SCATTERGRAM varlisttlow,high) 

varlisttLOWEST,high) 

varlist (low,HIGHEST) 

varlist(low,high) WITH varlist(l,h) 

Two-way scattergram plots of pairs of variables. 


SCATTERGRAM 


VAR003(0,10) WITH VAR006 (7,HIGHEST) 
VAR007 TO VAR009(LOWEST,99) 

OPTIONS: 

1 . 

Includes missing data. 


2 . 

Deletes missing data, listwise. 


3. 

Suppresses variable iabels. 


R . 

Suppresses plot grid lines. 


5. 

Prints diagonal grids. 


6 . 

A two-tailed test of significance, 
(with statistic 3) 


7 . 

Automatic scaling. 


8. 

If not enough core is available, plot: 
of the first n cases are produced. 

STATISTICS: 

1 . 

Pearson r. 


2. 

r-squared. 


3 . 

Significance of r. 


R . 

Standard error of the estimate. 


5. 

Intercept with the vertical axis. 


6 . 

Slope. 


SPECTRAL 


fUNIVARIATE | 


=varlist / 


WINDOW 


PARZEN 
<BARTLETT > 
TUKEY 
HAMMING 
UNIT 


[(lags)] / 


LAGS = <lag> / 

<unit> [(period)] 

Spectral Analysis of Time Series 
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SPSS - SUmARY TABLES 


SPECTRAL BIVARIATE= VAR1,VAR2 / 

WINDOW = HAMMING(10) / 

OPTIONS 4,5,18 

STATISTICS ALL 


OPTIONS: 1. 

2 . 

3. 

4. 

5. 

6 . 


11 . 
12 . 

13. 

14. 

15. 

16. 

17. 

18. 


Ignore missing value indicators. 

Listwise replacement of missing observations 
by the mean. 

Suppress variable labels. 

Plot auto-correlation and spectrum for each 
variable. 

Plot unweighted cross-correlation, phase and 
coherence. 

Plot weighted cross-correlation, phase and 
coherence. 

Not used. 

Not used. 

Write unweighted auto-covariances to 
alternate output file. 

Write weighted auto-covariances to alternate 
output file. 

Write spectrum auto-covariances to alternate 
output file. 

Write unweighted cross-covariances to 
alternate output file. 

Write weighted cross-covariances to alternate 
output file. 

Write co-spectrum and quadrature spectrum to 
alternate output file. 

Write phase and coherence to alternate output 
file. 

Write phase and amplitude to alternate output 
file. 

Write gains to alternate output file. 

Use the unbiased estimator of the covariance, 
rather than the positive-definite estimator. 


Only one "write" option may be selected from the 
univariate set, 9, 10, 11; similarly, only one from 
the bivariate set, 12 through 17. 

Each statistic requested is produced for each window- 
lag pair. 


STATISTICS: 1. 


For each variable, at each lag/band j: 

j=0, 1,.. . ,m 

auto-covariance, 

weighted covariance, 

auto-correlation, 

weighted correlation, 

power spectrum, 

spectral density, 

loglO of the power spectrum, 

period of the frequency band. 

For each variable: 

standardized bandwidth of the window, 
equivalent degrees of freedom, 
variance ratio 
variance. 

For each bivariate pair(x,y), at each 
lag/band j: 

cross covariance of x with lagged y, 
weighted covariance of x with lagged y, 
cross covariance of y with lagged x, 
weighted covariance of y with lagged x, 
co-spectrum, 
quadrature spectrum. 

Not used. 

For each bivariate pair, (x,y), at each 
band j: 

gain from x to y, 
gain from y to x, 
phase shift, 
amplitude, 
coherence, 
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period of the band, 
frequency of the band. 


SUMMARY TABLES dependent variable list [(label )] BY 

independent varlist (min,max) 

[BY control variable list (min , max BY...]/ 

Combines several crosstabulation of dichotomous variables 
into one single table. 

SUMMARY TABLES RULE 1 TO RULE9 BY SEX(1,2) 


OPTIONS: 


1. Ignore missing data indicator. 

2. Listwise deletion of missing data. 

3. Suppress printing of variable labels for 
dependent variables. 

N. Suppress printing of value labels for 
independent and control variables. 

5. Suppress printing of both variable labels and 
value labels for independent and control 
variables . 

6 . Print the count of non-zero cases in each 
cell. 

7. Cause percentages to be computed within rows 
instead of within categories of independent 
variable. 


STATISTICS: 


no statistics are available. 


SURVIVAL survival time variable BY group variable (lo,hi)/ 

INTERVALS = THRU upper limit BY interval width / 
[PLOT = (group), group2,...groupn ).../] 
[COMPARE = (groupl, group2, groupn) .../] 

Survival Analysis 

SURVIVAL ADMSURV BY RACE(1,3)/ 

INTERVALS; THRU 98 BY 2/ 

PLOT =0,2,3)/ C0MPARE=(1 ,2,3) 


OPTIONS: 1. 

2 . 

3- 


Ignores missing data indicator. 

Causes pairwise comparison of each of the 
groups named in a COMPARE specification. 
Causes only comparisons between groups to 
be performed. 


STATISTICS: 


there are currently no optional statistics 
available with SURVIVAL. 


TETRACHORIC VAR TABLES;varlist (low,high)/ 

CORRELATIONS;varlist WITH varlist/ 


Tetrachoric correlation coefficients 


TETRACHORIC VARIABLES=VAROO1 TO VAR009(1,2)/ 
CORRELATION=VAR001 TO 
VAR008 WITH VAR009/ 


OPTIONS: 


1. Includes missing data. 

2. Deletes missing data listwise. 

3. Suppresses variable labels. 

4. Invokes special handling for paired variable 
lists. 

5. Causes Hotelling's T**2 statistic for the 
equality of two vectors of means to be 
computed and printed. 


STATISTICS: 


none. No additional statistics available. 


SPSS EXAMPLES: 
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SPSS - EXAMPLES 

These examples are oriented to the CDC NOS operating 
system. 


An example of an SPSS job in which the data is on 
punched cards and only one procedure (viz. "FREQUENCIES") 
operates on that data: 


GRAB(SPSS) 

SPSS. 

7/8/9 
RUN NAME 
VARIABLE LIST 
INPUT FORMAT 
RECODE 

VALUE LABELS 

MISSING VALUES 
INPUT MEDIUM 
N OF CASES 
FREQUENCIES 
OPTIONS 
STATISTICS 
READ INPUT DATA 
20F685 
A0M266 


(END-OF-RECORD CARD) 

COMPLETE DECK FOR RUN FROM CARDS 
AGE,SEX,RACE,INCOME,EDUCATN 
FIXED(5X,F2.0,A1,3F1.0) 

SEX('M'=1) ( 1 F'=2) CELSE=0) 

SEX (1) MALE (2) FEMALE/ 

RACE (1) WHITES (2) BLACK (3) INDIANS 
AGE,RACE TO EDUCATN (BLANKJ/SEX (0) 
CARD 
10 

GENERAL=ALL 

1,4 

ALL 


(10 DATA CARDS IN ALL) 


FINISH 

6/7/8/9 (END-OF-FILE CARD) 


An example of an SPSS job in which the data is on 
punched cards where the number of cases is unknown 
and more than one procedure (eg. "CONDESCRIPTIVE", 
& "CROSSTABS") is used on that data: 


GRAB(SPSS) 

SPSS. 

7/8/9 
RUN NAME 
FILE NAME 
DATA LIST 

INPUT MEDIUM 
N OF CASES 
CONDESCRIPTIVE 
OPTIONS 
STATISTICS 
READ INPUT DATA 


(END-OF-RECORD CARD) 

COMPLETE DECK FOR RUN FROM CARDS 
EXAMPLE 

FIXED(2)/1 AGE6-7/2 SEX 8(A), 

RACE,INCOME,EDUCATN 9-11 

CARD 

UNKNOWN 

ALL 

3,4 

ALL 


7/8/9 

COMMENT 

*REC0DE 

CROSSTABS 

STATISTICS 

FINISH 

6/7/8/9 


(UNKNOWN NUMBER OF CASES) 

(END-OF-RECORD) 

END-OF-RECORD CARD TO SIGNAL END OF DATA 
SEX(’M' = 1) CF'=2) (ELSE = 0) 

TABLES=AGE BY SEX TO EDUCATN 
ALL 

(END-OF-FILE) 


An example of an SPSS job in which the input data is 
read from the permanent disc file "SPSSDATA": 


GET(DISK,SPSSDATA) 
GRAB(SPSS) 

SPSS. 

7/3/9 
RUN NAME 
FILE NAME 
VARIABLE LIST 
INPUT FORMAT 


(END-OF-RECORD) 

DATA ON DISC RATHER THAN ON CARDS 
EXAMPLE 

VAR01 TO VAR05 
FREEFIELD 
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COMMENT 
COMMENT 
COMMENT 
COMMENT 
VAR LABELS 

VALUE LABELS 
INPUT MEDIUM 
N OF CASES 
LIST CASES 
COMMENT 
COMMENT 

READ INPUT DATA 
SCATTERGRAM 
FINISH 
6/7/8/9 


'FREEFIELD' MEANS THE NUMBER MAY 
BE IN ANY FORMAT AS LONG AS THEY ARE 
SEPARATED BY BLANKS AND ORDERED WITH 
RESPECT TO VARIABLE LIST. 

VAR01,AGE/VAR02,SEX/VAR03,RACE/VAR0N 
INCOME/VAR05 EDUCATION 
VAR02 (1) MALE (2) FEMALE 
DISK 
100 

CASES=10/VARIABLES=VAR01 TO VAR05 
TO LIST FIRST 10 CASES OF THE 
SUBFILE GROUP 

VAR01 TO VAR05 

(END-OF-FILE) 


An example of an. SPSS job in which a "save" file 
(i.e. SPSS control cards and data) is created and 
stored in a permanent disc file "MYDATA": 


GRAB(SPSS) 

SAVE(SVFILE,MYDATA) 

7/8/9 (END-OF-RECORD) 

RUN NAME SAVE FILE RUN 

FILE NAME EXAMPLE 

VARIABLE LIST AGE,SEX,RACE,INCOME,EDUCATN 


(OTHER SPSS COMMANDS) 


READ INPUT DATA 


(DATA DECK) 

COMMENT ANOTHER SPSS PROCEDURE AFTER THE DATA 

T-TEST GROUPS = SEX( 1) / VARIABLES::INCOME 

SAVE FILE 
FINISH 

6/7/8Z9 (END-OF-FILE) 


An example of an SPSS job in which a "SAVE" file is 
retrieved from the permanent disc file and used as the 
data from a "CROSSTABS" procedure. 


GET(GTFILE,MYDATA) 


COMMENT 
COMMENT 
GRAB(SPSS) 
SPSS. 

7/8/9 
RUN NAME 
GET FILE 
COMMENT 
COMMENT 
COMMENT 
COMMENT 
CROSSTABS 


THE LOGICAL FILE NAME ON THE GET CARD MUST BE 
"GTFILE" 


(END-OF-RECORD) 

USING A SAVE FILE 
EXAMPLE 

THE PARAMETER IN THE SPECIFICATION FIELD OF THE 
"GET FILE" CARD MUST BE THE SAME AS THE PARAMETER 
USED IN THE SPECIFICATION FIELD OF THE FILE NAME 
CARD WHEN THE FILE WAS CREATED. 

AGE BY SEX BY RACE 


(OTHER SPSS COMMANDS) 

FINISH 

6/7/8/9 -(END-OF-FILE) 
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BHDP - BIOMEDICAL COMPUTER PROGRAMS 

BMDP INTRODUCTION = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = 

The BMDP computer programs are general purpose and 
designed to handle a wide variety of statistical problems. 
To adapt the programs to a specific problem, parameters 
must be specified that describe the characteristics of the 
data to be analyzed, the mode of analysis desired and 
choice of output. 

Important problem parameters include: 

- problem title (identifies the problem) 

- description of the input data (format, size and 
source) 


- description of the variables (names and limits) 
transformations of variables 


- modes of analysis (specific to the various problems) 

- choice of output (numerical results and plots) 

Variables can be given an alphanumeric name up to 8 
characters in length. A name can contain embedded blanks 
if the name is specified in quotes. 


Continuous variables can be categorized into intervals 
specifying CUTPOINTS. Catgories can be given names 
the GROUP paragraph. 


by 

in 


Data within BMDP programs on which statistical analysis is 
done are always processed as decimal values. 

Format . statements for BMDP use the FORTRAN F, X, A, / and 
repetition specifiers only. 


BMDP CONTROL . . ... 

THE BMDP control language is placed on input cards or 
records in free format, the only restrictions being that 
the control information for each new problem or subproblem 
must start on a new 80 column card or record. 

PARAMETER NAME 

A parameter name identifies a control language paragraph 
or sentence. Specific parameter names have defined 
meanings as given in the program descriptions. 

/PROBLEM TITLE='HEART STUDY'. 

PROBLEM and TITLE are parameter names. 

SUBSCRIPTING 

If parameters need to be used more than once in a 
paragraph in reference to different variables, subscripts 
may be used to identify the variable being referred to. 

CUTPOINTS(5) = 1,3,7. 

CUTPOINTS(2) = 5,9,15. 

These two sentences refer to variables 5 and 2 respect- 
tively. 


*BMDP-77, Biomedical Computer Programs, P-Series, Edited by W j" 
and M.B.Brown, University of California Press, Berkeley, 1977.’ 


Dixon 


68 







SENTENCES 


Each control language command or parameter is spedified in 
a sentence which must end with a period. Most sentences 
assign values. 

CASES = 213. 

FORMAT = • (5F7.0) 1 . 


Sentences can be directive. 

REWIND. 

NO REWIND. 

NO CORRELATION. 


PARAGRAPHS 


Sentences are grouped into paragraphs. The name is its 
first word (eg. INPUT, PROBLEM, SAVE, etc.). A '/' must 
appear either before the paragraph name or at the end of a 
paragraph. Although each type of paragraph allows only 
specific sentences, the order of these allowable sentences 
in the paragraph is arbitrary (except the PROBLEM 
paragraph which must appear first). 

/PROBLEM TITLE = 'HEART STUDY'. 

/INPUT VARIABLES ARE 5. 

CASES = 213. 

FORMAT IS ' (5F7.0) ' . 

/VARIABLE NAMES ARE 'WT IN 72', HEIGHT, AGE, 'PULSE RT' SEX 


. 213 data cards 
/FINISH. 

RESERVED WORDS 

IS, ARE, TO and BY are reserved words and can only be used 
as variable names if placed in quotes. 

SPECIAL FEATURES 

BLANKS 

Blanks nan be used freely to space program control 
information, except in the middle of names or numbers. 

EQUAL SYMBOL, IS or ARE 

The equal symbol '=' and the words 'is' and 'are' can be 
used interchangeably. 

CASES = 213. 

FORMAT IS '(5F7.0)'. 


REPETITION 

Repetition of the same numerical values in a list can be 
specified by the symbol '*'. 

1i2,3,3,3,V,6,6,6,6 is the same as 1,2,3*3,R,R*6 

IMPLIED LIST 

To facilitate specification of a list of numbers which are 
equally spaced, the procedure: 

beginning It TO ending It BY step 
may be used. 


For example the list: 

1,2,3,R,5,5.5,6,6.5,7 

May be given by 
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BMDP - MAIN CHART 


1 TO 5 BY 1,5.5 TO 7 BY .5 

TABBING 

A number in parenthesis preceding a value can be used as a 
subscript of the list element that follows 

NAMES = AGE, SEX, (15) WEIGHT, PULSE, 
give names to variables 1,2,15 and 16. 

INPUT and FORMAT 
are the same as 
INPT and FRMT 

PARAGRAPH NAME REPETITION 

The name of a paragraph can be repeated before each 
sentence if desired. 

COMPOUND ASSIGNMENT SENTENCES 

Elements in lists cf two or more assignment sentences in 
the same paragraph can be stated in a form to facilitate 
matching. For example: 

VARIABLE NAMES = AGE, WEIGHT, SEX, CHANGE. 

MINIMUMS = 0,0,1,-10. 

MAXIMUMS = 20,275,2,10'. 

can be written: 


VARIABLE NAMES, MINIMUMS, MAXIMUMS = 
AGE, 0 , 20, 

WEIGHT, 0 , 275, 

SEX, 1 , 2, 

CHANGE, -10 , 10. 


repetition, implied list, subscripting or tabbing 
cannot be used with this feature. 

REPEAT DECLARATIONS 

The same location in a list can be specified more than 
once, but only the last is retained. 

CONTRAST = 4*1, -7, 3*1. 

can be written 

CONTRAST = 8*1, (5), -7. 


MAIN CHART PARAGRAPH AND PARAMETER DECLARATIONS 

The following table lists the paragraph and parameter 
names which are common to most of the BMDP programs. In 
the individual program descriptions, those names which 
have a function in that program identical to that 
described in this table will be mentioned in name only and 
the user can refer to this table for comments. If an 
option is specific to a program, it will be explained in 
that program. Default values are listed in parenthesis at 
the end of comments. 

PARAGRAPH PARAMETER 

NAME NAME 


/PROB 

(required) 


TITLE='text'. 


COMMENTS 


Signifies the beginning of a new 
problem and is required to begin 
each problem. 

Text in quotes to describe the 
problem. Max. length 160 charact¬ 
ers. (blank) 
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/INPUT 

(required on 
first problem) 


VARIAB=0. 


FORMAT ='C' 


UNIT 


II. 


CONTENT = C. 


LABEL 


•C' . 


TYPE = name. 


Input describes the input data siz; 
format and source. Required for th 
first problem of a deck setup. 

States the number of variable; 
required unless input is from a sav 
file, (none/prev.) 

VARIABLES ARE 10. 

Specifies read input format. Musi 
be less than 800 characters long 
Required unless input is a savi 
file, (none/prev.) 

FORMAT is •(5F7.0.3F10.2)' . 

Specifies the number of cases to bi 
read. By default, read stops whei 
EOF is reached. (EOF/prev.) 

CASES = 300. 

Specifies logical unit number foi 
data or save file. Must be statec 
unless the data is with the control 
cards. Unit number is 5 by default, 
Cannot be 1,2, or 6. Should be les: 
than 10. (5/prev.) 

UNIT is 9. 

Specifies REWIND to rewind the input 
file. NO REWIND suppresses rewind. 
Input unit will rewind by default 
after input data is read. 
(YES/prev.) 

NO REWIND. 

Required if input is from a save 
file. The code name (1-8 char¬ 
acters) must correspond to the code 
given to the save file when origin¬ 
ally stored. (input controlled by 
format) 

CODE IS MYFILE. 

Specifies what type of save file is 
desired (eg. DATA, CORR, MEAN OR 
FREQ). If not stated, the first 
subfile whose code matches the one 
specified will be read, (none) 

CONTENT IS CORR. 

Allows identifying desired input if 
the same CONTENT (ie. DATA) is 
recorded more than once under the 
same code. A label can be up to 40 
characters long. (none) 

LABEL IS 1 COMPUTED IN BMDP4M'. 

Specifies the way the program 
should treat the input. Acceptable 
types vary from program to program. 
(DATA) 

TYPE = CORR. 

TYPE = FREQ. 
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BMDP - MAIN CHART 

/VARIAB 

(optional if input from 
BMDP save file) 

ADD = II. 


NAME = list. 


USE = list. 


LABEL = list. 


MISS = lllist. 


GROUP = V. 

MAX = #list. 
MIN = #list. 

BEFORET.or 

AFTERT. 


/GROUP or 
/CATEG 

(optional) 


CUTP (II ) = 

in, 112, _ 


CODE (II) = 

in, ii 2 . 


Specifies variable characteristics. 
This paragraph is not required if 
all defaults are acceptable. 

Specifies the number of variables to 
be added by transformations. Can be 
negative if the transformation 
reduces the number of 

variablesCcounting from the right). 
(0 /prev.) 

Assigns 1-8 character names to the 
variables. Names must be enclosed 

in single quotes if they contain 

embedded blanks.Names are optional. 

NAMES ARE HEIGHT, ’WT LBS'. 

Specifies which variables are actu¬ 
ally going to be used in this pro¬ 
blem. Optional. (all variables) 

USE = 1,3, 12 TO 16. 

Specifies variables (up to two var¬ 
iables of 4 characters each) which 
are to act as case labels, (none 
/prev.) 

LABELS are 4,5. 

Missing value list for the vari¬ 
ables. (none/prev.) 

MISSING ARE 5*999. 

Specifies a variable on thich the 
data will be group.(none/prev.) 

GROUPING is 3. 


Maximum value list.(none/prev.) 

Minimum value list.(none/prev.) 

MAX = 35, 12, (15)8. 

MIN = 0, 0, (15)1. 

Determines whether data checking 
(MISS, MAX & MIN) is to be done 
before or after transformations. 
(BEFORET/prev.) 


Specifies group or category cut- 
points and associated name for 
variables. This paragraph is 

optional. If any outpoints, codes 
or names are given, all from a 
previous subproblem are eliminated. 
If none are stated, any from a 
previous subproblem are still in 
effect. 


Outpoints for grouping or categor¬ 
izing a continuous variable. n-1 
outpoints determine n categories. 
(CODE/prev.) 


For categorical variable, specifies 
the values which are to act as 
grouping values for the data set. If 
no codes are given for a grouping 
variable, outpoints will be used. 
If neither given, the 10 smallest 
values are used. 
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/TRANSF 

(optional) 


NAME (S) = A 1-8 character name for outpoints 

C1,C2,- or codes. If not stated, outpoint 

or code values will be used for 
names when neccessary. 

CUTPOINTS (5) ARE 2, 1|. 

CODES (II) = 5,7,3 10. 

NAMES (5) = BABY, TODDLER, 
OLDER. 

NAMES (N) = A, B, C, D. 

RESET - RESET, all assignments in pre¬ 

vious GROUP paragraph are reset to 
preassigned values, (not reset) 

Provides capability for minimal 
computation. A FORTRAN subroutine 
capability exists for more exten¬ 
sive transformations. If paragraph 
consists only of TRANSF, the trans¬ 
formations used in the preceding 
problem will still be in effect. If 
any new transformations are defined, 
all old ones are eliminated! 
Variable may be referenced for 
transformation by name or by index. 

Up to 7 separated temporary var¬ 
iables may be used . 

Allowed transformations 

- Arithmetic 

simple assignment X(1)=X(10). 
binary operation X(1)=X(1)+X(2). 
functions X(1)=SQRT(X(3)). 

- Logical X(7)=X(1) LT X(2).,where 

X(7) is 1(true) if X(1) is less 
than X(2) and is O(false) 
otherwise. 

- Random number generation 

X(1) = RNDU{star ting value). 

Allowable forms 

- result = argument 

- result = argl operator arg2 

- result = function (argument) 

- where result can be: 

X(i) for i an integer 
a variable name 
a temporary variable 
a word USE 

- and the argument can be: 

X(i) for i an integer 
a variable name 
a previously defined 
temporary variable 
the words USE.KASE or XMIS 
a numerical constant 

- USE can be used to select 
cases for analysis. 


USE = X(1) LT X(2). 

- uses those cases for which 
X(1) is less than X(2). 


73 



BMDP - TRANSFORMATIONS 


/SAVE 

(optional) 

UNIT = #. 

NEW. 

CODE = C. 

CONTENT = 
name. 

LABEL = 

'text'. 


/END 

(required) 

/FINISH 


- KASE is the ease sequence //. 

X(5 ) = KASE. 

- places a sequential number in 
variable 5 which numbers the 
cases. 

- XMIS is the internal missing 
value code. 

TEMP = X(1) EQ 7. 

X(1 ) = XMIS IF TEMP. 

- states that the first variable 
is assigned the missing value 
whenever it has a value of 7. 

Causes information (data, statis¬ 
tical results, etc.) to be written 
to a save file. (Not required) 

Logical unit to which the save file 
is to be written. Cannot be 1, 2, 5 
or 6. Should be less than 10. 
(none) 

Required if this is a new save file, 
(not new) 

A 1-8 character label for the save 
file, (none) 

Specifies what is to be written (ie. 
DATA, FREQ, CORR, MEAN, etc). Check 
individual programs for options 
available. (DATA) 

Up to a 40 character label to 
identify specifically the save file, 
(blank) 

/SAVE NEW. 

UNIT = 9. 

CONTENT IS DATA. 

CODE IS WERNER. 

LABEL = 'WERNER DATA AS OF 
3-13-74’. 

Required as last paragraph. 


Terminates BMDP. Last card in deck. 
BMDP programs dated before 
August,1977, require the FINISH 
paragraph when the data are not on 
cards with the command deck or when 
the number of cases is specified. 
The FINISH paragraph must not be 
used when the number of cases is not 
specified and data is included on 
cards. The FINISH paragraph is not 
neccessary in BMDP programs dated 
August 1977 or later. 


TRANSFORMATION SUMMARY TABLE 

KIND OF HOW RESULT TYPICAL 
TRANSFOR- USED COMMENTS EXAMPLE 
MATIONS 
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X C1 >. 


Arithmetic a=b 
(simple 
assignment) 

Arithmetic LOG(a) 

C functions) 

SQRT(a) 
EXP(a) 
ABS(a) 

SIN(a) 
COS(a) 
ATAN(a) 

INK a) 


SIGN C a) 


NONB(a) 


Arithmetic a+b 
(binary) 

a-b 

a*b 

a/b 

a**b 

aMODb 


Logical a LE b 

Operation a LT b 

a GE b 

a GT b 

a NE b 

a EQ b 

a AND b 
a OR b 


Random RNDU(i) 

Number 

Generation 

RNDG(i) 


b placed into a 


loglO(a) 

square root of a 
to the power of a 
I at,absolute 
value of a 
trig sine of a 
trig cosine of a 
trig arc tangent 
of a 

integer part of a 

-1 if a<0 
0 if a=0 
1 if a>0 

blank field 
detection returns 
missing value if 
a is blank 

addition 

subtraction 

multiplication 

division 

exponentiation 

the remainder 
of a divided 
by b 

1 if a<b,0 if a>b 
1 if a<b,0 if a>..b 
1 if a>„b,0 if a<b 
1 if a>b,0 if a*.b 
1 if ayb ,0 if a/b 
1 if a=b,0 if ayb 
1 if ayo and byo 
a if byO, other¬ 
wise leave as is 

uniform (0,1) with 
starting value i. 

normal (0,1) with 
Starting value i. 


X (7) = 


WEIGHT = LOG(WEIGHT). 

X(2) = SQRT(X(5)). 

X(1 ) = EXP(X(2)). 
CHANGE = ABS(X(17)). 

X(1) = SIN(X(1)). 

X(2) = COS(X(2) ) . 

X(3) = ATAN(X(3)) . 

AGE = INT(AGE). 


X(5) = SIGN(X(1)). 


X (7 ) = N0NB(X(7)). 


X(3) = X(1)+X(2). 

WEIGHTS=WEIGHT-X(3). 
DIST=RATE*TIME. 

X(1 ) = 1.0/TIME. 

X(1 ) = HEIGHT*» 3- 
M0NTH=AGE MOD 12. 


USE = X(1) LE X(2). 
S(1) = X(2) LT X(3). 
USE = WEIGHT GE 50. 

T = W GT 25. 

USE = X(2) NE X(3)• 
X(1) = X(2) EQ X(3). 
USE = TEMPA AND T. 
USE = XMIS IF T. 


X (2) = RNDU(8693). 


X(5 ) = RNDG(8896). 


BMDP STATISTICAL PROGRAMS= = = = = = = = e = = = = = = = = = = = = = = = = = = - = = -- : ... 

BMDP1D SIMPLE DATA DESCRIPTION 

Computes summary statistics, lists all data or lists 
selected cases and stores data on a Save File. The 
program allows the user to create new grouping variables 
which are added to the data as additional variables if the 
data are saved on a Save File. 

OUTPUT 


For each variable; mean, 
ooef. of variation, largest 
and smallest standardized 
number of acceptable cases. 


std. dev., std error of mean, 
and smallest values, largest 
(z) -scores, range, and total 


For each grouping variable generated; mean, std. dev. and 
total II of cases grouped. 


75 






BMDP2D - FREQUENCY COUNTS 


CONTROL LANGUAGE USAGE 

PARAGRAPH PARAMETERS COMMENT 


/PROB 

/INPUT 

/VARIAB 

/GROUP 

/TRANSF y 

/SAVE 

/END 

/PRINT 


see MAIN CHART 


MISS. 


(Optional) MAX. 

MIN. 


DATA. 


Print oases containing missing 
values (no/prev.) 

Print oases with values exceeding 
maximum specified (no/prev.) 

Print cases with values less than 
minimum (no/prev.) 

Print all cases (no/prev.) 


BMDP2D FREQUENCY COUNT ROUNTINE 

Counts and lists all distinct values in ascending order 
and computes summary statistics of selected variables. 
Values may be truncated or rounded by different amounts 
for different variables. All computations are performed 
on the truncated or rounded values. 


For each selected variable, it of cases; max, min, range, 
1/2 interquartile range, mean, median, mode, std. dev., 
std. error of the mean, skewness and kurtosis. 

Each distinct value, its frequency, % of cumulative. 

A histogram for each selected variable. 

A line plot of the computed statistics. 


CONTROL LANGUAGE USAGE 


PARAGRAPH PARAMETERS 


COMMENT 


/PROB 

/INPUT 

/VARIAB 

/TRANSF 

/SAVE 

/END 

/COUNT 


j- see MAIN CHART 
TRUNC = (//)/). 


(optional) or 

ROUND=(if)tf . 

TRUNC=(3)10. 


R0UND=(3)1. 
ESTIM. 


truncate or round the specified 
variable(s) to the nearest unit 
given, (none) 

truncates variable 3 to next 
lower 10. 

(13.82 becomes IN) 

Requests a series of location 
statistics, (no) 


BMDP3D t-TEST AND T SQUARED ROUTINE 


Computes one and two sample t-tests along 
associated probability levels. Paired t-test 
creating a transformed difference variable and 
one sample t-test on this difference variable. 


with their 
is done by 
performing 


OUTPUT 


means, std. dev., std. error of means, max. and mins. 


a 20-interval histogram for each group. 
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F value for comparison of variances, and pooled seoarat, 
t-values. y 


two-tailed x values for each t and F. 

, Hotelling's T squarec 


optionally Mahalanobis 
significance. 


CONTROL LANGUAGE USAGE 


/PROB 

/INPUT 

/VARIAB 

/GROUP 

/TRANSF 

/SAVE 

/END 

/TEST 


> see MAIN CHART 


TITLE=’text’. 
VARIAB = list. 


Prints the within group correla¬ 
tion matrices if desired. 

(no/prev.) 

test title if desired.(blank) 
variables to be tested.(all 
variables) 

calculate Hotelling's T squared 
(no/prev.) 

Use only complete cases, (no/ 
prev.) 

specify specific groups to be 
compared if not all are desired, 
(all) 


BMDPND ALPHANUMERIC FREQUENCY COUNT ROUTINE 

ohar ^ters found in single column fields on the 
r™ 6 p r° Eram assumes that a single column 
Ind non , d f °r each variable in A1 format. All IBM 026 
and 029 keypunch characters are considered legal and the 

screening p^gr^ 13 ^ “riStly 

OUTPUT 

Frequency counts of all distinct characters. 

CONTROL LANGUAGE USAGE 


/PROB 
/INPUT 
/VARIAB 
/END 
/PRINT 

(optional) 


See MAIN CHART 


ALPHA = 

' Character' 


NUMERIC = 

'charac ter' 


BLANK r 
'character' 


PERIOD = 

' character' 


SIGN = 

' character' 


SYMB = 
'character' 


..No LABEL, CONTENT, CODE, TYPE 
, .NAME and LABEL only. 

A character to replace all 
alphabetic characters when 
printing data, (none) 

A character to replace all 
numerical characters, (none) 

A character replacing all 
blanks, (none) 

A character replacing all 
periods or decimals, (none) 

A character replacing all plus 
or minus signs, (none) 

A character replacing all 
symbols, (none) 
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BMDP5D - UNIVARIATE PLOTTING 

ILL = A character replacing all 

'character'. illegal values in the data. 

PRINT SYMBOL IS ’$' . 


BMDP5D UNIVARIATE PLOTTING 

Prints univariate plots. For each plot, cases belonging 
to one or more groups can be used. Cases from different 
groups can be identified by distinct letters. Plots can 
be printed for all groups in one plot or for each group 
individually. The size of the plot can be controlled. 

OUTPUT 

- A histogram and/or a cumulative histogram. 

- A normal probability plot, half-normal probability plot, 
and/or cumulative distribution plot. 

- The total frequency count, mean and standard deviation 
(for each group in a plot). 

CONTROL LANGUAGE USAGE 

PARAGRAPH PARAMETERS COMMENTS 


/PROB 

/INPUT y 

/VARIAB 

/GROUP 

/TRANSF 

/SAVE 

/END J 
/PLOT 

(required, 
may be 
repeated) 


see MAIN CHART 


VARIAB = list. 
TYPE = name. 


variables to be 
(ALL/prev.) 

One or more of: 
NORM, HALFNORM, 
/prev.) 


plotted. 

HIST, CHIST, 
DNORM, CUM(HIST 


GROUP = list. Groups used in one plot if not 
all. (ALL) 


SIZE = #1,#2. Length of horizontal & vertical 
axis in characters & lines for 
NORM, DNORM, HALFNORM, AND CUM. 
(80,50/prev .) 

SCALE = in, 112 . Frequency scale of HIST or CHIST 
SCALE = 10,2 start printing 
after a frequency of 10 is 
achieved and print one symbol 
every time 2 or more data points 
are in the cell.(0,1) 

MAX = /(list. . Upper limits for the variables, 
one per variable (none/prev.) 

MIN = It list. Lower limits, (none/prev.) 


BMDP6D BIVARIATE PLOTTING 

Prints bivariate scatter plots. Cases belonging to one or 
more groups can be used for each plot. The plots may have 
frequencies plotted for each point, or distinct letters 
may be used to indicate group membership. Plots can be 
printed for all groups in one plot for each group in 
separate plots. 

OUTPUT 
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- Scatter plots of each specified pair of variables using 
cases of desired groups. 

Mean, std. dev., equation of simple linear regression of 
each variable on the other. 

CONTROL LANGUAGE USAGE 


/PROB 
/INPUT 
/VARIAB 
/GROUP 
/TRANSF 
/SAVE 
/END 
/PLOT 
(required, 
can be repeated) 

YVAR = List. 


see MAIN CHART 


XVAR = List. 


PAIR, or 
CROSS. 


GROUP = List. 


COMMENTS 


Variables for X-axis 
(1st var/prev.) 

Variables for Y-axis 
(1st var./prev.) 

PAIR tabulates n XVAR 

by m YVAR. Cross tabulates 

each XVAR by XVARS (PAIR). 

Groups to use in plots if all 
are not wanted.(all in one plot) 


FREQ. 

CHAR = GROUP. FREQ causes frequency plots. 

VAR. GROUP and VAR causes group or 
variable membership to be 
plotted. 


STAT. 

SIZE = (/list. 

MIN = II list. 
MAX = (/list. 
XCOM. 

YCOM. 

COM. 


Requests statistics to be 
printed. (NO STAT/prev.) 

Length of the X & Y axis in 
characters and lines. 

(70, 42/prev.) 

Lower limits, one per var. 
(none/prev.) 

Upper limits, one per var. 

(none/prev.) 

All XVARS are plotted on 
common X-axis.(NO) 

All YVARS are plotted on 
common Y-axis.(NO) 

All plots are plotted in a 
common frame.(NO) 


SYMB - f C1C2'...Symbols used for group or 

variable identification. (1,2,.. 
if CHAR = FREQ, A, B,... 
otherwise) 


DESCRIPTION 

VARIANCE 


OF STRATA WITH HISTOGRAMS AND ANALYSIS OF 


Groups the data into a specified number ( 10) of groups 
based on the values of a grouping variable or on the order 
of t entry of the data. For each variable the program 
prints histograms for each group, side by side. The number 
of intervals in the histograms may be specified or the 
number intervals can be computed by the program. A one- 
or two-way ANOVA is reported. 
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BMDP8D - HISSING VALUE CORRELATION 

OUTPUT 

- Histogram for each variable and group. 

Mean, std. dev., count, max, min. for each variable and 
group. 

One-way or two-way ANOVA for each variable. 

Optional data listing, winsorized mean & confidence 
intervals by group and correlations. 

CONTROL LANGUAGE USAGE 

PARAGRAPH PARAMETERS COMMENTS 


/PROB 

/INPUT 

/VARIAB 

/GROUP 

/TRANSF 

/SAVE 

/END 

/PRINT 

(optional) 


see MAIN CHART 

CORR. 

DATA. 


ORDER. 


Print correlations of histo¬ 
grams for each group & all 
groups, (no/prev.) 

Print input data matrix. 

(no/prev.) 

Print input data in ascending 
order on the histogram vari¬ 
ables. (no/prev.) 


WINSOR. 


/HIST GROUP = List, 

(required) 

VARIAB = List. 


Print winsorized means, 
(no/prev.) 

Variables stratifying data into 
groups, (none) 

Variables for histograms and 
ANOVA. (all/prev.) 


INCR = Values. Interval width in histograms, 

one for each variable.(max-min/ 
INTER). 

INTER = Values. Number of intervals in histo¬ 
gram, one for each variable. 
(0/prev.) 


BMDP8D 


MISSING VALUE CORRELATION 

Computes correlation four different ways from data which 
contains missing values. 

- ALL VALUE: the mean of each variable is computed from all 
existing values of that variable; then deviations from 
these means are used. 

COVPAIR; each element of the covariance matrix is comput¬ 
ed from the existing pairs of values involved, which are 
then used to compute the correlations. 

CORPAIR; each element of the correlation matrix is comput¬ 
ed from existing pairs of values involved. 

- COMPLETE; cases with any missing values are not used. 


OUTPUT 

- Variance and mean of each variable. 
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Optional frequency table of existing pairs, covariance and 
correlation matrices, sum of weights, matrix of means and 
variances. 


CONTROL LANGUAGE USAGE 


PARAGRAPH PARAMETERS COMMENTS 


/PROB 

/INPUT 

/VARIAB 

/TRANSF 

/SAVE 

/END 

/PRINT 

(optional) 


see MAIN CHART 
CONTENT = 
List. 
COVA. 

CORR. 


FREQ. 


MEAN. 


VAR. 


SUMWTS. 


/CORR TYPE= list, 

(optional) 


SEG. 


DATA, COVA or CORR.(DATA) 

Print covariance matrix. 

(no/prev.) 

Print correlation matrix. 
(YES/prev .) 

Print frequency table. 
(YES/prev.) 

Print matrix of means. 

(NO/prev.) 

Print matrix of variances. 
(NO/prev.) 

Print weight matrix. 

(NO/prev.) 

One or more of ALLVALUE, 

COVPAIR,COR PAIR, or COMPLETE. 
(CORPAIR/prev.) 

Matrices too large are proces¬ 
sed in segments.(no/prev.) 


BMDP9D MULTIDIMENSIONAL DATA DESCRIPTION 


Provides cellwise descriptive statistics for cases when 
the data are classified multidimensionally. 

OUTPUT 


For each variable, cell frequencies, cell means and cell 
std. devs. 

For each variable across all cells, chi-square test of 
equality of cell frequencies, F test of cell means and 
Bartlett's test for homogeneity of cell variances. 

Optional plots of cell mean shifts. 


CONTROL LANGUAGE USAGE 

PARAGRAPH PARAMETERS COMMENTS 


/PROB 

/INPUT 

/VARIAB 

/TRANSF > 

/SAVE 

/GROUP 

/END 

/TABULATE 

(required) 


see MAIN CHART 


GROUP = List. Variables used for indexing 
(required) or grouping.(none/prev .) 

VARIAB = List. Variables for cellwise des¬ 
cription. (all variables). 

MARGIN s'Cl'C2.Literal map indicating 
marginals at which variables 
are desribed. (all groups). 

PLOT. Request to plot means of each 

cell, (no/prev.) 


BMDP1F 


TWO-WAY CONTINGENCY TABLES 








BI1DP2F - CONTINGENCY TABLE ANALYSIS 

Forms two-way contingency tables from raw data or from 
cell frequencies. 

OUTPUT 

Independence of the row and column classifications. 

Fisher's exact probability for 2x2 tables. 

Measures of association and correlation between the 
classifications as two sets of ranks. 

Optional original frequency, % total table frequency, % of 
row total, % of column total, expected value, etc. 

CONTROL LANGUAGE USAGE 

PARAGRAPH PARAMETERS COMMENTS 


/PROB 

/INPUT >• see MAIN CHART 
/VARIAB 

/CATEG . 


/TRANSF 

/SAVE CONTENT = list. 

/END J 

/PRINT OBS. 

(optional) 

EXP. 

ROWP. 

COLP. 

TOLP. 

DIF. 

STAND. 

ADJ. 

SMOOTH. 

FREEMAN. 

/STAT CHISQ. 

(optional) 

CONT. 

LARCHI. 

FISH. 

TETRA. 

CORR. 

SPEAR. 

GAM. 

LAM. 

LSTAR. 

TAU. 

UNC. 

MCN. 


Required if any ROW, COL, or 
COND variable has more than 10 
values. 


Print observed frequency table. 
(YES/prev.) 

Print expected cell values. 
(NO/prev.) 

Print row percentages.(NO/prev.) 
Print column percentages. 

(NO/prev.) 

Print percentages of total 
frequency.(NO/prev.) 

Print differences.(NO/prev.) 

Print standardized deviations. 
(NO/prev.) 

Print adjusted standardized 
deviations.(NO/prev.) 

Print smoothed values.(NO/prev .) 
Print table of Freeman-tukey 
deviates. (NO/prev.) 

Print Chi-square and Yate's 
correlation.(yes/prev.) 

Print C,phi,Cramers V, alfa, etc. 

(no/prev.) 

Print likelihood ratio test. 

(no/prev.) 

Print Fisher’s exact prob. 

(no/prev.) 

Print tetrachoric correlation, 
(no/prev.) 

Print product-moment correla¬ 
tion . (no/prev.) 

Print Spearman's corr. coef. 
(no/prev.) 

Print gamma, Kendall’s taub, 
and Stuart's tauc Somer's D. 

(no/prev.) 

Print lambda,and lambda asym. 
(no/prev.) 

Print Goodman and Kuskal lambda 
asym.(no/prev.) 

Print Goodman and Kuskal tau asym. 
(no/prev.) 

Print measures of uncertainty 
U and Uasym. (no/prev.) 

Print McNemar's test of symmetry 
(no/prev.) 


One or both DATA or OBS. 
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/TABLE 


COL=list. 


Identification of col. variables 
(none) 

ROW=list. Identification of row variables. 

(Required (none) 

may be re- PAIR. OR PAIR tabulates variable pairs in 

peated) CROSS. order. CROSS tabulates each row 

var. vs. each col. var. (PAIR). 

C0ND=li3t. Names or numbers of table condi¬ 

tioning variables, (none) 

COUNT=list. Variables used to sum freq. 

counts in ROW, COL and/or COND 
tables, (none) 

MIN = //. Minimum expected frequency for 

collapsing tables. (0) 

DELTA=fl. Value added to each cell fre¬ 

quency before calculating 
statistics. (0) 


BMDP2F TWO-WAY CONTINGENCY TABLES — EMPTY CELLS AND THE 

IDENTIFICATION OF DEPARTURES FROM INDEPENDENCE. 

Analyzes two-way frequency tables. Differs from BMDP1F in 
that: 

- Cells may be excluded from the table by declaring them 
empty. This allows a test of independence using all cells 
not declared to be empty and is identical to the usual 
Chi-square when there are no empty cells. 

- Cells or patterns of cells that contribute to a departure 
from independence or quasi-independence may be identified. 

This program does not compute measures of association, corre¬ 
lation, or prediction because of empty cells. 

OUTPUT 

- Observed frequency table and associated Chi-square. 
Expected values for each cell. 

- Stepwise chi-square, cell elimination criteria,prevoiusly 
eliminated cells. 

Final summary table. 

Optional tables. 

CONTROL LANGUAGE USAGE 

PARAGRAPH PARAMETER COMMENTS 



see MAIN CHART 


. Required if any ROW, COL or COND 

variable has more than 10 values. 
C0NTENT=name. One only of DATA or OBS.(DATA) 


/TABLE ROW,COL,PAIR or CROSS. 

(Required, see gMDPIF 

may be COND,COUNT,OBS,DELTA 

repeated) FIT. Print table of fitted cell 

values.(no) 

EMPTY (i) = Row and column cell in indices 

list. empty a priori, (i) refers to the 

table number. EMPTY (3) = 1,1,2, 
3.(none) 

C0NV=//. Convergence criterion to find 

fitted value. (.05) . 

ITER = #. Max iterations for fitting quasi¬ 

independent model. (20) 
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BMDP3F - MULTIWAY CONTINGENCY TABLES 


/STEP 

CRIT=name. 

Criterion for cell elimination 

(optional, 

(required) 

CHISQ chi-square,LRCHI likelihood 

may be re¬ 


ratio, USTAN largest absolute 

peated) 


std. quasi-residual, QUADJ 
largest absolute adjusted std. 
quasi-residual, STAN max. abso¬ 
lute adjusted std. residual, ADJ 
max. absolute adjusted std. 
residual, DIF max. absolute 
residual, (none) 


MAXST=tf. 

Maximum // of steps.(degrees of 
freedom) 


ALPHAS. 

Level of significance for testing 
quasi-independence.(.05) 


ELIM=//. 

Criterion limit to eliminate 
cells. (1.0) 


REEN=#. 

Criterion limit to reenter cells. 
(1.0) 


PRINTS. 

Steps at which tables below are 
printed.(none) 


FIT. 

Print table of fitted values.(no) 


CHISQ. 

Print resulting chi-square when 
cell is eliminated.(no) 


LRCHI. 

Print table of likelihood ratio 
chi-square.(no) 


QUSTAN. 

Print table of quasi-std. devi¬ 
ations . (no) 


QUADJ. 

Print table of quasi-adjusted 
std. devs.(no) 


STAN. 

Print table of standard residuals 
(no) 


ADJ. 

Print table of adjusted std. 
residuals.(no) 


DIF. 

Print table of residuals.(no) 


BMDP3F MULTIWAY CONTINGENCY TABLES 

Forms and analyzes multi-way contingency tables. The 
purpose of the analysis is to obtain a description of the 
relationships between factors of the table, either by 
forming a model for the data or by testing and ordering 
the inportanoe of the interactions between the factors. 
The analysis is based on fitting a log-linear model to the 
cell frequencies. 

Input may be the usual variables-by-case format, cell 
indices and frequency counts, or a multiway frequency 
table. 


- The original frequency table and marginal totals. 

- A test for each k that the k-factor interactions are zero. 

An ANOVA-like table of likelihood ratio Chi-square stats. 

Estimated arameters and tests of fit for a user defined 
model. 


COMMENTS 


CONTROL LANGUAGE USAGE 
PARAGRAPH PARAMETER 


/PROB 


/INPUT 

^ see MAIN CHART 

/VARIAB 


/CATEG 


/TRANSF 


/SAVE 

C0NTENT=name 

/END J 


/INPUT 

TABLE = mst. 


Required if any ROW,COL or COND 
var. has more than 10 different 
values. 

One only of DATA or OBS.(DATA) 
Number of categories for each 








index if a frequency table is 
input fastest index first.(none) 

/TABLE 

INDEX=List. 

Variables used as indices for 

(Required, 


table.(none) 

may be 

SYMB=List. 

Symbols assigned to INDEX vari¬ 

repea ted) 


ables, in output.(none) 


COUNT=List. 

Variables that are frequency 
counts.(none) 


COND=List. 

For every level of these vari¬ 
ables, a table is formed.(none) 


OBS. 

Print observed frequency table. 

(yes/prev.) 


MARG=#. 

Print marginal subtables up to 
order k. (O/prev .) 


DELTA=Value 

Add Delta to each frequency be¬ 
fore analysis.(O/prev.) 


INTERS. 

Max II of iterations.(20/prev.) 


CONV = (/. 

Criterion for convergence of 
model fitting.(.1/prev .) 


SIMUL. 

Calculate tests-of-fit of models 
of full order.(yes/prev .) 


ASSOC. 

Calculate tests of partial and 
marginal association.(yes/prev.) 

/FIT 

ALL. 

Calculate all models if 2 or 3- 
way table.(no/prev.) 

MODEL=List. 

Hierarchical model to be fitted. 

(Optional, 
may be 


(none) 

repeated) 

EXP. 

Print table of fitted (expected) 
values.(no) 


STAN. 

Print std. residuals.(no) 


FREEMAN. 

Print Freeman-Tukey deviates.(no) 


LAMB. 

Est. log-linear parameters of 


BETA. 

Est. multiplicative parameters of 
MODEL.(no) 


ORTH=List. 

Fit orthogonal polynomials for 



SYSMBOLS spec ified.(none) 


PROB=». 

Print LAMB, BETA, ORTH only if 
test-of-fit of model is greater 
than this value . (. 05) 


ADD. 

Form models that contain MODEL as 
a subser.(no) 


DELETE. 

Form models contained by MODEL, 
(no) 


BMDP1M CLUSTER ANALYSIS ON VARIABLES 

Clusters similar variables by using an initial measure of 
association between pairs of variables to form a cluster 
of the two most similar variables and then using an amal¬ 
gamation rule to form further clusters. 

OUTPUT 

- Summary table for the clustering process. 

- Tree diagram of clusters. 

Table of similarity or distance matrix scaling. 

Explanation of tree diagram. 

CONTROL LANGUAGE USAGE 

PARAGRAPH PARAMETER COMMENTS 


/PROB 

/INPUT 

/VARIAB 

/TRANSF 

/SAVE 


TYPE=name One only of DATA, COVA, CORR, ANG, 

i SIMI, DIST. (DATA/prev.) 

f see MAIN CHART 
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BMDP2H - CLUSTER ANALYSIS ON CASES 

/PRINT SHADE. 

(Optional) 

CORR. 

/PROC VARIAB= 

(Optional) List. 

MEAS=name. 


LINK=name. 


Print shaded correlation matrix, 
(no/prev.) 

Print correlation matrix. 

(no/prev.) 

Variables to be clustered.(all in 
use/prev.) 

Association measure: CORR 
Correlation, ABSCORR absolute 
correlation, ANG arccos (Corre¬ 
lation), ABSANG arccos (abs corr) 
(corr) 

Amalgamation rule, one only: 
SINGLE, COMP, AVE.(SINGLE/prev.) 


BMDP2M CLUSTER ANALYSIS ON CASES 

Clusters cases on the basis of one of four available 
distance measures. 

OUTPUT 

A diagram indicating clustering of cases. 

- Optional standardized input data, initial distance case 
matrix and vertical clustering diagram. 

CONTROL LANGUAGE USAGE 

PARAGRAPH PARAMETER COMMENTS 


/PROB 
/INPUT 
/VARIAB 
/TRANSF V 
/SAVE 
/END 
/PROC 
(Optional) 


see MAIN CHART 


SUMOFSQ 


SUMOFPsp. 

CHISQ. 

PHISQ. 

STAND. 

WEIGHT=name. 


/PRINT DATA. 

(Optional) DIST. 


REG. or VERT. 


AMALG. 


SHADE. 


Square root of sums of squares 
of differences of distances be¬ 
tween clusters.(yes/prev.) 
pth root of sum of pth powers of 
diff. as distance measure.(2) 
Chi-square as distance.(no) 
Cbi-square/n as distance.(no) 
Standardize the input data.(yes) 
Variable containing case weights 
(unit weights), (no weights/prev.) 
Print std. data, (no/prev.) 

Print initial std. distances. 

(no/prev.) 

Print horizontal or veritcal 
triangle tree diagram.(VERT/ 
prev.) 

Print table of distances, case 
weights and averages by step 
(yes/prev.) 

Print shaded form of distance 
matrix.(no/prev.) 


BMDP3M BLOCK CLUSTERING 

Simultaniously clusters cases and variables of a data 
matrix. Three types of clusters are thus discovered: 
clusters of data values (the blocks), clusters of cases, 
and clusters of variables. 

CONTROL LANGUAGE USAGE 

PARAGRAPH PARAMETER COMMENTS 


/PROB ") 
/INPUT I 
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/VARIAB INTERs#list. 

/TRANSF 

/SAVE 

/END /see MAIN CHART 
/PASS TOTALS. 

(Optional) 

THRESH=#list. 


Number of intervals, one per 
variable (first 35 values) 


Number of passes on rows and 
columns. (10) 

Number indicating threshold at 
each pass, (increase from 2/ 
TOTAL to 1) 


BMDPRM FACTOR ANALYSIS 

Performs a factor analysis of either a correlation matrix 
or a covariance matrix. The input may be raw data, 
correlation matrix, covariance matrix, factor loadings or 
factor score coefficients. 


OUTPUT 


- Means, std . dev., coefficients of variation, maximums and 
minimums, standard scores. 

- Correlation matrix, covariance matrix. 

- Factor loadings and their plots. 

- Squared multiple correlations and Kaiser's measure of 
sampling adequacy. 

- Factor correlations and factor structure. 

- Factor scores, factor score coefficients, factor score 
covariances and factor score plots. 

CONTROL LANGUAGE USAGE 


PARAGRAPH PARAMETER. COMMENTS 


/PROB 

/INPUT 

/VARIAB SCORE = list. 

/TRANSF WE IGHTsvar. 

/SAVE 

/END ^/see MAIN CHART 
/INPUT TYPEsname. 

(Required on 

first pro- CONTENTsname. 
blem) 

FACT0R=#. 


/SAVE 
(Optional) 
/PRINT 
(Optional) 


SHAPE=name. 

CONTENTmame. 

CASE=#. 

FSCOREs#. 


HILEVs#. 

LOLEVs#. 


STAN. 

COVA. 


CORR. 


INV. 

PART. 


FSTR. 


Variables used for factor scores, 
(none) 

Variable containing case weight, 
(no weight/prev.) 

One of DATA, CORR, COVA., LOAD, 
FSCF.(DATA or CONTENTs) 

One of DATA, CORR, COVA, LOAD, 
FSCF. (DATA) 

Number of factors when TYPEsLOAD 
or FSCF. (NONE) 

One of SQUARE, LOWER. (SQUARE) 

One or more of DATA, COVA, UFLD, 
RFLD, FCOR. (DATA) 

Number of cases for which data is 
printed. (5/prev.) 

Number of factors for printing 
factor scores, (all) 

High level factor loading. 
(.5/prev.) 

Low level factor loading. 
(.25/prev.) 

Print standard scores.(no/prev .) 
Print covariance matrix. 

(no/prev.) 

Print correlation matrix. 

(yes/prev.) 

Print inverses.(no/prev.) 

Print partial correlations. 

{no/prev.) 

Print factor structure matrix 
(oblique rotation).(no/prev.) 
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BMDP6M - CANONICAL CORRELATION 


/PLOT 

(Optional) 


/FACTOR 
(Optional) 


/ROTATE 

(Optional) 


SHADE. 

RESI. 

INITIAL:#. 

FINAL://. 
FSCORE = //. 

CONST=value. 

NUMBER://. 


METHOD=name. 


ITERATE://. 

EPS=//. 

COMMUN=name. 


FORM:name. 


METHOD:name. 


GAMMA:value. 


CONST=value. 

MAXIT://. 

NORMAL. 


Print shaded correlation. 

Print residual correlations. 

(no/prev.) 

Number of unrotated factors for 
factor loadings plots.(A) 

Number of rotated factors (A) 
Number of factors for factor 
score plots. (3) 

Value limiting factors to 
eigenvalues >CONST.(1/prev.) 
MAXIMUM // OF FACTORS OBTAINED. 

(// of variables) 

One only of PCA principle compo¬ 
nents, MLFA maximum likelihood, 
LJIFFY Kaiser's second generation 
Little Jiffy, PFA principle factor 
analysis.(PCA). 

Maximum // of iterations . (25/prev .) 
Convergence criteria.(.001/prev.) 
One of the following communality 
options: UNALT unaltered, SMCS 
squared multiple correlations, 
MAXROW, maximum absolute row 
values of the correlation matrix. 
(UNALT for PCA, SMCS otherwise/ 
prev .) 

One of the following for factor 
form: CORR factor correlation 

matrix, COVA covariance matrix, 
0C0RR correlation matrix about 
the origin, 0C0VA covariance 
matrix about the origin . (CORR/ 
prev.) 

NONE no rotation, VMAX varimax, 
DQUART direct quartimin, ORMAX 
quartimax, EQMAX equamax,0RTH0G 
orthogonal with gamma below, 

DOBLI direct oblimin with gamma, 
ORTHOB orthoblique. (VMAX except 
with LJIFFY/prev.) 

Parameter for METH0D:0RTH0G or 
DOBLI. (1 for 0RTH0G, 0.0 for 
DOBLI) 

Convergence criteria for rotation 
( .00001/prev.) 

Max. // of iterations for ratation. 
(50/prev.) 

Request Kaiser's normalization, 
(yes/prev.) 


BMDP6M CANONICAL CORRELATION ANALYSIS 

Computes canonical correlation analysis for two sets of 
variables. Input can be data, a covariance or a correl¬ 
ation matrix. Computations can be done in single or 
double precision. 

OUTPUT 


- Means, std. devs., coefficients of variation, skewness 
kurtosis, mins and maxs. 

- Correlations, covariances and canonical correlations. 

- Eigenvalues associated with each pair of canonical 
variables. 


Coefficients of canonical variables and canonical vari¬ 
ables. 


— Correlation of variables with canonical variable and plots 
CONTROL LANGUAGE USAGE 
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PARAGRAPH 


PARAMETER 


COMMENTS 


/PROD 

/INPUT 

/VARIAB 

/TRANSF L 
/SAVE 
/END 
/PRINT 
(Optional) 


/PLOT 

(Optional) 


/CANON 
(Required) 


SHAPE=name. SQUARE, LOWER. (SQUARE) 

TYPE=name. One of: DATA,COVA,CORR.(DATA) 

WEIGHT=V. Variable containing case weight, 

(none/prev.) 

see MAIN CHART 

CONTENT=list. One or both DATA,COVA. (DATA) 


CASE = 11. 
MATR=list. 


XVAR=names. 

YVAR=names. 
SIZE = /)1 ,//2. 

FIRSTmames. 

SECOND=names . 

TOL = 11. 

PREC=name. 

ZERO. 

CONST = />. 

NUMBERS#. 


Number of cases to print data. 
(5/prev.) 

CORR correlations, COVA covari- 
iances, CANV canonical variables, 
COEF canonical var. regression 
coefficients, LOAD canonical vari¬ 
able loadings. (CORR,LOAD) 

Names or subscripts of x-axis 
vars. (none) 

Y-axis variables, (none) 

#1 is number of columns. II 2 is 
number of rows.(50,50/prev.) 

First set of variables.(none) 

(required) 

Second set of variables, (none) 
(required) 

Tolerance of matrix inversion. 

(.01/prev.) 

SINGLE single precision, DOUBLE 
double precision.(DOUBLE/prev .) 
Assume means zero. (NO ZERO/prev.) 
Limit that canonical correlation 
must be greater than.(0.0/prev .) 
Max. II of canonical variables. 

Ul in smaller of two variable sets 
FIRST and SECOND) 


BMDP7M STEPWISE DISCRIMINANT ANALYSIS 

Performs stepwise multiple group discriminant analysis. 
OUTPUT 


Group means and std. devs. 

F statistics, Wilk lambda and Magalanobis' D square. 

- Classification functions, classification matrices and % 
correct classification. 

- Coefficients for canonical discriminant functions, 
canonical correlations and canonical variables. 

- Plots of the first two canonical variables. 

CONTROL LANGUAGE USAGE 


PARAGRAPH PARAMETER 


/PROB V sei 
/INPUT 
/VARIAB 
/GROUP 

/TRANSF 


/SAVE 

/END J 
/PRINT 
(Optional) 


MAIN CHART 
GROUP=#. 

USErlist. 

PRIOR=#list. 

NCANs#. 

STEP. 

CLASS=#list. 


COMMENTS 


Number of groups if >10 and no 
GROUP paragraph.(none/prev.) 

Groups for within-group cross- 
product matrix.(all groups/prev.) 
Prior probabilities used to com¬ 
pute posterior probabilities. 

(1.0/prev.) 

Number of canonical variables 
saved.(# of groups - 1/prev.) 

Print detailed stats of each step, 
(yes/prev.) 

Print classification functions 
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BMDPAM - MISSING DATA 


POST. 

POINT. 

WITH. 

/PLOT CANON. 

GROUP=list. 

CONTR. 

/DISC TOL = #. 

(Optional) 

ENTER = #1, #2. 
REMOVE=#1 ,112. 
HETHOD=1 or 2. 

LEVEL:#list. 


FORCE=#. 

STEP:#. 
CONTR:valuea. 

JACK. 


when number of variables entered 
equals # on this list. (none/prev .) 
Posterior probabilities and D 
squared are printed.(yes/prev.) 
Values of canonical variables are 
printed.(yes/prev.) 

Print pooled with-group covariance 
matrix, (no/prev.) 

Plot of first two canonical vari¬ 
ables. (yes/prev.) 

Names or subcripts of groups 
used for canonical var. plots. 

(all groups/prev.) 

Canonical variables are based on 
specified contrasts.(no/prev.) 
Number (.001,1) used to guard 
against singularity.(.01/prev.) 
Limiting value of F to enter. 

(A . A/prev.) 

Limiting value of F to remove. 
(3.99,3.99/prev.) 

Method of deletion; 1 use F to 
remove, 2 delete to obtain better 
prediction.(1/prev.) 

List of levels for forcing vari¬ 
ables into equation. Zero level 
not entered. Lower levels entered 
before higher levels. Levels list 
to correspond to variable indices, 
(next higher than any in list) 
Levels up to which variables are 
to be forced into the equation. 
(0/prev.) 

Max. # of steps allowed.(twice 
# of variables/prev.) 

Values specifying a contrast used 
to direct stepping. (contrast to 
equality of group means/prev.) 
Compute posterior probabilities, 
(no/prev.) 


BMDPAM DESCRIPTION AND ESTIMATION OF MISSING DATA 

Describes the pattern of missing values for multivariate 
data and provides estimates of the missing data. Esti¬ 
mates can be means; regression values on the available 
variable with which it is most highly correlated; re¬ 
gression on a set of variables highly correlated or re¬ 
gression on all available variables. 

OUTPUT 

Pattern of missing data. 

- Squared multiple correlations of each variable with all 
others. 

Missing data estimates. 

Bivariate plots. 

Correlation matrix, eigenvalues, covariance matrix. 

CONTROL LANGUAGE USAGE 

PARAGRAPH PARAMETER COMMENTS 


/PROB 

/INPUT 

/VARIAB 

/GROUP ^ see MAIN CHART 
/TRANSF 

/SAVE J C0NTENT:list. One or more of DATA,COVA CORRDICH. 
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/VARIAB 

WEIGHT. 

(DATA/prev.) 

Variable containing case weights. 

(Optional) 


(none) 


GROUPrv. 

Variable used for grouping.(none/ 

/PRINT 

CASE=#. 

prev.) 

<1 of cases to print data. (5/prev .; 

(Optional) 

MATR=List. 

Matrices to print: CORR correlat¬ 



ing, EIGEN eigenvalues, PAT pat¬ 
tern of missing data, FREQ sample 
size 4 % missing for each pair of 

/PLOT 

XVAR=names. 

vars. DIS Mahalanobis' distance 
from each case to center. COVA 
covariances, SUMWT pairwise sums 
of weights, INVCV inverse of cor¬ 
relation matrix, printed as a 
correlation matrix, CORRDICH cor¬ 
relations computed using the dich¬ 
otomies: 1.0 if the value is ac¬ 

ceptable and 0.0 otherwise.(CORR, 
EIGEN, PAT, EST, FREQ, DIS, CORR- 
DICH/prev.) 

X-axis variables.(none) 

(Optional) 

YVAR=names. 

Y-axis variables.(none) 

/EST 

SIZE=#,#. 

Length of X 4 Y axis.(50,50/prev.! 

METHODsname. 

MEAN means, SINGLE regression on 


(Optional) correlated variable, TWO-STEP or 


STEP stepwise regression estimate 
REGRES regression on all variables 
NONE no estimates.(none/prev .) 


ENTERS. 

F to enter . (A/prev .) 

TYPE=name. 

Algorithm for computing univariate 
statistics; COMPLETE complete 
cases only, ALLVALUE all existing 
values(ALLVALUE/prev.) 

VLIM=0. 

Maximum number of missing vari¬ 
ables for each case.(variables 
used -1/prev.) 

CLIM=%. 

Maximum % of cases missing for 
any variable.(50/prev .) 

T0L=tf. 

Tolerance (0,1) for matrix in¬ 

version . ( .0001/prev.) 

ZERO. 

Assume means are zero.(no/prev.) 

MIN. 

Treat values below minimum limit 
as missing.(no/prev.) 

MAX. 

Treat values above maximum limit 
as missingtno/prev.) 

GROUP. 

Estimate missing grouping vari¬ 
able (no/prev.) 


BMDP1R MULTIPLE LINEAR REGRESSION 

Computes a multiple linear regression equation on all data 
and on groups or subsets of data. Equations with or 
without an intercept can be chosen. 

OUTPUT 

Means, std. dev., mins and max for variables. 

Multiple R and std. error of estimate. 

ANOVA table on regression. 

- Optional, covariance and/or correlation matrix. Scatter 
plots, normal and detrended normal probability plots of 
residuals. Partial residual plots. Residual, predicted 
value and data for each case. 

CONTROL LANGUAGE USAGE 

PARAGRAPH PARAMETER COMMENTS 
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BMDP2R - STEPWISE MULTIPLE LINEAR REGRESSION 


/PROB 
/INPUT 
/VARIflB 
/GROUP > 
/SAVE 
/END 
/PRINT 
(Optional) 


/PLOT 

(Optional) 


/REGRES 
(Required, 
may be 
repeated) 


TYPE=name. 

MAIN CHART 
CONTENT=list. 

COVA. 

CORR. 

DATA. 

ANOVA. 

STEP. 

PARTIAL. 

COEF. 

FRATIO. 

SUMMARY. 

RREG. 

VARIAB=list. 

NORMAL.or 
DNORMAL. 
PREP=list. 

SIZE=#1,#2. 

RESID. 

TITLE='text'. 
DEPEND=name. 

INDEPslist. 

TOL=#. 

ENTER = # 1 , //2. 
REM0VEs#1,#2. 

METHODsname. 


TYPEsname. 

LEVELs#list. 


FORCEs#. 

STEPs#. 


One of: DATA,COVA,CORR. (DATA) 


One or both of DATA,COVA. (DATA) 

Print covariance matrix.(no/prev.) 
Print correlation matrix.(no/ 
prev.) 

Print input data, residuals and 
predicted values, (no/prev.) 

Print ANOVA for each var.(yes/ 
prev.) 

Print full output at each step, 
(yes/prev.) 

Print partial correlations, 
(yes/prev.) 

Print regression coefficients. 

(yes/prev.) 

Print F ratios, (no/prev.) 

Print multiple R, R squared, F-to- 
enter and F-to-remove.(yes/prev .) 
Print corr. between reg. coef. 
(no/prev.) 

Variables for plotting. 

(none/prev.) 

Print normal or detrended nor¬ 
mal plots.(no/prev.) 

Variables for partial residual 
plots.(none/prev.) 

Length of X & Y axis. (50,50/ 
prev.) 

Plot residuals and squared re¬ 
siduals against predicted. 

(no/prev.) 

Regression title <80 cols.(blank) 
Dependent variable; required 
first.(none/prev.) 

Independent variables.(all but 
dependent/prev.) 

Tolerance to guard against singu¬ 
lar! ty .(. 01 /prev . ) 

Limits to F-to-enter.(N,N/prev.) 
Limits to F-to-remove.(3.9,3.9/ 
prev.) 

Stepping algorithm F: based on 
same as F with variable inter¬ 
change, R based on increase in 
multiple R, R5WAP same as R with 
variable interchange.(F/prev.) 

ZERO zero intercept, NONZERO, 

FLOAT treat intercept as a vari¬ 
able to be entered only if sig¬ 
nificant. (NONZERO/prev.) 

Number list, corresponding to 
variable order, specifying level 
for forcing entry. Zeros will not 
be entered. Lowest level entered 
first. (highest + 1) 

Level through which variables will 
be forced into equation, (none/ 
prev.) 

Maximum steps, (twice # of vari- 
ables/prev.) 


BMDP2R STEPWISE MULTIPLE LINEAR REGRESSION 

Estimates the parameters of multiple linear regression in 
a stepwise manner. Variables can be forced into the 
equation and equations without intercepts can be used. 

OUTPUT 

- Covariance and/or correlation matrix. 
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Coefficients, std. coef., std. errors, F ratios, partia 
correlations. 


- Summary tables for F-ratios, partial correlations an 
coefficients. 

- Predicted values, residuals and data for each case. 

Scatter plots of predicted and observed values versu 
selected variables. 

CONTROL LANGUAGE USAGE 


PARAGRAPH PARAMETER COMMENTS 


/PROB 
/INPUT 
/VARIAB 

/TRANSF 
/SAVE 
/END 
/PRINT 
(Optional) 


TYPE=name. 
WEIGHT=name. 

see MAIN CHART 
CONTENTrlist. 

COVA. 

CORR. 

DATA. 

ANOVA. 

STEP. 

PARTIAL. 

COEF. 


FRATIO. 

SUMMARY. 


RREG. 

/PLOT VARIAB=list. 

(Optional) 

NORMAL.or 
DNORMAL 
PREP=list. 

SIZE=#1,#2. 

RESID. 


/REGRES 
(Required, 
may be 
repeated) 


TITLE='text' 
DEPEND=name. 

INDEP=list. 

TOL=#. 

ENTER=#1,#2. 
REMOVES#!,#2 


METHODsname. 


TYPEsname. 


LEVEL=#list. 


Any one of DATA,COVA,CORR(DATA) 
Variable containing case weights, 
(none/prev.) 

One or both of DATA.COVA(DATA). 

Print covariance matrix.(no/prev. 
Print correlation matrix.(no/prev 
Print input data, residuals and 
predicted values.(no/prev.) 

Print ANOVA for each variable. 

(yes/prev.) 

Print full output at each step, 
(yes/prev.) 

Print partial correlations, 
(yes/prev.) 

Print regression coefficients. 

(yes/prev.) 

Print F ratios.(no/prev.) 

Print multiple R, R squared, F-to- 
enter and F-to-remove.(yes/prev.) 
Print corr. between reg. coef. 
(no/prev .) 

Variables for plotting. 

(none/prev.) 

Print normal or detrended normal 
plots.(no/prev.) 

Variables for partial residual 
plots.(none/prev.) 

Length of X&Y axis.(50,50/prev.) 
Plot residuals and squared re¬ 
siduals against predicted. 

(no/prev.) 

Regression title <80 cols.(blank) 
Dependent variablejrequired 
first.(none/prev.) 

Independent variables.(all but 
dependent/prev.) 

Tolerance to guard against singu¬ 
larity. ( ,01/prev.) 

Limits to F-to-enter.(N,N/prev.) 
Limits to F-to-remove.(3.9,3.9/ 
prev.) 

Stepping algorithm F: based on 
F-to-remove & F-to-enter,FSWAP 
same as F with variable inter¬ 
change, R based on increase in 
multiple R,RSWAP same as R with 
variable interchange.(F/prev.) 

ZERO zero intercept, NONZERO, 

FLOAT treat intercept as a vari¬ 
able to be entered only if sig¬ 
nificant. (NONZERO/prev.) 

Number list, corresponding to 
variable order, specifying level 
for forcing entry. Zeros will not 
be entered. Lowest level entered 
first. (highest + 1) 
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BMDP3R - NONLINEAR REGRESSION 

FORCE=#. Level through which variables 

will be forced into equation. 

{none/prev.) 

STEP=#. Maximum steps, (twice II of vari¬ 

ables/ prev.) 


BMDP3R NONLINEAR REGRESSION 

Least squares fit to a non linear function. Five most 
frequently used non linear functions are built in, others 
may be specified with FORTRAN statements. 

OUTPUT 

Parameter values and residual sum of squares. 

Asymptotic correlations and std. dev. for est. parameters. 
Predicted and observed value for the dependent variable. 
Residual value. 

Observed values for independent variables. 

SPECIFYING SPECIAL REGRESSION FUNCTIONS IN FORTRAN 

Any function can be fitted to the data by specifying the 
function and its derivatives in the FORTRAN subroutine 
FUN. The FORTRAN subroutine containing the correct 
regression function and derivatives must then become part 
of the program BMDP3R for that problem run If the 
derivatives of the desired nonlinear expression are 
difficult to obtain or specify, the user is advised to use 
BMDPAR which does not require the derivatives to be 
specified. The following is an example subroutine FUN for 
the regression equation 

f=p(1)*EXP(P(2)*t)+p(3)*EXP(P(4)*t): 

SUBROUTINE FUN(F,DF,P,X,N,KASE,NVAR,NPAR,IPAS,XLOSS) 
DIMENSION DF(NPAR),P(NPAR),X(NVAR) 

IMPLICIT INTEGER*4 (I-N) 

DF C1)=EXP(P(2)*X(3)) 

DF(2)=P(1)*X(3)*DF(1) 

DF(3)=EXP(P(4)*X(3)) 

F = P(1)*DF(1)+P(3)*DF(3) 

RETURN 

END 

CONTROL LANGUAGE USAGE 

PARAGRAPH PARAMETER COMMENTS 


/PROB 
/INPUT 
/VARIAB 
/TRANSF L 
/SAVE 
/END 
/PLOT 
(Optional) 


see MAIN CHART 


VARIAB=list. 
NORMAL.or 
DNORMAL. 


SIZE=#1 ,11 2, 
RESID. 


/REGRES TITLE='text'. 

(Required) NUMBERS. 

DEPEND=name. 


TOL=#. 

WEIGHTmame. 


Variables for plots.(none/prev.) 
Print normal or detrended normal 
probability plots of residuals. 

(no/prev.) 

Length of X A Y axis.(50,50/prev .) 
Plot residuals and squared resid¬ 
uals against predicted.(no/prev.) 
Regression title <160 cols.(blank) 
Select built-in function.(See 
manual) 

Dependent variable required, 
(none/prev.) 

Singularity tolerance.(.001/prev .) 
Variable containing the case 
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PARAM;//. 

CONST = //. 

ITER = </. 
CONVER;//. 

HALV=//. 

MEANSQ;//. 

PASS;//. 

LOSS. 

/PARAM INIT;//1 ist. 

(Optional) 

NAME;names. 

MAX;// list. 
MIN;// list. 


weights.(none/prev.) 

// of parameters in the regression 
function required.(none/prev.) 

The number of constraints on para¬ 
meter s .( none/prev . ) 

Max // of iterations. (50/prev .) 
Convergence tolerance.(. 00001 / 
prev.) 

Max It of increment halving.(5/ 
prev.) 

// used to compute asymptotic 
standard deviations.(none/prev.) 

It of passes of data at each itera¬ 
tion . (1/prev.) 

Loss function specified in FORTRAN 
subroutine.(no/prev.) 

Starting values for regression 
parameters.(0,0,.../prev.) 

Names of regression parameters. 

Up to 8 characters.(P(subscript)/ 
prev.) 

Maximum values for regression 
parameters.(none/prev.) 

Lower limits for regression 
parameters.(none/prev.) 


BMDPAH REGRESSION ON PRINCIPLE COMPONENTS 

Computes a regression analysis for each dependent vari¬ 
able on a set of principle components computed from the 
independent variables. 

OUTPUT 


Means, std. devs. and correlation or covariance matrix. 

- Eigenvalues and eigenvectors (principle component co- 
elricients). 

inHlnfnH 0 "*. 00 ? f [, ioients for principle components and 
independent variables. 

- Residual sums of squares and F ratios at each step. 

- Optional principle component scores, scatter plots of raw 
data, and/or normal and detrended normal prob. plots of 
residuals. 


CONTROL LANGUAGE USAGE 
PARAGRAPH PARAMETER 


/PROB 
/INPUT 
/VARIAB 
/TRANSF l 
/SAVE 
/END 
/PRINT 
(Optional) 


see MAIN CHART 
SCORE. 


/PLOT VARIAB;list. 

(Optional) NORMAL.or 

DNORMAL. 


/REGRES 
(Required, 
may be 
repeated) 


SIZE=#1,02. 
RESID. 

TITLE;'text' 
DEPEND;name. 

INDEP=list. 


STAND. 


COMMENTS 


Print standardized principle 
component scores for each case, 
(no/prev.) 

Variables included.(none/prev.) 
Print normal or detrended normal 
probablility plots of residuals, 
(no/prev.) 

Length of X A Y axis.(50,50/prev.) 
Plot residuals.(no/prev.) 
Regression title <80 cols.(blank) 
Dependent variable required 1st. 

(none/prev.) 

Independent variable list.(those 
in USE statement/prev.) 

Regression on standardized vari- 
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BMDP5R - POLYNOMIAL REGRESSION 


EIGEN. 

CORR. 

LIMIT=#1,#2. 


ables. (no/prev.) 

Enter components by size of 
eigenvalues, (no/prev.) 

Enter components by size of abso¬ 
lute value of correlation with 
dep.var. (yes/prev.) 

Smallest correlation and eigen¬ 
value to cause variable to enter. 
(.01,. 01/prev.) 


BMDP5R POLYNOMIAL REGRESSION 

Fits a polynomial in one variable to the dependent 
variable. The form of the regression equation is: 

y(i)=B(0)+B(1)«x+B(2)*x**2+...+B(r)*x(i)**r+e(i) 


OUTPUT 


Correlation matrix for regression coefficient estimates. 
Estimates of coefficients and their std. errors. 


A summary table of goodness-of-fit statistics. Optionally 
residual and fitted polynomial values, scatter plots of 
observed and predicted or residual values vs. independent 
variable, normal and detrended norm prob, plots of resi¬ 
duals, coefficients in the expansion of each polynomial P 
(x) in powers of x. 

CONTROL LANGUAGE USAGE 

PARAGRAPH PARAMETER COMMENTS 


/PROS 

/INPUT 

/VARIAB WEIGHT=V. 

/TRANSF see MAIN CHART 

/SAVE 

/END 

/PLOT RREG. 

(Optional) 

DEGREE=#list. 


/PLOT 

(Optional) 


/REGRES 

(Required) 


ORTHPOL. 

NORMAL.or 
DNORMAL. 
SIZE=#1,#2. 
DEGREE=#list. 

TITLE='text'. 
DEPEND=name. 
INDEP=list. 
DEGREEr#. 


Variable for case weights.(none/ 
prev.) 


Print correlation matrix of reg. 
coefficients, (no/prev.) 

Print residuals at stated degrees 
of regression.(none/prev.) 

Print orthogonal poly, coeffi¬ 
cients. (no/prev.) 

Print normal plots.(no/prev.) 

Length of X & Y axis.(50,50/prev.) 
Print plots at specified degrees 
of regression.(none/prev.) 
Regression title <80 cols.(blank) 
Dependent variable.(none/prev.) 
Independent variables.(none/prev.) 
Degree of final regression equa¬ 
tion. (3/prev.) 


BMDP6R PARTIAL CORRELATION AND MULTIVARIATE REGRESSION 

Computes the partial correlations of a set of variables 
removing the linear effects of a second set of variables. 
The program can also be used for regression, especially if 
multiple dependent variables are present. Computation can 
be performed in either single or double precision. 

OUTPUT 

Means, std. devs., coeff. of variation, skewness,kurtosis, 
min.,max. 


96 







- Correlations, covariances, partial correlations and par- 
tial covariances. 

- Regression coef. , std. coef., std . errors for coef. 

- Covariances and correlations for reg. coefs. 
t-tests and significance levels for reg. coefs. 

- Squared multiple correlations and significants levels. 

- Residuals from partial correlation analysis. 

Bivariate plots of variables and residuals. 

CONTROL LANGUAGE USAGE 

PARAGRAPH PARAMETER COMMENTS 


One of: DATA, COVA, CORR. (DATA) 
SQUARE or LOWER for CORR or COVA. 
(SQUARE/prev.) 

One or more of DATA,COVA,COVAPART. 
(DATA) 

Number of cases of data printed. 
(5/prev.) 

One or more of CORR, COVA, PART, 
COVAPART, RESI, COEFF, TTEST, 

STANC std. reg. coefs., RREG cor¬ 
relations of the reg.coefs.(CORR. 
PART) 

X-axis variables, (none) 

Y-axis variables, (none) 

Size of plots. (50,50/prev.) 

Print normal prob. plots.(no/prev.! 
Names of dependent variables re¬ 
quired. (none) 

Independent variables required. 

(none) 

Tolerance for matrix inversion. 

(.01/prev.) 

DOUBLE double precision, SINGLE 
single precision. (DOUBLE/prev.) 
Assume zero means.(no/prev.) 


/PROB \ 
/INPUT 
/VARIAB 
/TRANSF 
/SAVE 
/END 
/PRINT 
(Optional) 


see MAIN CHART 
TYPE=name. 
SHAPE=name. 

CONTENT=list. 


MATR=names. 


/PLOT 
(Optional) 


/REGRES 

(Required) 


XVARslist. 
YVAR=list. 
SIZE=#1 , 112. 
NORMAL. 
DEPEND=list. 

INDEPrlist. 

T0L=/(. 

PREC=name. 

ZERO. 


BMDP9R ALL POSSIBLE SUBSETS REGRESSION 


Identifies 'best' subsets of predictor variables. Best is 
defined in terms of the same R-squared, adjusted R- 
squared or Mallow’s Cp. 


OUTPUT 


Means, std. devs., coef. of variation, minimums, maxi- 
mums, minimum std. scores, maximum std. scores, skewness 
and kurtosis. 

Correlation and covariance matrix. 

Bivariate plots. 

For best M subsets, selection criteria, regression coef. 
and F-ratio. 

For best subset, R-squared, F ratio and significance, reg. 

=nd f t'i Std ‘ std - errors, t-tests, signifi- cance 

and tolerance, etc. 


CONTROL LANGUAGE USAGE 
PARAGRAPH PARAMETER 


COMMENTS 
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BMDPAR - DERIVATIVE-FREE 110*1 LINEAR REGRESSION 


/PROB 

/INPUT 

/VARIAB 

/TRANSF V si 
/SAVE 
/END 
/PRINT 
(Optional) 


/PLOT 

(Optional) 

/REGRES 

(Required) 


SHAPE=name. 
TYPE=name. 
WEIGHT=v. 

le MAIN CHART 
CONTENT=list. 

C ASE = #. 

MATR=list. 

XVAR=list. 
YVARrlist. 
NORMAL. 
DEPEND=name. 

INDEP=list. 

METH=name. 


PEN=#. 
NUM=#. 
TOL = //. 
ZERO. 


SQUARE, LOWER. (SQUARE) 

One of: DATA, COVA, CORR. (DATA) 
Variable to weight eases (none/ 
prev.) 

One or both of: DATA, COVA.(DATA) 

Number of eases data printed.(5/ 
prev.) 

One or more of CORR, COVA, RESI, 
RREG. (CORR/prev.) 

X-axis variables, (none) 

Y-axis variables, (none) 

Normal plots, (no/prev.) 

Dependent variables, at least 3. 

(none) 

Independent variables, at least 
3- (none) 

One of: RSQ R-squared, ADJ adj. 
R-squared, CP Mallows Cp, NONE 
ordinary regression on all sub¬ 
sets. (CP/prev.) 

Number to add to Cp as penalty 
for including a variable.(2/prev.) 
Number of 'best' subsets to be 
chosen. (10 max.) (5/prev.) 
Tolerance for matrix inversion. 

(.0001/prev.) 

Request zero intercept.(no/prev.) 


BMDPAR DERIVATIVE-FREE NON LINEAR REGRESSION 

Estimates the parameters of a non linear function by least 
squares using a pseudo-Gauss-Newton algorithm (Ralston and 
Jennrich,1977). PAR is appropriate for a wide variety of 
functions that are not linear in the parameters, and for 
which derivatives are difficult to specify or costly to 
compute. The non linear function must be specified by 
FORTRAN statements input into the code of the program 
prior to running the program. The derivatives are not 
specified. 

OUTPUT 

- Weighted means and std. devs. for each variable. 

Residual sum of squares at each iteration. 

The minimal residual sum of squares and estimates of the 
parameters after the final iteration. 

Final asymptotic correlation matrix of the parameter esti¬ 
mates. 

- Estimated final mean square error. 

The asymptotic std. devs. of the parameter estimates 
(zeros indicate estimation problems.) 

A table containing, for each case; The residual y-f, the 
observed y, the predicted value, the standard error of the 
predicted value and the values of the independent 
variables and the case WEIGHT variable. 

SPECIFYING THE REGRESSION FUNCTION IN FORTRAN 

Any function can be fitted to the data by specifying it in 
the FORTRAN subroutine FUN. The specific FORTRAN 
subroutine FUN containing the regression function de¬ 
sired a problem must then become a part of the program 
BMDPAR for that run. The following is an example sub¬ 
routine FUN for the regression equation. 
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f=p(1)*EXP(p(2)*t) + p(3)*EXP(p(l|)*t) : 


SUBROUTINE FUN(F,PX,N,KASE,NVAR,NPAR IPASS xmnni 
IMPLICIT DOUBLE PRECISION (A-H O-Z) ’ ’ > 

IMPLICIT INTEGER*!) (I-N) 

DIMENSION P(NPAR),X(NVAR) 

RETURN 5 * DEXP (2) *X (3 5 )+P(3)*DEXP(P(A )*X (3}) 

END 


CONTROL LANGUAGE USAGE 

PARAGRAPH PARAMETER COMMENTS 


/PROB 
/INPUT 
/VARIAB 
/TRANSF 
/SAVE 
/END 
/REGR 
(Required) 


/PARAM 

(Required) 


/PLOT 

(Optional) 


see MAIN CHART 

NUMBz//. 

DEPENDzv. 

PARAM = //. 

TITLEz'c'. 
NUMBz#. 

WEIGHTzv. 

CONSTz//. 

PRINTzList. 

CONVz#. 
ITERz#. 

HALV = It . 

TOL = #. 

MEANSQzff. 

PASSz1.or2. 

LOSS. 

INIT=#list. 

MAXz/flist. 

MIN=#list. 

NAMEzList. 

DELTA=#list. 


FIXEDznames. 
or subscripts 
CONSTz#list. 


LIMITzfflist. 
VARzList. 


RESID. 

NORM. 

DNORM. 

SIZE = //1 , tf2. 


Number of analyses for which pre¬ 
dicted values are saved in BMDP 
Dependent variable. Required 
(none/prev.) 

# of parameters. Required.(none/ 
prev.) 

riai e »°5 ana ^ sis « <160 ools.(blank) 
Code // to select among several' re¬ 
gression functions, (none/prev.) 
Variable containing case weights. 

(no variable/prev.) 

it of inequality constraints. (none/ 

prev.) 

Data printed for each case, (all 
variables/prev.) 

Convergence criteria.(.00001/prev.) 
Max. # of iterations.(50/prev.) 

Max. It of increment halvings. 
(5/prev.) 

Check against round-off errors. 

(10 /prev.) 

If specified, used to compute 
asymptotic std. devs. 

? of passes of data for each 
iteration. (1/prev.) 

Loss function specified in FORTRAN 
subroutine, (no/prev.) 

Starting values for parameters. 
Required, (none/prev.) 

Upper limit of parameters.(none/ 
prev.) 

Lower limits of parameters.(none/ 
prev.) 

/prev S ) f ° r parameters -( p (subscript) 

Step size for computing additional 
starting values of the function. 
tntt^J ibscded, (.1* init ()) when 
INIT ( 1 )= o, .01 otherwise/prev.) 
Parameters whose estimates are 
). lx ??. st ioitial values . (none) 
Coefficients of inequality con¬ 
straint, may be repeated.(o 0 
/prev.) ’ 

Upper A lower limits for CONST 
may be repeated. (- , /prev.) 

Variables against which predicted 
and observed values and residuals 
are plotted, (none/prev.) 

Residuals and residuals squared 
plotted against predicted values. 

(no/prev.) 

Print normal probability plot of 
residuals.(no/prev.) 

Print detrended normal prob. plot 
(no/prev.) H 

Width A height of plot. (50,50/ 
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BMDP1S - MULTIPASS TRANSFORMATIONS 


prev.) 


BMDP1S MULTIPASS TRANSFORMATION 

Can be used when information from the data is required to 
compute transformations. 

OUTPUT 

For each pass, the transformed data and any statistics 
computed. The transformed data can be saved on a save 
file. 


CONTROL LANGUAGE USAGE 


PARAGRAPH PARAMETER COMMENTS 


/PROB ] 

/INPUT L see MAIN CHART 
/VARIAB f 

/SAVE PASS=#. 


/END J 

/PASS TITLE= 1 text' 

(Required) 

/COMPUTE MEAN. 


HARMONIC. 

GEOMETRI. 

EXTR. 

VARIABrlist. 

/PRINT MEAN. 

(Optional) 

HARMONIC. 

GEOMETRI. 

EXTR. 

DATA. 

VARIAB=List. 


Sequence t of PASS in which to 
create BMDP save file, required. 

(none/prev.) 

Alphanumeric pass title <160 cols, 
(blank) 

Compute mean and std.dev.(yes/ 
prev.) 

Compute harmonic means.(no/prev.) 
Compute geometric means.(no.prev.) 
Compute max. 4 mins, (no/prev.) 
Variables used in computations, 
(all variables) 

Print means and std. devs.(yes/ 
prev .) 

Print harmonic means.(no/prev .) 
Print geometric means.(no/prev .) 
Print max. 4 mins, (no/prev.) 

Print transformed data.(no/prev.) 
Variables involved.(same as COM¬ 
PUTE) 


BMDP3S NONPARAMETRIC STATISTICS 

Performs one of more of the following nonparametric sta¬ 
tistical tests: 


Sign test. 

Wilcoxon signed-ranks test. 
Kendall rank correlation coeff. 
Spearman rank correlation coef. 
Friedman two-way ANOVA. 

Kendall coef. of concordance. 
Mann-Whitney U test. 
Kruskal-Wallis one-way ANOVA. 


OUTPUT 


Test statistics requested. 

Level of significance or probability of occurrence for 
each test statistic except Kendall and Spearman rank 
correlation coef. 

CONTROL LANGUAGE USAGE 

PARAGRAPH PARAMETER COMMENTS 


/PROB 

/INPUT 

/VARIAB 

/GROUP 
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/TRANSF ^ 
/SAVE 

/END J 

see MAIN CHART 


/TEST 

TITLE='text' 

Test title <160 cols.(blank) 

(Required) 

VARIAB=list. 

Variables to be tested.(all in 
USE/prev.) 


SIGN. 

Sign test, (no/prev.) 


WILCOX. 

Wilcoxon signed-ranks.(no/prev. 


KENDALL. 

Kendall rank correlation coef. 
(no/prev.) 


SPEAR. 

Spearman's rank corr. coef.(no/ 
prev.) 


FRIED. 

Friedman two-way ANOVA.(no/prev 


KRUWAL. 

Kruskal-Wallis one-way ANOVA. 
(no/prev.) 


BMDP1V ONE-WAY ANALYSIS OF VARIANCE AND COVARIANCE 

Performs a one-way ANOVA, or a one-way analysis of covar¬ 
iance if covariates are specified. Group sizes may be un¬ 
equal. Parallel analysis can be performed using several 
dependent variables. 

OUTPUT 

For ANOVA: 

Group means, ANOVA table and/or pairwise t-test for group 
means and t-tests for contrasts of group means. 

For analysis of covariance on each dependent variable: 

- Regression coefs., std. errors and t values. 

- Group means, adj. group means, and std. errors. 

ANOVA table with test of slope. 

Pairwise t-test for adjusted group means and t-tests for 
contrasts. 

- Regression ooef. for each covariate in each group. 

Optional: 

- Correlation matrix for the reg. coefs. and/or adj. group 
means. 

Scatter plots. 

- Group statistics and variance, covariance and/or corre¬ 
lation matricies. 

CONTROL LANGUAGE USAGE 

PARAGRAPH PARAMETER COMMENTS 


/PROB 

/INPUT 

/VARIAB 

/GROUP 

/TRANSF 

/SAVE 

/END 

/PRINT 

(Optional) 


see MAIN CHART 


MAX. 

MIN. 


MEAN. 

COVA. 

CORR. 


Print max. value of vars. by 
groups, (no/prev.) 

Print min. value of vars. by 
groups, (no/prev.) 

Print mean value in groups.(no/ 
prev.) 

Print covariance matrix.(no/prev) 
Print correlation matrix, (no/ 
prev.) 
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BHDP2V - REPEATED MEASURES ANOVA « A'lCDVA 



TOT. 

Print total covariance matrix.(no/ 
prev.) 


BET. 

Print between cov. matrix.(no/prev) 


WIT. 

Print pooled within groups cov. 
(matrix, (no/prev.)) 


RREG. 

Print corr. between coefs. of co¬ 
variates. (no/prev.) 


MCOR. 

Print corr. matrix for adj. group 
means, (no/prev.) 


TITLE=’tex t'. 

Subproblem title <160 cols.(blank) 

/DESIGN 

DEPEND=list 

Dependent variables, (all except 

(Required, 


GROUPING variable/prev.) 

may be 

INDEP=list. 

Independent variables, (none) 

repeated) 

GROUP=list. 

Groups in subproblem.(all/prev.) 


CONTR=list. 

Coefficients for linear contrasts. 


PLOT. 

Plot each dependent variable 
against the independent variables, 
residuals, squared residuals,and 
predicted, (no/prev.) 


BMDP2V ANALYSIS OF VARIANCE AND COVARIANCE, INCLUDING REPEATED 
MEASURES 

Performs ANOVA or analysis of covariance for a general 
repeated measures model. For each subject, the trial 
factors must have a complete factorial structure with no 
missing observations. The group factors can be crossed 
or nested; group indices are read as data and determine 
the nesting relationships without further specification. 

Unequal cell sizes are allowed but not empty cells. 

OUTPUT 

- ANOVA table. 

- Cell means and std. devs. for dep. var. and covariates. 
Adjusted cell means and reg. coef. 

Design information. 

CONTROL LANGUAGE USAGE 

PARAGRAPH PARAMETER COMMENTS 


/PROB 
/INPUT 
/VARIAB 
/GROUP 
/TRANSF 
/SAVE 
/END 
/DESIGN 
(Required 


see MAIN CHART 


RESI. Save residuals and predicted. 

(yes/prev.) 

Alphanumeric design.(none/prev.) 
Subscripts of grouping variables, 
(none/prev.) 

Subscripts of dependent variables, 
(none/prev.) 

Subscripts of dependent variables. 
Subscripts of covariates.(none/ 
prev.) 

Names of trial factors. Up to R 
characters.(R,S,T,V,W,Y,Z/prev.) 
ORTH. Requests orthogonal decomposition 

of the trial main effects.(NO/ 
prev .) 

POINT (It) = #list. Spacing for levels of trial factor 


FORM='o'. 
GROUP=list. 


LEVEL=list. 


DEPEND=list. 
COVA=list. 


NAME=list. 
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SYMM. 

INCL=#list. 


RESI=name. 


PRINT. 
TOL = //. 

EXAMPLES OF DESIGN: 
LONG FORM 


#. may be repeated. (1,2,3,...) 
Test for compound symmetry.(no/ 
prev.) 

Terms to be included or excluded - 
from the model for the grouping 
factors, (include all main effects 
and iterations for grouping fac¬ 
tors) 

Method of computing residuals and 
predicted values MEAN or ORTH. 
(MEAN/prev.) 

Residuals and predicted values 
printed.(no/prev .) 

Tolerance limit. (.01/prev.) 


SHORT FORM 


GROUPzI. DEPEND=2,3,A. LEVEL=3. 
GROUP=1,2.DEPEND=3 TO 8.LEVEL=2,3. 
GROUP= 1. DEPEND = 3,5,7. C0VA=2 1) 6 
LEVEL = 3. 

GROUPs1. DEPENDs3,4,6,7,9,10. 

LEVELs 3,2. C0VAs2,2,5,5, 


FORMs 'G ,3(Y)'. 
FORMS 1 2G,2(3(Y))'. 
FORMs'G,3(X,Y)'. 

FORMs'G,3(X,2(Y))'. 
8 . 


BMDP3V GENERAL MIXED MODEL ANALYSIS OF VARIANCE 

This program uses the maximum likelihood (ML) and re- 
stricted maximum likelihood (REML) approaches to the fixed 
°° effi ° i ? nts model. The program is designed to 
handle mixed models of quite arbitrary form without 
requiring the balance demanded by P2V. On the other hand, 
the maximum likelihood approach is not as well studied or 
as widely available as the more classical approach of P2V 
Therefore, it is recommended that this program be used in 
conjunction with more classical approaches and for mixed 
Son Eunni,7 S not satlsf y 1 hg the balance requirements of 

rcV OP bMUUoV. 


Reports the estimates and std. devs. of the fixed effects 
parameters, variance components, estimates of the cell 
means, and residuals. 


The log-likelihoods for the model and for sub models (that 
are specified) are computed. 


CONTROL LANGUAGE USAGE 


PARAGRAPH PARAMETER 


/PROB j 

/INPUT N. see MAIN CHART 

/VARIAB WEIGHTrV. 


/GROUP 
/TRANSF 
/SAVE 
/END 
/DESIGN 
(Required) 


DEPEND=v. 

COVArList. 
FIXED=List. 


FNAME=List. 
RANDsList. 


COMMENTS 


Variable containing weights.(none/ 
prev.) 


Dependent variable, required.(none/ 
prev.) 

Covariates.(none/prev.) 

Variables defining fixed effects 
components, may be repeated.(none/ 
prev.) 

Names of fixed components.(FIX 
(subscript)/prev.) 

Variables that define a random 
effects component, may be repeated 
(none/prev.) 
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BMDP1L - LIFE TABLE ANALYSIS 

METH=name. 
PARAM = //list. 

MAXIT=#. 

/HYPOTH FIX=#list. 

(Optional) 

RAND=#list. 

/GROUP 


One only of ML or REML. 

Initial estimates for parameters, 
(all 0.0 except =1/prev.) 

Max. It of iterations. (30/prev .) 
Fixed components set to zero. 

(none) 

Random components set to zero. 

(none) 

Required for fixed effects. See 
MAIN CHART. 


BMDP1L LIFE TABLES AND SURVIVAL FUNCTIONS 

Estimates the survival (time-to-response) distribution of 
patients who have been observed over varying periods of 
time. The data can arise from experiments in which all 
subjects enter the study either concurrently or at differ¬ 
ent times. In the latter case, all patients are observed 
for unequal amounts of time. A major advantage of the 
techniques used in this program is that all patient fol¬ 
low-up information (although unequal in time) is used in 
estimating the survival probabilities. 

OUTPUT 

- Estimates of the survival function using the actuarial- 
type life table(Cutler-Ederer). 

Estimates of the survival function based on the product 
limit (Kaplan-Meier) estimate. 

- The Mantel-Cox and Breslow statistics to test the equality 

of survival distributions for different patient groups 
(these statistics are the analogues of non parametric rank 
statistics that are appropriate to censored data.) 

CONTROL LANGUAGE USAGE 

PARAGRAPH PARAMETER COMMENTS 


/PROB 
/INPUT 
/VARIAB 
/GROUP 
/TRANSF 
/SAVE 
/END 
/EST 

(Optional) 
may be 

repeated) one PERI0D=(). 

only WIDTHr/). 

CUTP=01ist. 

VAR=List. 


see MAIN CHART 


METH=n 

PRINT. 


PLOT=list. 
SIZE=#1 ,112. 
GR0UP=v. 


STAT=list. 


/FORM UNIT=name. 

TIME=V. 
STAT=V. 
RESP=#list. 


LOSS=#list. 


Required if GROUPing var. is 
specified and it has more than 
10 different values. 

One only of LIFE,PROD.(LIFE/prev.) 
Print survival distribution esti¬ 
mate . 

Provide the intervals for the life 
table when METH=LIFE. (PERIOD=10 
prev.) 

Variables printed with product- 
limit estimate, (none) 

One or more of SURV,LOG,CUM,HAZ, 
DEN. (SURV/prev.) 

Width and height of plots.(100, 
50/prev.) 

Variable used to classify cases 
into groups.(none) 

One or both of MANTEL, BRESLOW. 

!none/prev.) 

One only of DAY, WEEK, MONTH, 

YEAR. (MONTH) 

Variable that contains survival 
times in UNITS, (none-) 

Variable containing status of 
patient, (none) 

Values of STATus variable that 
are responses.(first code of 
STATus) 

Values of STATus variable that 
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ENTRY=List. 


are losses.(none) 

Month, day and year of patient 
entry, (none) 

TERM=List. Month, day and year of patient 

termination, (none) 

The following are used when a life table is used as input. 


NEN = v. 


NDEAD=V. 

NLOST=V. 

NWITH=V. 

INTER=V. 


Variable containing II of parients 
entering interval.(sum of NDEAD + 
NLOST + NWITH) 

Variable containing II of patients 
dying in interval.(none) 

Variable containing II lost in 
interval, (none) 

Variable containing II withdrawn 
during interval.(none) 

Variable containing lower limits 
for intervals. (0,1,2,...) 


BMDP EXAMPLES===============================: 

EXAMPLE 1: 


/PROBLEM TITLE IS 'VARIABLE DEFINITIONS AND SAVE FILE CREATION 
EXAMPLE USING BMDPID 1 . 

/INPUT VARIABLES ARE 16. 

CASES ARE 50. 

FORMAT IS '(FA.0, 2FN.1, 2FA.0, FA.1, 3FA.0, 5FN.0, AN)'. 
/VARIABLE LABEL IS ID. 

NAMES ARE HREATIN, CALCIUM, PHOSPHRS, CLUCOSE, BUN, 
URICACID, CHOLSTRL, PROTEIN, ALBUMIN, 

BILIRUBNM, ALKAPHOS, LACTDEHY, SGOT, BRTHPILL, 
ILLNESS, ID. 

MIN IS (7) 100.0. MAX IS (7) 500.0. 

/CATEGORY OUTPOINTS (1) ARE 1.5, 2.5. 

/PRINT MIN. MAX. 

/SAVE NEW. UNIT IS 9. CONTENT = DATA. CODE = WERNER. 

LABEL IS 'WERNER DATA AS OF 3-13-7M'. 


/END. 


50 data records 

/FINISH. 


EXAMPLE II: 

/PROBLEM TITLE IS 'EXAMPLE T-TEST ON SELECTED GROUPS USING 
BMDP3D'. 

'/INPUT VARIAB ARE 2. CASES ARE 13. FORM IS ’(2F2.0)'. 
/VARIAB NAMES ARE CATEGORY, VALUE. GROUP = CATEGORY. 
/GROUP CODES = 1 TO N. NAMES ARE A, B, C, D. 

/TEST GROUPS ARE A, B, D. 

/END. 

. 13 data records 

/FINISH. 

EXAMPLE III: 


1C 



BMDP - EXAMPLES 


/PROBLEM TITLE IS 'DISCRIMINANT ANALYSIS ON THREE GROUPS USING 
BMDP7M'. 

/INPUT VARIAB = 5. CASES = 150. FORM = '(AF3.1, F3.0)'. 
GROUPS ARE 3. 

/VARIABLE NAMES ARE 'SEPAL L', 'SEPAL W, 'PETAL L', 'PETAL W, 
IRIS. 

GROUPING IS IRIS. 

/GROUP CODES ARE 1 TO 3. 

NAMES ARE SETOSA, VERSICOL, VIGINIC. 

/PRINT CLASSIFICATION IS 1 TO A. 

/END. 

. 150 records 

/FINISH. 
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EXAMPLE IV: 

/PROBLEM TITLE 'FREQUENCY COUNTS USING BMDP2D ON CARD DATA IN 
SEPARATE FILE'. 


/INPUT UNIT IS 7. 

FORMAT IS '(FA.0, 2FA.1)'. 

/VARIABLES NAMES ARE HREATING, CALCIUM, 

/END. 

/FINISH. 



PHOSPHRS. 


EXAMPLE V: 

/PROBLEM TITLE 'STEPWISE REGRESSION BMDP2R USING COVARIANCE 
MATRIX ON A SAVE FILE.'. 


/INPUT UNIT IS 9. 

CONTENT IS COVA. 

CODE IS WERNER. 

LABEL IS 'WERNER COVARIANCE MATRIX FROM BMDP8D'. 

/PRINT COVARIANCE. 

/REGRES TITLE IS 'PREDICT BUN'. 

DEPEND IS BUN. 


/END. 

/FINISH. 




* 
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ACCEPT (SIR), 27 
ACCEPT REC IF (SIR), 12 
across record functions (SIR),5 
ADD (SIR), 27 
ADD REC (SIR), 18 
AFTER REPORT (SIR), 37 
AFTER RETRIEVAL (SIR), 21 
AGGREGATE (SPSS), 49 
algebraic expressions (SIR), 12 
altering a schema (SIR), 17 
AMOD (SIR), 4 
ANOVA (SPSS), 50 
ANOVA-mixed model (BMDP),103 
-one way (BMDP) , 101 
-repeated measures (BMDP),102 
-SPSS, 50 

arithmetic expressions (SIR),12 
arithmetic operators (SIR),3 
ASSIGN MISSING (SIR), 21 
assigning password (SIR), 9 
assigning security (SIR , 9 
AT END (SIR), 37 


batch retrieval commands (SIR),20 
batch retrieval example (SIR),42 
BEFORE REPORT (SIR),38 
bivariate plotting (BMDP),78 
BLANK LINES,(SIR), 21 
block clustering (BMDP),86 
blocks in retrievals(SIR) , 20 
BMDP SAVE FILE (SIR),34 
BMDP-computation (BMDP), 74 
-control 1anquaqe (BMDP), 68 
-examples (BMDP), 105 
-introduction (BMDP), 68 
-main chart (BMDP), 70 
-reserved words (BMDP), 69 
-statistical programs (BMDP),75 
-transformatlonsfBMDP), 74 
BMDP1D-simple data desc.(BMDP) ,75 
BMDP1F-two-way continqency{BMDP),81 
BMDPIL-life tables (BMDP), 104 
BMDP1M-clus ter analysis(BMDP) ,85 
BMDP1R-linear regression(BMDP),91 
BMDP!S-multipass trans. (BMDP),100 
BMDP1V-one way ANOVA (BMDP),101 
BMDP2D-frequency counts(BMDP),76 
BMDP2F-two-way with indep.{BMDP) ,83 
BMDP2M-C1us ter analysis(BMDP),86 
BMDP2R-stepwise m.l. reg.(BMDP),92 
BMDP2V-repea ted meas .ANOVA(BMDP),102 
BMDP3D-t-tests (BMDP),76 
BMDP3F-multiway cont.tables(BMDP),84 
BMDP3M-block cl ustering(BMDP),86 
BMDP3R-nonlinear regression(BMDP),94 
BMDP3S-nonparametric stats (BMDP),101 
BMDP3V-mixed model ANOVA (BMDP),103 
BMDP4D-alphanumeric counts(BMDP) ,77 
BMDP4M-factor analysis (BMDP),87 
BMDP4R- pri nc i pie comp. reg. (BMDP),95 
BMDP5D-univariate plots (SPSS),78 
BMDP5R-polynomial reg. (BMDP),96 
BMDP6D-bivariate plots (BMDP),78 
BMDP6M-canonical corr. (BMDP) ,88 
BMDP6R-multivariate reg.(BMDP),96 
BMDP7D-histograms & ANOVA (BMDP),79 
BMDP7M-dis c riminant analysis (BMDP), 

BMDP8D-correlation (BMDP), 80 
BMDP9D-multidimensional data (BMDP), 

BMDP9R-a 11 subsets reg.(BMDP),97 
BMDPAM-missing data est. (BMDP),90 


BMDPAR-non linear reg.(BMDP),9 
branching (SIR), 33 
BREAKDOWN (SPSS), 50 

CALL (SIR), 32 
CANCORR (SPSS), 51 
canonical correlation (BMDP),8 
case accessing (SIR), 22 
case definition (SIR), 9 
case deletion (SIR) , 32 
CASE ID (SIR), 9 
CASE IS, (SIR), 32 
case selection (SIR), 32 
case structure listing (SIR), 
CDATE (SIR), 3 
changing a schema (SIR), 17. 
character set (SIR) , 1 
CHECK (SIR), 178 
CIR (SIR), 10 
CLEAR (SIR), 27 
CNT (SIR), 4 
CNTR (SIR), 5 
command format (SIR), 2 
COMMENT (SIR), 6 
comments in SIR (SIR), 6 
common information record (SIR' 
10 

COMMON LIST (SIR) , 10 
computation (BMDP), 74 
COMPUTE (SIR), 12, 22 
CONDESCRIPTIVE (SPSS), 52 
CONDESCRIPTIVE (SIR), 34 
contingency table analysis 
(BMDP), 82 

continuation lines (SIR), 2 
control 1anguaqe (BMDP), 68 
control words (SPSS), 48 
COPY (SIR), 28 
COUNT (SIR), 6 
CROSSTABS (SPSS), 52 
CTIME (SIR), 3 
CWRITE (SIR), 28 

DATA LIST (SIR), 12 

data types (SIR), 2 

DATE (SIR), 37 

date functions (SIR), 3 

DATE VAR LIST (SIR), 12 

DATET (SIR), 3 

DELETE (SIR), 28 

DELETE CASE (SIR), 19,32 

DELETE CASE IF (SIR) , 19 

DELETE REC (SIR), 19 

DELETE REC IF (SIR), 19 

DELETE RECORD (SIR), 32 

delimiters (SIR), 1 

designing retrievals (SIR),20 

DETAIL (SIR), 38 

detailed schema listing (51R), 1 

DISCRIMINANT (SPSS), 53 

DISK RESTORE (SIR), 9 

DO REPEAT (SIR), 8 

DOCUMENT (SIR), 6 

dumping SIR records (SIR),41 

dumping to cards (SIR), 36 

dumping to tape (SIR), 39 

EDIT (SIR), 8 
editor return (SIR), 30 
editor usaqe (SIR), 26 
EDITOR $ (SIR), 31 
-ACCEPT (SIR), 27 
-ADD (SIR), 27 
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EDITOR (SIR) cont. 

-CALL (SIR), 32 
-CLEAR (SIR), 37 
-COPY (SIR), 28 
-CWRITE (SIR), 28 
-DELETE (SIR), 28 
-END (SIR), 28 
-EP (SIR), 31 
-FINISH (SIR), 28 
-GET (SIR), 28 
-HELP (SIR , 28 
-LIST (SIR), 28 
-OUTPUT (SIR), 29 
-PAGE (SIR), 29 
-PREAD (SIR), 29 
-PURGE (SIR), 29 
-PWRITE (SIR), 29 
-qUIT (SIR), 28 
-READ (SIR), 29 
-REPLACE (SIR), 29 
-RESEQUENCE (SIR), 30 
-REWIND (SIR), 30 
-RUN (SIR), 30 
-SAVE (SIR), 30 
-TAB (SIR), 31 
-WRITE (SIR), 31 
-X (SIR), 31 
END (SIR), 22,28,38 
END LEVEL (SIR), 38 
END REPEAT (SIR), 8 
END REPORT (SIR), 38 
EP (SIR), 31 
ERROR LIMIT (SIR), 8 
EVICT REC (SIR), 20 
examples-BMDP {BMDP).105 
-REPORT (SIR), 43 
-SIR batch retrieval (SIR),42 
-SIR schema defn.(SIR),41 
-SPSS (SPSS), 66 
execution options (SIR), 45 
execution parameters (SIR), 27 
EXISTS (SIR), 6 
EXIT (SIR), 32,38 
EXIT BLOCK IF (SIR), 22 
EXIT IF (SIR), 38 

FACTOR (SPSS), 54 

factor analysis (BMDP), 87 

FILE NAME (SIR), 8 

file statistics (SIR), 16 

FINISH (SIR), 8,28 

fixed decimals (SIR), 14 

fixed format (SIR), 2 

FOOTING (SIR), 22 

FOR EACH CASE (SIR), 22 

FOR EACH REC (SIR), 22 

formats (SIR), 13 

formatting (SIR), 25 

formatting commands (SIR),22 

FREQUENCIES' (SPSS), 44 

FREQUENCIES (SIR), 34 

FST (SIR), 4 

FSTR (SIR), 5 

functions (SIR), 3 

G3SLS (SPSS), 55 
general rules (SIR), 1 
GET (SIR), 28 
GET FILE (SIR), 8 
GUTTMAN SCALE (SPSS) ,55 


INDEX (cont.) 

HEADING (SIR), 21 ,23 
HELP (SIR), 28 
histograms (SIR), 34 
histograms (SPSS), 54 
histograms (BMDP), 78 

IF (SIR), 12, 23, 32 
index key size (SIR) , 10 
initialization (SIR), 2 
initializing variables (SIR) 23 
INPUT FORMAT (SIR), 13 ’’ 

INPUT MEDIUM (SIR) , 7 
input record size (SIR), 10 
integer constants (SIR), 2 
INTEGER LIST (SIR) , 13 
integer variables (SIR), 13 

ln /?rnf tive ret Fieval example 
(SIR), 44 

interactive usage (SIR), 26 
introduction (BMDP), 68 

JFACTOR (SPSS), 56 
JULC (SIR), 4 
JULN (SIR), 4 
JUMP (SIR), 33 
JUMP TO LEVEL, 38 

key size (SIR), 10 

labels listing (SIR), 16 

labels for variables (SIR) 25 

LEN (SIR), 4 

LEVEL (SIR), 38 

LIFE TABLES (BMDP), 104 

linear regression (SPSS), 61 

UsT(Si r Rr28 i0n (BHDPj * 91 

LIST SCHEMA (SIR), 16 
LIST STATS (SIR), 16 
listing SIR files (SIR), 41 
locked schema (SIR) , 17 

logical computation (SIR ),12 
logical operators (SIR), 3 
LSTR (SIR) ,5 

main chart (BMDP) , 70 
MANOVA (SPSS), 57 
MAX (SIR), 4 

MAX INPUT COLS (SIR), 10 
MAX KEY SIZE (SIR),10 
MAX REC COUNT (SIR),10 13 
MAX REC TYPES (SIR), in 
MAXR (SIR), 5 
MEAN (SIR), 4 
MEANR (SIR), 5 

K'l&l!" ("«• 

MI NR (SIR), 5 

MISSING VALUES (SIR), 13 

MISSING VALUES-retrievals (SIR), 

mixed model ANOVA (BMDP). 103 
MODIFY CASE (SIR), 19 
MODIFY CASE IF (SIR), 19 
MODIFY REC (SIR), 19 
MODIFY REC IF (SIR), 19 
MODIFY SCHEMA (SIR), 19 
modifying a schema (SIR),17 


108 




INDEX (cont.) 


MOVE (SIR), 28 

MOVE VAR LIST (SIR), 23 

MULT RESPONSE (SPSS), 58 

mult.Iin.regression (SPSS),61 

mult.1in.regression (BMDP),91 

multi-pass trans. (BMDP), 100 

multi-way tables (BMDP), 84 

multiple summary records (SIR), 22 

N OF CASES (SIR), 11 

n-way tables (SPSS) ,52 

non linear regression (SPSS), 58 

non parametric tests (SPSS), 59 

NONLINEAR (SPSS), 58 

NONPAR CORR (SPSS), 59 

NPAR TESTS (SPSS), 59 

NREAD (SIR), 33 

NUMBERED (SIR), 9 

NUMBR (SIR), 4 


ON ERROR-report (SIR), 39 
ON PAGE -report (SIR), 39 
one way ANOVA S COVA (SPSS),50 
one way ANOVA & COVA (BMDP , 101 
ONEWAY ANOVA (SPSS), 60 
order of retrieval commands (SIR), 
21 

order of schema commands (SIR), 11 

OUTPUT (SIR), 29 

output specifications (SIR), 25 

PAGE (SIR), 29,37 
PAGE EJECT (SIR), 23 
PAGESIZE (SIR), 9 
paragraphs (BMDP), 69 
PARTIAL CORR (SPSS) , 61 
PASSWORD (SIR), 9 
PEARSON CORR (SPSS), 61 
PERFORM PROCS (SIR) , 22 
PLOT (SIR), 35 

plotting-bivariate (BMDP), 78 

plotting-univariate (BMDP), 78 

polynomial regression (BMDP), 96 

PREAD (SIR), 29 

PRESET (SIR), 23 

PRINT BACK (SIR), 7 

procedure calls (SIR), 32 

procedures (SIR), 34 

process rec (SIR), 23 

programming retrievals (SIR),20 

prompting (SIR), 33 

PURGE (SIR), 29 

PURGE SIR FILE (SIR), 40 

PWRITE (SIR), 29 


query commands (SIR), 26 
QUIT (SIR), 28 

RAND (SIR), 6 
random numbers (SIR), 6 
READ (SIR), 29 
READ INPUT DATA (SIR), 18 
READ SECURITY (SIR), 9 
real numbers (SIR) , 2 
REC SECURITY (SIR), 14 
RECODE (SIR), 24 
record accessing(SIR), 22 

record definition (SIR), 11 
record deletion (SIR), 32 
RECORD IS (SIR), 33 
record ordering (SIR), 14 


record processing (SIR), 23 
RECORD SCHEMA (SIR), 13 
record selection (SIR), 24 
record size (SIR) , 10 
RECS PER CASE (SIR), 11 
RECTYPE COLS (SIR), 11 
REGRESSION (SPSS), 61 
regular schema listing (SIR),16 
REJECT CASE IF (SIR) , 24 
REJECT REC IF (SIR), 12, 24 
relational operators (SIR),3 
RELIABILITY (SPSS), 62 
repeated measures ANOVA (BMDP), 
102 

repeating commands (SIR ),8 
repeating summary records (SIR), 
22 

REPLACE (SIR), 29 
REPLACE REC (SIR),18 
report generator (SIR), 36 
REPORT-AFTER REPORT (SIR),37 
-AT END (SIR), 37 
-BEFORE REPORT (SIR), 38 
-DETAIL (SIR), 38 
-END (SIR), 38 
-END LEVEL (SIR), 38 
-END REPORT (SIR), 38 
-EXAMPLE (SIR), 43 
-EXIT (SIR), 38 
-EXIT IF (SIR), 38 
-JUMP TO LEVEL (SIR), 38 
-LEVEL (SIR), 38 
-ON ERROR (SIR), 39 
-ON PAGE (SIR) , 39 
required commands (SIR), 6 
RESEQUENCE (SIR), 30 
RETRIEVAL (SIR), 24 
retrieval blocks (SIR), 20 
retrieving variables SIR),23 
RETURN (SIR), 30 
REWIND (SIR), 30 
reworking a schema (SIR), 17 
RUN (SIR), 30 
RUN NAME (SIR) , 9 

SAVE(SIR), 30 

save file -BMDP (SIR), 34 
-SPSS (SIR), 36 
scaled integers (SIR) , 2 
SCALED VAR LIST (SIR), 14 

SCATTERGRAM (SPSS), 63 
schema listing (SIR), 16 
schema case definition (SIR) ,9 
schema modification (SIR), 16 
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